Markus Ganser 


ON THE ELECTRO-CHEMO-MECHANICAL 
COUPLING IN SOLID STATE BATTERIES 
AND ITS IMPACT ON MORPHOLOGICAL 
INTERFACE STABILITY 


SCHRIFTENREIHE DES INSTITUTS 
FUR ANGEWANDTE MATERIALIEN BAND 90 


SICHT zusishing 


Markus Ganser 


On the Electro-Chemo-Mechanical 
Coupling in Solid State Batteries 
and its Impact on Morphological 
Interface Stability 


Schriftenreihe 
des Instituts fur Angewandte Materialien 
Band 90 


Karlsruher Institut fur Technologie (KIT) 
Institut fur Angewandte Materialien (IAM) 


Eine Ubersicht aller bisher in dieser Schriftenreihe erschienenen Bande 
finden Sie am Ende des Buches. 


On the Electro-Chemo-Mechanical 
Coupling in Solid State Batteries 
and its Impact on Morphological 
Interface Stability 


by 
Markus Ganser 


SICHT Fikisha 


Karlsruher Institut fur Technologie 
Institut fur Angewandte Materialien 


On the Electro-Chemo-Mechanical Coupling in Solid State Batteries 
and its Impact on Morphological Interface Stability 


Zur Erlangung des akademischen Grades eines Doktor-Ingenieurs 
von der KIT-Fakultät für Maschinenbau des Karlsruher Instituts für 
Technologie (KIT) genehmigte Dissertation 


von Markus Ganser, M. Sc. 


Tag der mündlichen Prüfung: 30. April 2020 
Referent: Prof. Dr.-Ing. Marc Kamlah 
Korreferent: Prof. Dr.-Ing. Thomas Böhlke 


Impressum 
NC Scientific 
Publishing 
Karlsruher Institut für Technologie (KIT) 
KIT Scientific Publishing 


StraBe am Forum 2 
D-76131 Karlsruhe 


KIT Scientific Publishing is a registered trademark 
of Karlsruhe Institute of Technology. 
Reprint using the book cover is not allowed. 


www.ksp.kit.edu 


© OO This document - excluding the cover, pictures and graphs - is licensed 
D4 under a Creative Commons Attribution-Share Alike 4.0 International License 
(CC BY-SA 4.0): https://creativecommons.org/licenses/by-sa/4.0/deed.en 


© OXO) The cover page is licensed under a Creative Commons 
Rare Attribution-No Derivatives 4.0 International License (CC BY-ND 4.0): 
https.//creativecommons.org/licenses/by-nd/4.0/deed.en 


Print on Demand 2021 - Gedruckt auf FSC-zertifiziertem Papier 


ISSN 2192-9963 
ISBN 978-3-7315-1047-5 
DOI 10.5445/KSP/1000123310 


On the Electro-Chemo-Mechanical Coupling 
in Solid State Batteries and its Impact on 
Morphological Interface Stability 


Zur Erlangung des akademischen Grades 
Doktor der Ingenieurwissenschaften 
der KIT-Fakultat fiir Maschinenbau 
Karlsruher Institut fiir Technologie (KIT) 


genehmigte 
Dissertation 


von 


Markus Ganser, M. Sc. 


Tag der miindlichen Priifung: 30. April 2020 
Referent: Prof. Dr.-Ing. Marc Kamlah 
Korreferent: Prof. Dr.-Ing. Thomas Böhlke 


Zusammenfassung 


Festkörperbatterien mit einer Lithium Metall Elektrode gelten in Bezug 
auf die erreichbaren Energiedichten als die nächste Generation in der 
Hochenergie-Batterietechnologie. Lithium Metall tendiert jedoch zu 
Dendritenwachstum, das ein Zellversagen zur Folge haben kann. Der 
Einsatz von Festelektrolyten ist ein vielversprechender Ansatz um die 
Lithium Metall Grenzfläche zu stabilisieren, die Mechanismen und die 
diesbezüglichen Materialanforderungen sind jedoch bisher weitgehend un- 
bekannt. Im Rahmen dieser Arbeit wird ein elektro-chemo-mechanisches 
Modell einer binären Festkörperbatterie hergeleitet, um die morphologis- 
che Grenzflächenstabilität von Lithium Metall zu analysieren. 

Das Modell beinhaltet eine Transporttheorie, welche die Elektrostatik, 
den Massentransport und die geometrisch nichtlineare Mechanik, darunter 
das Anschwellen des Elektrolyten beschreibt. Die Butler-Volmer Reak- 
tionskinetik wird derart verallgemeinert, sodass sie mechanische Span- 
nungen an der Grenzfläche berücksichtigt. Zudem wird eine Kinematik 
vorgestellt, die die Anlagerungs- und Ablöseprozesse als auch mechanische 
Deformationen während des Ladens und Entladens konsistent beschreibt. 

Basierend auf der vorgestellten Theorie werden Einflussgrößen auf 
die morphologische Stabilität von Polymerelektrolyten untersucht. Die 
Ergebnisse dieser Analyse zeigen sowohl einen positiven Zusammenhang 
zwischen den mechanischen Eigenschaften und der Stabilität, als auch, 
dass bessere Transporteigenschaften die Anforderungen an die Mechanik 
reduzieren. Die Resultate numerischer Studien sowie analytischen Be- 
trachtungen lassen zudem den Schluss zu, dass Dendritenwachstum ins- 


besondere durch inhomogene Grenzflächenwiderstände entsteht. 


Abstract 


Solid state batteries with a lithium metal electrode are considered the 
next generation of high energy battery technology. Unfortunately, lithium 
metal is prone to harmful protrusion or dendrite growth which causes 
dangerous cell failure. Although solid state electrolytes have shown a 
promising potential to stabilize the interface morphology, a clear under- 
standing of the related mechanisms and material requirements is still 
missing. Within this work the problem of protrusion growth is tackled by 
deriving a novel electro-chemo-mechanical theory tailored for binary solid 
state batteries which is then used to discuss the impact of mechanics on 
interface stability by extensive numerical studies. 

An ion transport theory in finite strain continuum mechanics is devel- 
oped that captures electrostatics, component transport and nonlinear 
mechanical interaction such as swelling, stress driven diffusion and defor- 
mation. An enhanced version of Butler-Volmer reaction kinetics enriched 
by the influence of mechanics is carefully derive. Furthermore, a model is 
proposed that consistently describes material deposition during charging 
and coupled deformation of metal electrode and solid electrolyte. 

The comprehensive theory then is used to analyze morphological stabil- 
ity for polymer electrolytes. Quasi static simulations reveal not only that 
increasing the mechanical stiffness favors interface stability but also that 
the stiffness required for stability decreases with increasing transport 
properties. The result of a transient numerical study and an analytical 
network model show that inhomogeneous interface conductivities are a 


likely reason for harmful protrusion growth. 
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1 Introduction 


1.1 Basics of Lithium Batteries 


Batteries are electrochemical storage devices that convert chemical energy 
into electric energy. Depending on the type of the electrochemical process 
and the materials used, a number of classes of battery storage devices 
are distinguished. Usually, one differentiates between primary batteries 
in which discharge is possible only once due to irreversible processes, 
and secondary batteries which can be operated reversibly in charge and 
discharge direction. The sate of the art technology in today’s secondary 
batteries is the lithium ion (Li-ion) technology, which dates back to 1970 
by Whittingham at Exxon [I]. The first reversible lithium battery was 
commercialized in the late 1980s by Moli Energy with a lithium metal 
anode, a MoS cathode material and a liquid electrolyte [2]. Although 
produced in millions, this early attempt of using a lithium metal anode 
caused tremendous safety problems, e.g. fire accidents caused by dendrite 
formation, and had to be recalled. In 1991, Sony was finally able to 
build a reliable Li-ion secondary battery cell by using ionic lithium with 
graphite as intercalation material and a liquid electrolyte. The basic 
concept of this electrochemical storage device has not been changed since 
then. The graphite material has become the quasi standard and the 
research focus was shifted to the counter-electrode cathode material. 
Recently, batteries with lithium metal anode in combination with solid 
state electrolytes obtained much attention and might be the successor of 
Li-ion batteries. These solid state batteries (SSB) feature metallic lithium 


and solely solid components within the electrochemical cell to increase 
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energy density and safety beyond the classical Li-ion technology with 
liquid electrolyte. Furthermore, SSB might pave the way to far-future 
technologies such as lithium-sulfur and lithium-air where also elemental 
lithium is utilized. 

Although this work is mainly motivated by the challenges of the lithium 
technology, we emphasize that many ideas can be transfered to alternative 


materials such as sodium based batteries. 


1.1.1 Fields of Application 


The availability of high energy secondary battery systems enables con- 
nectivity and electrification in daily life. Over the last decades, consumer 
electronics tremendously changed the way we communicate and work. 
Forecasts predict that the number of mobile phone users will rise to 4.8 
billion people worldwide by 2020 [3]. Thereof, the more sophisticated 
smartphones have a share of 60%. This success is closely related to 
durable and affordable batteries with high energy densities. Despite the 
large amount of devices, the capacity of a standard smartphone in 2018 
is quite small with roughly 10 Wh, which explains the low fraction in the 
extrapolated market share, see Fig. [1.1] The rising market of the Internet 
of Things [4] with several billions of small, connected and battery driven 
devices will further boost the need for high energy batteries. 

As can be seen in Fig. most of the predicted battery demand will 
arise from the mobility sector and especially the automobile market with 
a single vehicle requiring 5 to 15 kWh [5]. Traditional gasoline and diesel 
combustion engines exhaust greenhouse gases such as carbonate dioxide 
and nitrogen oxides that are known to accelerate climate change and 
harm the individual human health. Several countries therefore enforce the 
reduction of exhaust emission by subsidizing new environmental friendly 
technologies or even by banning cars with high exhaust output from 


the streets [6]. A transportation system with an electrified powertrain 
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Figure 1.1: The annual battery demand is predicted to increase by 32% per year from 
2016 to 2030 according to a recent study by McKinsey [7]. This is strongly driven by 
the mobility sector such as battery electric vehicles (BEV). 


for motorcycles, cars, buses, trucks and even flying vehicles is therefore 
a logical and necessary step into the future of mobility. Political and 
economic factors such as oil price development and infrastructure invest- 
ments further influence the future of mobility concepts. The requirements 
for electrochemical energy storage in transportation systems are high. 
The energy density, which is directly related to the driving range, has to 
reach critical values, especially because the users are accustomed to the 
high driving ranges of fossil fuels with their excellent energy densities, 
e.g. gasoline with 1700 Te [2]. For comparison, current Li-ion cells 
exhibit an energy density of up to 250 T and lithium metal batteries 
have the potential of the double density with up to 440 Ya Pl. Beyond 
that, it is mandatory that electrochemical storage devices provide a high 
level of safety. Hazards due to uncontrolled chemical reaction within 
the battery e.g. due to production faults [8] or accidents [9] have to be 
ruled out as far as possible. Other than mobile devices with a life span of 


approximately two years, batteries in cars also have to last at least one 
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decade before their capacity decreases below 80% [5]. Another important 
aspect for the acceptance of batteries in transportation devices is fast 
charging to enable reasonable waiting times for intermediate recharging 
stops. All of that has to be achieved for an affordable price, the major 
factor for the attractiveness of an electrified vehicle. Almost all major 
automotive manufactures provide or have announced fully or partially 
electrified vehicles, which can be distinguished in hybrid (HEV), plug-in 
hybrid (PHEV) and battery electric vehicles (BEV). The ultimate goal 
in the development of a BEV is to provide the same comfort and similar 
ranges compared to combustion driven vehicles to gain high costumer 
acceptance. 

Other applications for battery systems are stationary energy storages in 
combination with a smart grid infrastructure. Together with local energy 
generation such as photovoltaic or wind plants, they can potentially enable 
self-sustaining energy concepts for both single houses and communities 
and eliminate dependence on fossil fuels. As can be seen, secondary 
batteries play an increasingly important role in today’s innovation, which 


explain the intense research activities. 


1.1.2 Design and Working Principle 


The working principle is very similar for a conventional battery with a 
liquid electrolyte and a solid state battery with a solid electrolyte and 
a lithium metal electrode. Fig. shows the concept of both types 
including a description of the electrochemical and mechanical processes. 
We briefly describe the idea of a conventional battery before introducing 
the concept and advantages of a solid state battery with a lithium metal 
electrode. 

An ordinary battery has three functional layers, two porous electrodes 
and an electrically isolating porous separator in between. The three 


layers are infiltrated by a liquid electrolyte to enable an ion conducting 
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Figure 1.2: Conventional battery with a liquid electrolyte and graphite particles (grey) 
on the left and a solid state battery with a solid electrolyte and a lithium metal anode 
on the right. The dominant effects in a conventional battery are purely electrochemical 
whereas solid state batteries require the additional description of mechanics. 


pathway between the two electrodes. A standard liquid electrolyte in Li- 
ion batteries is based on hexafluorophosphate (LiPFg) as conducting salt 
and organic compounds, which are stable against the highly reactive ionic 
lithium. The porous electrodes are a composition of storage particles, 
binder and conductive carbon. The storage particles have both access to 
an electric pathway (with the conductive carbon) and to an ionic pathway 
(via the electrolyte). The storage particles determine the majority of the 
electrochemical response of the electrode and therefore, when speaking 
of electrode materials, one usually refers to the storage particles. The 
binder is necessary to obtain the mechanical integrity of the porous 
electrode. A commonly used combination of electrode materials is, for 
example, lithium graphite Li,Cg in the anode and lithium iron phosphate 
Li,_xFePO, in the cathode where x ranges between 0 and 1. The battery 
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is completely charged for x = 1 meaning that the majority of lithium is 
intercalated in the graphite. 

It can be seen that each Li-ion in a porous electrode is accompanied 
by several other compounds, be it the active storage material or the 
components to maintain electronic and ionic conductivity such that only 
a fraction of the electrodes contributes to energy storage. As we will see, 
one possibility to increase energy density is by substituting the porous 
anode with metallic lithium. 

The fundamental principle of energy storage of a secondary battery is 
hidden in the electrode materials. The electrodes serve as the storage 
of lithium. Depending on the chemical composition of the electrode 
material, the electrochemical energy of the ions varies heavily. For 
example, a lithium metal electrode has a higher energetic level compared 
to the Li¡..FePO, electrode. This means that Li-ions prefer to form 
a compound with Li,.,FePO, to being in the metal phase. If the two 
electrodes are connected by an electronic and ionic pathway and no 
external voltage is applied, the ions perform a spontaneous reaction 
from the higher energetic level to the lower one. This is a consequence 
of the second law of thermodynamics stating that reaction follows the 
direction which minimizes the total energy. The process can be reversed 
by applying an external voltage and thus charging the battery. The 


corresponding reaction equation reads 


Lit FePO, “E LiFePO, (1.1) 


charge 
which is a reduction-oxidation process and can be split into two steps. 


During discharge (charge) lithium is oxidized (reduced) at the metal 


electrode, 


discharge 
Li "=" Lit+e-. (1.2) 


charge 
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On the macroscopic level, this means that the amount of metal decreases 
(increases) layer by layer or, in other words, lithium is stripped from 
(deposited on) the metal electrode. In the second reaction, the Li-ion is 


reduced (oxidized) to form lithium iron phosphate 


Reir SEE LiFePO,. (1.3) 
charge 

The involved ions and electrons are transfered between the electrodes 
through the electrolyte and the external circuit, respectively. The process 
of oxidation is also known as anodic reaction and reduction is known 
as cathodic reaction. By definition, the electrodes are named during 
discharge and therefore the lithium metal is denoted as anode and the 
iron phosphate as cathode. Alternatively, one can utilize the reduction 
potential (defined against a hydrogen electrode) and refer to the lithium 
metal with its relatively small reduction potential as negative electrode 


and the iron phosphate as positive electrode. 


1.1.3 Dendrites and Solid State Electrolytes 


Enabling batteries with metal electrodes instead of an intercalation 
material might be one of the next milestones in increasing the energy 
density in electrochemical energy storage. These systems do not require 
a storage material at the negative electrode (e.g. graphite, see Fig. |1.2) 
and therefore offer best in class energy density due to the lowest possible 
anode potential in combination with the best theoretical gravimetric 
(3860 =) and volumetric capacity (2 2%) 21[T0]. However, pure lithium 
metal is highly reactive and has a low melting temperature of 180.6°C 
[L]. It is therefore rarely applied and only few characterizations have 
been executed to describe the complex mechanical response. The Young’s 
modulus of lithium is reported to be E = 1.9 — 8GPa [IIHT3]. It has a 
low mechanical strength (the reported strengths ranges from 1 MPa [12] 
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to 105 MPa [14]), shows visco-plastic behavior [15], is likely to creep and 
has an anisotropic micro-structure [14]. 

Unfortunately, a lithium metal electrode is highly unstable against 
traditional liquid-based electrolytes such as LiPFg. Heterogeneous nucle- 
ation, most likely due to a heterogeneous solid electrolyte interface (SEI), 
causes an uncontrolled dendrite growth in form of mossy and needle-like 
structures [I6HI8] that might cause the battery to fail by short-circuits. 
For example, Steiger et al. [19] reports filaments of a diameter of approx- 
imately 500nm and complex growing patterns. Several approaches have 
been presented to prevent dendrite growth in liquid electrolytes [TO], e.g. 
using micro/nanostructures [20], concentrated liquid electrolytes [21] as 
well as coatings and additives [22] to modify the SEI. 

An alternative approach is to substitute the liquid electrolyte and 
the porous separator with a solid electrolyte to enable metal anodes. 
The requirements on solid electrolyte materials in this context are a 
sufficient ionic conductivity of Li-ions, negligible electronic conductivity, 
a wide electrochemical stability window and mechanical stiffness which 
is said to suppress harmful dendrites growth. Besides addressing the 
issues of dendrites, solid electrolytes also have the potential to increase 
safety, specifically by decreasing flammability compared to liquid based 
systems [23]. A broad class of solid electrolytes, ranging from inorganic 
ceramics to polymers, are currently under investigation [24] [25]. Inorganic 
electrolytes, such as oxidic or sulfidic ceramics, belong to the class of single 
ion conductors. Usually, they offer a very high ionic conductivity with up 
to 1 S at a wide temperature range [25] and mechanical stiffnesses of up 
to 20 GPa [26]. However, depending on the composition, they have issues 
with stability against moisture, processing and cost [27]. Furthermore, 
in combination with lithium metal, they are prone to dendrite growth 
in the grain boundaries [28], suffer from adhesive problems [29] and are 


therefore not yet ready for the mass market. 
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[>] 
[polymer 
Bite JES lithium 
X 
100 um 
Figure 1.3: Protrusion growth in polymer electrolytes: (a) Photograph of a protrusion 


growth in a PEO polymer at high current density (taken from Brissot et al. ). (b) 
Three 3D renderings from a X-ray microtomopgraphy of a PS-PEO polymer (taken 
from Harry et al. [31]). 


This work focuses on the class of polymer electrolytes that act as a solid 
solvent for a salt such as LiTFSi. In these electrolytes, both anions and 
cations are mobile leading to the classification as a binary conducting 
electrolyte. The most prominent polymer electrolyte is poly(ethylene 
oxide) (PEO). It is thermodynamically stable against lithium (up to 
100°C [32]), less flammable than liquid electrolytes and is usually op- 
erated at elevated temperatures in order to exceed the glass transition 
temperature (60°C) to obtain sufficient ionic conductivity (0.1 5 [83]). 
The ion transport is in that case supported by segmentation motion of 
the PEO chains. However, at temperatures of e.g. 80°C, PEO is viscous 
and, hence, tends to lose its mechanical integrity. Nevertheless, it is 
partially able to suppress classical mosslike dendrites but suffers from 
lithium intrusion on the length scale of micrometers [30] [34] [35] as shown 
in Fig. [1.3p. 

Several polymer modifications have been suggested in order to increase 
electrochemical properties as well as mechanical strength [32] [36] 37], e.g. 
via adding ceramic nanoparticles [38], cross linking [39/41] or using gel 
electrolytes such as double polymer networks [42], zwitterionic copolymers 
[43] or ionic liquid-immobilized polymers [44]. Another approach are 
block/triblock copolymers [45H48] such as PS-PEO which follow the idea 
that a soft block (PEO) forms lithium conducting channels and a hard 
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block of polystyrene (PS) provides the mechanical strength. Nanochannels 
of ~10nm are able to conduct ions but block the growth of the 
beforementioned mosslike dendrites. Unfortunately, the microscopic 
structure of block copolymers reduces the transport properties since non- 
conducting material (PS) is added and also the segmentational motion 
of PEO is hindered [32] [49]. 

Although block copolymers decrease the likelihood and speed of harmful 
protrusion growth [31], the so-called globule and multi-globular structures 
remain a challenge (see Fig. [1.3p). Recent experimental studies 
correlate their appearance with sub surface structures [51] and impurities 
[50] in the lithium metal anode. However, the exact mechanisms causing 
the undesired protrusion growth and the measures necessary to suppress 
them are questions of ongoing research. 

Exploring undesired dendrite growth in experiments as well as synthe- 
sizing, processing and testing new kinds of solid electrolytes is complex, 
time-consuming and costly. Therefore, it is beneficial to derive mathemat- 
ical models that accurately describe relevant effects of electrochemistry, 
mechanics and their coupling in order to increase the understanding of 


the different phenomena in solid state batteries. 


1.2 State of the Art of Battery Modeling 


Different theoretical model description have been proposed to assess 
the complex response in battery systems. One can distinguish models 
for single ion conductors (e.g. [52]) and more sophisticated models for 
binary electrolytes. Binary electrolytes are usually described by the 
well-established concentrated solution theory [53], initially derived for 
liquid electrolytes and later applied to binary solid state electrolytes [54]. 
The ion transport of non-ideal electrolytes is parametrized by four mea- 


surable transport parameters, namely conductivity, transference number, 
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diffusivity and thermodynamic factor [55H57]. Ferrese et al. and 
subsequent works combine the conventional concentrated solution 
theory with a mechanical model to describe the in-plane redistribution 
of lithium metal during cycling. However, effects of stress driven dif- 
fusion, as e.g. known in active storage particles 67|, are not taken 
into account even though quite stiff materials are assumed. Bucci et al. 
[68] propose a multi-component transport model of charged species in a 
deformable host. Similar ideas are found e.g. in hydrogels, where small 
molecules diffuse in a network of long fibers and where a change of local 
component concentration leads to heterogeneous swelling or shrinkage of 
the host material and consequently to mechanical stresses, which can in 
turn provide an additional driving force for diffusion [69]. 

Another ingredient for the description of an electrochemical system is 
the rate in which ions are deposited on and stripped from an interface. 
The Butler-Volmer equation is a well established form of reaction kinetics 
and relates the current density, to a jump of the electric potential along 
an interface [70]. The most prominent extension of the Butler-Volmer 
equation with respect to mechanics was proposed by Monroe and Newman 
[71], describing a dependence of the exchange current density on the 
mechanical stress of electrode and electrolyte. This serves as a starting 
point for several numerical studies, e.g. [58] [72H74] on morphological 
stability. Only some experimental evidence is available that indicates a 
correlation of interface kinetics and mechanics [75] [76]. 

Various studies have been presented to gain better understanding of 
dendrite growth in battery systems. Chazalviel [77] relates morphological 
stability of lithium with liquid electrolyte with the electrochemical trans- 
port properties of a binary electrolyte arguing that an anion depletion 
and therefore large electric fields are the cause of dendrite growth. His 
model predicts a correlation of dendrite growth time and the time until 
an anion depletion occurs, the so called Sand’s time [29] [78]. Although 
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this behavior is linked to high current densities, similar findings have 
been observed for low current densities [78]. 

Other studies relate reaction kinetics at an SEI with the morphological 
stability of a lithium metal electrode. Barton and Bockris [79] discuss 
the impact of surface energies, which turns out to be of high relevance 
for dendrite radii in the range of 5 nm [80], but not for classical dendrites 
[8I]. The model was extended by non-linear reaction kinetics [82], reeval- 
uated with respect to thermodynamic reference points of PEO [83] and 
extended with mechanics [72]. The work of Monroe and Newman is 
frequently cited and states that stable deposition can be achieved if the 
shear modulus of the separator is twice as high compared to a lithium 
metal electrode, and serves as starting point for more recently published 
studies [73] [84] [85]. However, Barai et al. [86] recognized an unphysical 
mechanical stress state in the analysis of Monroe and Newman [72] and 
conclude that, according to this type of analysis, no stable deposition is 
possible, only the growth speed can be slowed down. A follow up work 
discusses the impact of the non-elastic material properties of lithium 
for deposition [74]. Raj and Wolfenstine [87] describe the nucleation of 
dendrites within a ceramic electrolyte as a function of fracture strength[?] 
Tikekar et al. apply a perturbation analysis to study structured 
electrolytes neglecting elastic stiffness. This work was extended in a 
subsequent works Tikekar et al. [89] [90] by mechanics stating that modest 
Young’s moduli are sufficient for morphological stability. Their result, 
however, was questioned by McMeeking et al. who determine an 
analytic expression for the maximal wavelength of an heterogeneous metal 
electrode in contact with a single ion conductor as a function of elasticity 
and charging rate. Natsiavas et al. [92] deduce with an asymptotic analy- 


sis a positive effect of prestress onto surface roughening during cycling. 


1 Nucleation of Li-ions to lithium metal requires an supply of electrons which is not 
considered by Raj and Wolfenstine [87]. 


12 


1.3 Objectives and Overview 


Klinsmann et al. predict the requirements of ceramic electrolytes to 
resist cracking due to lithium growth. 

Transient models of deposition and stripping during charge and dis- 
charge provide insights into the growing patterns of heterogeneous lithium 
growth. Ferrese et al. combined in the first step the concentrated 
solution theory with an explicit model of boundary velocity and extended 
the model in a second step with the mechanical response of a stiff sep- 
arator [59] [60]. Further approaches to understand dendrite growth on 
different length scales are phase field [94] [95], stochastic and atomistic 
models [97]. 


1.3 Objectives and Overview 


In this work, we seek a deeper understanding of the morphological 
interface stability of a solid state battery with a lithium metal electrode. 
To do so, we have to understand the multi-physics of the interface between 
lithium metal and solid electrolyte as well as the response of the adjacent 
materials, the solid electrolyte and the metal electrode. 

We aim for a sound description for a binary solid electrolyte. Three 
properties of ions within the electrolyte have to be considered. First, the 
ions carry electric charge and thus are affected by an electric field. Second, 
diffusion of the different species has to be described and third, the size of 
the ions has to be accounted for. The former two involve non-idealities 
due to complex multi-component transport and require a concentration- 
dependent parametrization as known from the concentration solution 
theory [53]. Redistribution of ions of certain volume within a host 
material yields to swelling that has to be accounted for. Furthermore, 
the formulation has to capture the effect of deformation caused by e.g., 


an intrusion. 
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Another key aspect to understand protrusion growth is the determination 
of the deposition and stripping rate at the metal electrode. From an 
electrochemical point of view, this is usually done via the Butler-Volmer 
reaction kinetics describing the interface current density. However, it is 
not trivial to answer the questions of whether and how mechanics influence 
the reaction at an interface. We aim for a well defined formulation of the 
Butler-Volmer equation considering both electrochemistry and mechanics. 

Furthermore, we want to describe the deposition and stripping process 
at the metal electrode involving geometrical changes. Together with an 
adjacent solid electrolyte, we expect a deformation of both electrolyte 
and the electrode if deposition occurs inhomogeneously. The complex 
interplay of deformation and electrochemistry shall be described in the 
course of this work. 

Eventually, we want to deduce requirements for the electro-chemo- 
mechanical properties of the solid electrolyte such that morphological 
interface stability is achieved. Two comprehensive studies shall give 
insights into the correlation of mechanics and transport. The findings 
shall help to design robust solid electrolytes which are able to prevent 
harmful protrusion growth and pave the way for the next generation of 
energy storage. 

This work is structured as follows: We briefly introduce some mathe- 
matical preliminaries in the context of finite strain continuum mechanics 
in Chapter |2| We then gradually increase the complexity of the battery 
modeling. 

Chapter [3| presents the derivation of an electro-chemo-mechanical 
transport model of a binary solid electrolyte. Using conservation of mass, 
momentum, energy and entropy in combination with the second law of 
thermodynamics and Gauss’ law, we derive a generic framework for multi- 
component transport in solid electrolytes. We then introduce constitutive 
models for the solid electrolyte and derive a simplified model based on 


the assumption of local electroneutrality in the spirit of concentrated 
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solution theory but extended by mechanical driving forces induced e.g. 
by swelling processes. 

Thereafter, in Chapter [4] we introduce the electrochemical trans- 
port properties that are used throughout this work and analyze the ion 
transport in a solid electrolyte and its electro-chemo-mechanical coupling 
mechanisms with three examples. A blocking electrode setup shows the 
need of a model that solves the Poisson equation. The second exam- 
ple deals with galvanostatic charging and determines the influence of 
elasticity on the transport performance. The third example describes 
an externally deformed solid electrolyte and highlights the need of the 
comprehensive electro-chemo-mechanical model. 

In Chapter |5| we present a rigorous derivation of an extended form of 
the Butler-Volmer equation based on the transition state theory. We put 
special emphasis on the mechanics of deposition and stripping. Further- 
more, we present a novel concept which generalizes the Butler-Volmer 
reaction kinetics by means of electrode, electrolyte and barrier energies. 
This allows us to deduce standard formulations for reaction kinetics in 
liquid and solid battery systems found in literature as well as a consistent 
description of how mechanics influence the reaction kinetics. 

Using the electro-chemo-mechanical transport model and the novel 
interface description, we carry out a first interface stability analysis in 
Chapter [6] We analyze the impact of the elastic stiffness at a perturbed 
interface and study the correlation of the required mechanical stiffness 
for stability and transport parameter using the parametrization of the 
two polymer electrolytes PEO and PS-PEO. Moreover, we vary selected 
transport and geometric parameters to determine their influence on the 
morphological interface stability. 

Chapter [7] deals with the description of the metal electrode. The 
deposition and stripping of material is determined by the interface current. 
Besides the formulation of finite strains within the metal electrode, we 


discuss a second mapping to describe the growth of material. Two possible 
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compatibility conditions between electrode and electrolyte are presented. 
Furthermore, we describe Ohm’s law in finite strain. 

In Chapter we carry out a second stability analysis by which 
the growth of lithium adjacent to a solid electrolyte is simulated. We 
discuss several hypotheses for the trigger mechanism of protrusion growth. 
Furthermore, we investigate the influence of the elastic stiffnesses on 
the growth velocity for both polymer electrolytes PEO and PS-PEO. To 
rationalize the findings, we present a analytical network model for an 
in-plane current that causes protrusion growth. 


We conclude this work and give a brief outlook in Chapter [9] 
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In this chapter, we present mathematical preliminaries to lay the foun- 
dation for the description of solid electrolytes (Chapter |3) and metal 
electrodes (Chapter [7). Since solid electrolytes consist of a mechanically 
stiff host with mobile components, we utilize concepts of solid mechanics 
and mixture theory. We employ the concept of finite strain continuum 
mechanics [98] [99] to describe the kinematics of the host material of the 
solid electrolyte and the metal electrode. Hereby, the mathematical equa- 
tions are written in a reference configuration associated to a undeformed 
material state, which is also known as the Lagrangian point of view. The 
motion of mobile components is commonly modeled within a mixture 
theory [100] in the Eulerian point of view where a fixed control volume 
acts as observer. We are going to combine both ideas by associating the 
fixed control volume of the Eulerian representation with the deformable 
host material. Furthermore, we link the motion of mobile components 
with the Lagrangian description. We note that this chapter is partially 
based on the work of Ganser et al. [101]. 


2.1 Kinematics of a Solid Material 


2.1.1 Deformation Gradient 


Following Holzapfel [98] and Truesdell and Noll [99], we employ the frame- 
work of finite strain continuum mechanics, and describe the deformation 


of a solid material at time t € T by the mapping + between a reference 
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deformed geometry 


undeformed geometry 


al. o TEP CB: 


Figure 2.1: The mapping + correlates undeformed (material) and deformed (spatial) 
geometry Bo and $+, respectively. Each property in the spatial description such as the 
volume element dv, area element da or normal vector n is related via the deformation 
gradient F to the material description of the corresponding properties dV, dA and N. 


configuration By C R3 with material coordinates X € Bo and deformed 
configuration B, C RÌ with spatial coordinates x € B,, such that 


Bo xT => Bt 
Pt = (2.1) 
(X,t) > z= p(X), 
as seen in Fig. Accordingly, the material and spatial velocity of the 


material are given by 


Vo = AX) oX) and w=éder a). (22) 


The displacement of the material is u = x—X. To describe local stretches 
and rotations, we make use of the deformation gradient F' defined by 
Op. Ox Ox; 


F 
OX; 


18 


2.1 Kinematics of a Solid Material 


The determinant J = det F thereby describes the mapping between ma- 
terial volume elements dV C Bo and their deformed spatial counterparts 
dv C Bi, i.e. 


— = det F =: J > 0, (2.4) 


where the positivity constraint is required to avoid interpenetration of 
matter. The deformation gradient can be uniquely decomposed into a 


rotation and a stretch according to 
F = RU, (2.5) 


where R is a orthogonal rotation tensor and U is the right stretch 
tensor. The former describes rigid rotation (RTR = I) and the latter is 
positive definite, symmetric (U = U E) and measures the local stretch 


(or contraction) along the orthogonal eigenvectors. 


2.1.2 Volume and Surface Transformations 


Note that all quantities and equations can be formulated both in the 
referential and the spatial setting. Since the physical interpretation is 
more intuitive in the spatial setting, we will usually formulate governing 
equations in the spatial form and transform them to the material for- 
mulation, which is the natural setting for a numerical treatment. From 
here on, whenever possible, we use lower case symbols for the spatial 
quantities (e.g. £, Ca) and upper case symbols or diamond superscript 


(e.g. X, c2) for the material quantities. 
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With the help of the deformation gradient F, one can describe the 
mapping of tangential and normal vectors between material configuration 


(N,T) and spatial configuration (n,t) as 
t=FT and n=FTN, (2.6) 


see e.g. Holzapfel [98]. The integrals in the material and spatial config- 
uration coincide for properly transformed quantities. The volume and 
the surface integrals (with spatial and material surface element da C OP; 


and dA C Po, respectively) thus are 
f ydv= | YdV and / y - nda = Y -NdA, (2.7) 
Pt Po OP; aPo 


where y is a scalar or vector field defined on the domain P, and y a 
vector field defined on the boundary OP;, respectively. The corresponding 


material properties are 
Y=Jy and Y=-WrT". (2.8) 


The latter is known as Nanson’s relationship. 

Now consider a scalar quantity y defined on a deformed surface. Its 
scalar material counterpart Y will differ because the effective size of 
the material area element dA is different to the spatial area element da. 


However, the identity 


/ YdA= gda (2.9) 
OPo OP; 
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has to apply. To obtain a relation between ĝ and Y, we insert N-N =1 
on the left hand side and n-n = 1 on the right hand side. We then use 
Eq. (2.7) and multiply both sides with N to obtain 


Y =JnF"* . Nj= 1 5%: (2.10) 


F:(n@N 
The normals N and n are both normalized. Hence we have n 4 FTN 
and Eq. (2.10) cannot be simplified further. 

Note that material points in Py are usually associated with an unde- 
formed material. Since Po is a reference description, it is independent of 


time and thus the integral operator and time derivatives are interchange- 
able 


d d . 
— f YdV = —Y dV = Y dV. (2.11) 
dt Po Po dt Po 


The property of Eq. (2.11) does not hold for an integral formulated in 


the spatial configuration. 


2.2 Kinematics of Mobile Components 


We describe the motion of mobile components with respect to a control 
volume associated to a host material in both spatial and material config- 
urations. Fig. [2.2] gives an overview of the various transformation steps 


which will be discussed in the following. 


2.2.1 Spatial Description 


We assume that a species a is distributed as a continuum in a control 


volume P? C B. Postulating that PP coincides at time t with a control 
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Figure 2.2: Transformation steps: @ Host domain P; and species domain P¿* underly 
a different motion during a timestep ôt. O With the Reynold’s theorem, one considers 
the change of concentration within the host material in terms of a mass flux ją. O All 
quantities can be expressed with respect to a material, undeformed configuration Po. 


volume P; associated to a host material moving at velocity vo = ®t, 


volume and surface integrals coincide at time t: 


I ydva =| ydv and L Y- ndaa =| J: nda, (2.12) 
Pr Pı OPE OP; 


If the number of moles of species a in the control volume element dv is 


given by dN,, we can define the spatial molar concentration as 


(2.13) 


The spatial concentration cy describes the amount of species of type a 
per deformed volume dv. 

Now, let the species in P? move at velocity va, which is in general 
different from vo. If we consider the rates, we have to take the different 


time-dependencies of the integration domains Pf and P; into account 
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(0 in Fig. |2.2). To describe the relative motion as seen in € of Fig. [2.2 


we use Reynold’s theorem 


d d 
— ydva = — ydv +/ y ((Va — vo) +n) da 
a (2.14) 
= — E ==, -n 
dt Jp, yav Jr. an? Ja a, 
where we have introduced the quantity 
Ja = Ca (Va — Vo) (2.15) 


which will turn out to be the molar mass flux of species a. The second 
term on the right hand side of Eq. (2.14) denotes the influx of volume 
quantity y into the host domain due to a relative velocity of the species 


with respect to the host. 


2.2.2 Material Description 


The third transformation in Fig. [2.2] is necessary to obtain a continuum 
description of the mobile components with respect to the referential 
configuration of the host material. Therefore, we define a material 
concentration similar to Eq. and link it to the spatial concentration 
by Eq. (2.4), which yields 

dNa 


C= IV = Jca. (2.16) 


Furthermore, we use Eq. (2.8) to transform the spatial quantity j, from 


Eq. (2.15) into the material configuration with 


Ja = JF tja. (2.17) 
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2.2.3 Reynolds Theorem in Finite Deformations 


Combining Eq. (2.7) with Eq. (2.14) and using 4 = Ju _ gives us 


cz 


the transformation to link moving species within a reference description 
of matter via 
d 


. 1 
£ | yda= | Yav+ n Y (Ja: N) dA. (2.18) 
dt Jpo Po OPo Cà 


The divergence theorem f. ap, Y  NdA= Sp, Div (7) dV with the def- 


inition of the divergence operator Div (Y) = ada yields a fully 
volumetric description with 

d . . 1 

— ydvg = Y + Div | —Y 8 J, | dV 

dt Jpe Po ca 


(2.19) 


= | Y +Grad (žr) Ja + + yDiv (Ja) dV, 
c Ce, 


Po Q 


where we have applied the chain rule to obtain the second line. Eq. (2.19) 
relates the rate of change of a quantity y in a moving domain Pf to an 


undeformed reference state. 
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Binary Solid Electrolyte 


Within this chapter, we establish a consistent model for ion transport 
in a deformable host. Based on the work of Ganser et al. [IOI], we 
exploit concepts of finite strain continuum mechanics, rigorously consider 
electrostatic influences and use standard notation of electrochemistry. 
With the help of the thermodynamic framework of Helmholtz energy and 
a series of consecutive steps deriving the physics and electrochemistry, 
a clear and consistent picture of the interdependencies of electrochem- 
istry and mechanics is drawn, leading to a combination of ideas from 
multicomponent mass flow under the influence of an electric field with 
finite deformation theory for the mechanical response of the host. Both 
a treatment based on Gauss’ law for a system containing unbalanced 
charges and an electroneutral version for a binary electrolyte are pre- 


sented, thereby establishing a link to the well known concentrated solution 


| lon Transport ( 


(migration / diffusion) g | Ertemälfordes 


Mechanics 
(elasticity etc.) 


| Stress driven diffusion Geometry all A 


Figure 3.1: Representation of the coupling between multicomponent ion transport and 
mechanics in a solid state electrolyte. Mechanical response in the form of geometry 
changes and stresses is triggered by swelling or external forces. The ions are subject 
to this change of geometry and mechanical driving forces give rise to stress driven 
diffusion. 
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theory similar to Monroe and Delacourt by introducing local 
electroneutrality and by linking the commonly used transport parameters 
conductivity, transference number, diffusivity and thermodynamic factor 
to mobilities and chemical potentials. Non-idealities usually associated 
with a concentrated solution are thereby accounted for. 

Throughout the derivation, specific care is taken in considering coupling 
mechanisms between ion transport and mechanics. In the spirit of Bucci 
et al. [68], we include coupling mechanisms such as swelling, stress driven 
diffusion and the change of geometry on which the transport takes place, 
see Fig. [5.3] Similar couplings between transport and mechanics are 
also present in hydrogels, where small molecules diffuse in a network of 
long fibers and where a change of local component concentration leads to 
heterogeneous swelling or shrinkage of the host material and consequently 
to mechanical stresses, which can in turn provide an additional driving 
force for diffusion [69]. Similar relations between swelling and transport 
are also found in the context of battery active materials [61167]. 

In the following, we present the fundamentals of conservation laws 
in electro-chemo-mechanical systems, introduce Helmholtz energies and 
their relation to constitutive material properties, discuss the swelling 
mechanism and introduce thermodynamically consistent mass fluxes. 
Eventually, we conclude the chapter with a theory for binary solid elec- 
trolytes similar to those in concentrated solution theory [53] but extended 


by mechanical driving forces. 


3.1 Field Variables 


For the description of the fully coupled electro-chemo-mechanical system 
we distinguish between a solid host a = 0 and N species with indices 


a= 1...N that can move relative to the host. The electroneutral host 
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F 
— 


undeformed geometry deformed geometry 


X € Po C Bo 
o; 


1 
g N E |a 
Ja TE j 


concentration el. potential deformation concentration el. potential deformation 


material description spatial description 


Figure 3.2: An undeformed geometry provides the material description where all 
balance equations are solved. The deformed geometry as well as spatial field variables 
are linked to material counterparts via the deformation gradient F. 


provides the mechanical strength whereas the (ionic) mobile components 
act as charge carrier. 

As introduced in Section [2.1] we Characterize the kinetics of the host 
with a mapping pS”, see Eq. (2.1). The deformation gradient F° of 
Eq. (2.3) depicts local stretches and rotations. For the sake of simple 
notation, we will drop the index for the deformation gradient and continue 
with F in this chapter. We emphasize that the host displacement with 
velocity $?” = vo will serve as reference to which the species transport 
will be defined. 

Now we shift our attention to the solutes. We assume that N species 
are distributed each as a continuum in the solid matrix consistent with a 
mixture theory. If the number of moles of species a in an undeformed 
volume element dV is given by dN, and its molar mass is Ma, we can 
define together with Eq. (2.16) the partial material densities as 

= dNaMa _ 


o =Ma =1.N. al 
Pa WV Ca Va (3.1) 
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The spatial equivalent follows a similar definition and is linked to the 
spatial concentration by Eq. (2.13) and to the material density by 


Eq. (2.16), yielding 


dNa Ma 
= Ma Va =1..N. (3.2) 


Pa dv J 


We assume that at each point of the host material, species œ moves with 
velocity va. Recall that the velocity of the host is defined via Eq. (2.2) 
and generally va 4 +. The overall densities of the mixture, including the 
host material in the undeformed and deformed configurations, respectively, 


are given by 


N N 
Prot = 5 Pa Ptot = > Pa- (3.3) 
a=0 a=0 


Since species œ = 1..N can carry charges, whereas the host material is 
assumed to be charge neutral, the material and spatial charge densities 
for each species are 

o 


Pog = Fz, Pga E26; (3.4) 


where Za is the charge number associated with ion a. The total charge 


density is computed via 


N N 
p=. He > Pa (3.5) 
a=1 a=1 


for the material and spatial configuration, respectively. 
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3.2 Balance Principles 


3.2.1 Balance of Mass 


For the balance of mass, we distinguish between the host (a = 0) and 
the chemical species (a = 1...N). The host material is assumed to have 


constant mass which yields [98] 


pe =0. (3.6) 


To derive the balance of mass for a single species a, we consider at time t 
an arbitrary control volume Pf C B, with species velocity Va. The 
only change in mass of a in that volume occurs due to a mass source rq 


associated with chemical reaction and we can write 


d 
— Cadua = | Tadva Va #0. (3.7) 
t PE 


To rewrite this in terms of the control volume P; tied to the motion of 
the host, we use the Reynolds theorem (2.14) and obtain 


e Cadu + a Ca (da — vo) nda = [ Tadu Va #0. (3.8) 
t t 

The second term on the left hand side denotes the mass flux 
Ja = Ca (Va — Vo) through the boundary OP; as introduced in Eq. (2.15). 
Hence, we describe naturally the mass flux of a species relative to the 
velocity of the host vg. That distinguishes our description from multi- 
component transport in liquids which usually utilizes the molar mixture 
velocity [103], barycentric velocity [104] or the solvent velocity [102]. 
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Now, we apply a transformation to the material setting, the divergence 
and localization theorem, to obtain the local form of the balance of species 


mass and total mass, respectively, as 
& + Div (Ja) = Ra Va #0 => fa + Div (Sto) =0. (3.9) 


The material species and total molar flux are J, = JF”' Jos see 
Eq. (2.17), and Jiot = E Ja. The material concentration source is 
Ra = Jra and since a chemical reaction does not produce mass, we have 
oe Ra = 0, which gives the second equation. 


a=1 


3.2.2 Balance of Charge and Gauss’ law 


As we assume that charge is exclusively transported by ions without 
involvement of electrons, the balance of charge can be directly derived 


from the balance of mass 
pe +Div(Jq,)=Ry, Va => ¿¿+Div(J)=0, (3.10) 


with charge flux Jq, = FzaJa and source Ry, = Fz Raq. Since charge 
cannot be generated in reactions, we have Spa Ra, = 0, which — 
together with the definition of current J, = a J,,,— yields the total 
balance of charge on the right hand side. Note that Eq. (3.10) is not 
independent of the balance of mass as seen in Eq. (3.9). As a result, only 
N — 1 mass fluxes Ja are independent for a given total current flux Jq. 
The Gauss law relates the spatial electric displacement d to the total 
charge density p¿. The spatial and material representations are 
div(d) =p, # Div(D) = pz (3.11) 
with the material electric displacement D = JF~'d. Further we follow 
McMeeking and Landis [105] and introduce a surface charge w° = [D] - N 
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as the jump of electric displacement. The electric field e is related to the 


magnetic field by curl (e) = ib = 0 due to Maxwell’s equations, where 


we have assumed that the change of magnetic field b in time can be 
neglected. To satisfy Maxwell’s equatior"] we define the electric potential 


® such that 


e = —grad (9) = -F "Grad (®) = FE (3.12) 


with E = —Grad(®) the material electric field. The electric field is 


related to the electric displacement by 
d = eye + p = €0€;€, (3.13) 


where p is the polarization. The quantities ey and e, are the vacuum and 
relative permittivity, respectively. By inserting Eq. (3.12) in Eq. (3.11), 


we obtain the well known Poisson equation 
—div (energrad (®)) = pg. (3.14) 


In the case of electroneutrality with pg = 0 and in the absence of chemical 
and mechanical driving forces, the electric potential is decoupled from 


individual mobile components and Ohm's law 
jq = Ke = —Kgrad (®) (3.15) 


in terms of the total current density ją becomes valid. In such circum- 
stances the current is driven by an electric field through a conductivity «K. 
An extended version of Eq. that is valid for solid binary electrolytes 
will be presented in Section [3.5] 


1 Here we make use of the identity curl (grad (£)) =0 VE. 


31 


3 Theory of Ion Transport in a Binary Solid Electrolyte 


3.2.3 Balance of Linear Momentum 


To state the balance of linear momentum, we sum the rates of the partial 
linear momenta of host and species and equate them with mechanical 
(e.g. gravity), electrical and chemical body forces of each species bM, bE, 


bE and effective tractions tM, t®, tC acting on the host material and 


obtain 

d N N 

Sy] Padadva = >| BM bE + bode | tM +t da. 
g a=0 1 Pr a=0 Pt IP, 


(3.16) 


We now have to carefully apply the time derivative on the left hand side. 
Therefore, we apply Reynolds theorem in finite deformations (Eq. (2.19)) 
and obtain for the left hand side of Eq. (3.16) 


d 
al Padadva =| L (p Va) + Div (Ma VaJa) dV 
di Je mt 


= / PV. + Grad (Ma Va) - Ja + MaVaRa dV, 
Po —— MMS, 
BY 


a 


(3.17) 


where we have utilized the balance of mass (Eq. (3.9)) for the second 
line. The body force BY corresponds to an inertia force of component a 
within the host material. 

The total electric body force BE, = Nos JbE is expressed as a Piola- 
type Maxwell stress tensor BE, := Div (PP). The same procedure is used 


for the chemical body forces. Further, we follow McMeeking and Landis 
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105| and state that the mechanical, electric and chemical tractions are 


balanced via a jump in mechanical stress, yielding 
/ +t + toda = TM TES TOdA 
OP; OPo 


==) [PM] - NdA (3.18) 


= / Div (P™) av. 
Po 


Consequently, an expression for the balance of momentum of the host 


material is 


N 
PRP = BY, + Div | PMPP +P] -XBY (3.19) 
A 
P a=1 


where we have introduced the effective total mechanical body force 
N N 
Bie = X BY = > db, (3.20) 
a=0 


a=0 


and the total Piola stress P. The inertia, advection and reaction contri- 
bution of the mobile components can be regarded as an additional body 


force BY acting on the host material. 


3.2.4 Balance of Angular Momentum 


The balance of angular momentum 


d 
ge=M (3.21) 
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states that the rate of change of the angular momentum 


L= >| LX PaValVa (3.22) 


is equal to the total torque 


N 
M= 7 x x BEOM + o dy, +/ x x PM ga (3.23) 
a=0 PE OP; 


with the body force bE™ = == bE + b + b™ and traction force 
tECM — M + ¿E + tC. Further, an angular momentum due to dipoles has 
to be considered. Following McMeeking and Landis [105], we use o® = 
a = €sj kkj» where &;;x is the Levi-Civita symbol. Although a deviatoric 
contribution of the chemical stress is unlikely, we set 0° = EH Ong and 
obtain by combination with the electrostatic contribution o&° = oË + of. 


The rate of angular momentum computes to 


N 
d d ; (1 Š 
Ge = 2 f 5 (X x paVa) + Div (4 (X x p2Va) ® Ja) dV 


N 
d 
= ~/ EijkXj È (oV) + Divm (MVE Jm) + MacindimVe Ja dV, 
(3.24) 


where we have used the index notation in the second line. The applied 


torque in index notation is 


N 
M — 27 coto | 
a=0 Pt 


x x (0 n)da +f cross dv 
OPt 


Pt 
N 

= | Eijk Xj (apor + Divm (ae) + €:jkOmj PH Fmi + €ijk PRS FjdV. 
a=0 Po 


(3.25) 
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Combining Eq. (3.24) and Eq. (3.25) gives 
a d 
5 eik Xi E (PEVE) + Divm (Ma Vè Jm) — BE™ — Divm (it) dv 
Po 


>, Eijk (fa En — MVE J?) dv. 
Po 

(3.26) 
The left hand side is equal to zero due to the balance of linear momentum, 


see Eq. (3.17), and BEC = BM + Div (PS + P”); The right hand side 


results in the symmetry condition 


N N 
- Y MaJa 8 Va = FPT -X Ma Va 8 Ja. (3.27) 


a=1 a=1 


Without component motion this reduces to 
PF’ = FP’, (3.28) 


which is the standard result for balance of angular momentum in finite 


strain solid mechanics. 


3.2.5 Balance of Energy 


Consider the balance of energy in an arbitrary volume P; C B+. The 
change of total internal energy E'"*, as well as kinetic energy ER, is 
equal to the power supply in P, due to heat LĦ, mechanics LM and 
electricity LE: 

d 


z (qn + EX) = LE + LM + [E (3.29) 
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The left hand side of Eq. (3.29) describes the rate of internal energy ea 


and kinetic energy 4 PaVa * Va of each species 


N 
Em = >; Palade, 
a=0 E 


2 1 
ES? = >, ¿Pava "Vadva: 
a=0 E 


Note since each component œ moves with its own velocity vg, we have to 
again apply Reynolds transport theorem through Eq. (2.14) and obtain 


d int ( d f / 
Zp = Paeadv + Pata (Va — vo) nda 
dt 2 dt Pr OPt 
a d 
z f qq (Pata) + MaGrad (ea) - Ja + Ma Rata — MataeadV, 
a=0 Po 


(3.31) 


where we have used Eq. (2.19) together with Eq. (3.9) to obtain the 
second version, involving the transformation to the material configuration, 
the divergence theorem and balance of mass. The same procedure applies 


to the kinetic energy 


d A 1 
-p| = BY .V,-—V,:-V.M, 32 
ai 5 K Va = ¿Va? VaMaRadV, (3.32) 


a=0 Po 


where we recall the definition of BY in Eq. (3.17). On the right hand 


side of Eq. (3.29), we observe energy production due to heat sources 
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and heat fluxes. In terms of heat we treat host and components as one 


system and get the contribution 


i =f w- | qu. nda, (3.33) 


where we have introduced a homogenized heat flux j E and an external 
heat source rf, e.g. an electromagnetic heating. Note that heat generated 
due to chemical reactions is an internal process of conversion within the 


internal energy. The reference description of Eq. (3.33) is 
L” si RË — Div (1%) dv. (3.34) 
Po 


The rate at which mechanical work is done on the host material and the 


species is 


N 
EM =v f (be - va) a+ | (EM. vo) da, (3.35) 
a=0" rt % 


where only the host material can support traction forces. The individual 
body forces act on the host material as well as on each mobile component. 


The result converts to 


ey, bi (va — vo) de + > f bi! vodu + f t™ . voda, 


OP: 
= > FB PE ane (Bi, + Div (P)) -Vo +P: Fav, 
Po a=0 2 
(3.36) 


where we have made use of Eq. (3.18), the divergence theorem and the 
notation of double contraction, meaning A:B = A;,;B;; to obtain the 


version with respect to the reference configuration. 
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The electrical energy supply due to transport of charged species at an 


electric potential ® by means of a current j, reads 


E =- bj, nda= | e: jg — Pdiv (jq) dv. (3.37) 


Transformation to the reference formulation yields 


IE -f + Jq + Div (D) av 
Po 


(3.38) 
=| z. Sg + bay + f SJD] -NdA 
Po OPo 


where we have used Div (Jo + D) = 0, a consequence of Eq. (3.10) 
and Eq. (3.11). The last term can be identified as a surface charge 
wSC = [D] - N, see e.g. [105]. Since all free charges are carried by ions 
and the host is electroneutral, we assume that the change of surface 
charge is zero and hence wSC = [DJ]: N = 0 on ôPo. 


Combining Eq. (3.19), (3.32), (3.34), (3.36), (3.38), we obtain the local 


material form for the balance of energy 


N 
m Pike =R" — Div (3%) +P:F+E-D+E-J, 
7 N N 
+Y Maea- Y (Ba + Grad (Maea)) Ja (3.39) 
a=1 a=1 


a 1 
u 5 (Maea u 5 MaVa va) Ra 
a=1 


with the extended specific body force B, = -2F* (BM + BY). In 
addition to the thermal and mechanical contributions, we see a power 
supply due to electric interaction, e.g. Joule heating, species transport 


and chemical reaction. 
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3.2.6 Balance of Entropy 


We proceed similarly with the balance of entropy and sum over the 
component entropies Na. We locally assume equal absolute temperature 


O = 0, for each constituent and obtain 


dw pH j” podv 
— PaNadva = dv - nda +f —— dv 3.40 
dt vf p, O ap, 9 Pp, O ) 


where dy is the dissipation density. Applying Reynold’s theorem through 
Eq. (2.14) on the left hand side, divergence theorem on the right hand 


side and switching to the material frame gives 


N 
d o o i 1 H : H 1 H 
Pat pov + 3 (ce Div (J )) +57 Grad (@) 
e (3.41) 
— $0 [Grad (Mana) Ja — Manada + Mana Ra] 
a=1 


as an expression for the rate of total entropy. 


3.3 2nd Law of Thermodynamics 


We now introduce the total Helmholtz energies Y and the species 
Helmholtz energy Ya from the Legendre transforms of the internal 


energies eg 


N N 
D(F, c, O, D) =) Rta = >> p (ea — Ona). (3.42) 
a=0 a=0 
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Combining Eq. (3.39) and Eq. (3.41), making use of Eq. (3.42) and 
solving for p§dv, we obtain the so-called full dissipation inequality, also 


known as the Clausius-Planck inequality [98] as 


N av. au \ . 
podv = (En aje ' (» se) 


a=0 
. $ : 
olas Ow 
M, = al 
> ( ara =) dl OD 
N ~ 
=>» [Ba + Grad (Matpa) + FzaE — (Mana) Grad ()] Ja 
a=1 
— ¿Grad (9) - J — ( Matha — 5MaVa ` Va ) Ra 
>0. 


(3.43) 


The left hand side has to be non-negative with pjdy > 0 due to the 
second law of thermodynamics. Following the Coleman-Noll argument 
[106], we assume independence of the rates 0, F, Cc; D and therefore 
require the parenthetical terms in Eq. to vanish. This yields the 


consistency conditions 


x a Ow 
5 Na = tar: P = ap? 
09 OF 
a=0 . i (3.44) 
_ oy e 
Ha := Mata EA aD’ 


where we have introduced the chemical potential (per mole) as 
Ha := Maa of species a, which is also sometimes referred to as the 
convective potential [107], a measure for the energy required to insert 


a mole of species a into a fixed volume dV. Consequently the total 
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Helmholtz energy of Eq. (3.42) is a combination of the host energy and 


the sum of the chemical potentials of the components 
_ N 
b= poyo + y Rha. (3.45) 
a=1 


Further, since Y is composed from of the individual Ya, the definition in 
Eq. (3.44) has an implicit structure and can be resolved via 


9 Y N 
= B_ 
Mata = ace, 2,250 > 2% ace == 
j ~ (3.46) 
OMe Ovo 
> = — 
= 2% Oc ae’ 


where the last representation can be seen as a Gibbs-Duhem like equation 
at constant temperature and deformation. 

We continue with the evaluation of the dissipation inequality (3.43). If 
all equations of Eq. are fulfilled and Joules heating is reformulated 
to 


N 
2- J = — Y Grad (8) - FzaJa, (3.47) 
a=1 
we can reduce Eq. to 


0<- > (Grad (Ba FaF P20) + Mano Grad (9) Ja 
a=1 


(3.48) 
1 P 1 
= Grad (6) J" - > (n - ¿MaVa: va) Ra. 


Eq. (3.48) constitutes a thermodynamic constraint for the molar fluxes 
Ja (and therefore species velocities va), the heat flux J" and the 
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reactions Ra. Under the assumption of each process being indepen- 
dent, the first term on the right hand side of Eq. (3.48) relates molar 
flux to chemo-electrical driving forces, inertial and body forces. Similarly, 
when independent of other terms, the second term leads to the well 
known Fourier’s law of thermal conduction and the third term guides the 
direction of reactions. 

Assuming uniform temperature (Grad (9) = 0), quasistatic conditions 
for the host (2% ~ 0), the absence of body forces (BM = 0) and further 
assuming that inertia related forces are small (BY ~ 0), we obtain the 


thermodynamic constraints 


N N 
0<-) Grad (wa): Ja and = 0< -S paRa, (3.49) 
a=1 


a=1 


where we have introduced the electrochemical potential 
Wa = Ha + Fza®. (3.50) 


The first of Eq. (3.49) provides a constraint for all mass fluxes Ja as 
discussed in the upcoming sections. The second of Eq. (3.49) constitutes 
a requirement for chemical reactions within the bulk, e.g. degeneration 
of electrolyte. This will not play a role in the following sections, where 


Ra = 0 will be assumed. 


3.4 Material Description 


To complete the electro-chemo-mechanical description of ion transport, 
we have to specify a material model for the solid electrolyte. Therefore, 
we will introduce the swelling mechanism, define the Helmholtz energy 


including mechanical, electrical and mechanical effects and derive the 
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Figure 3.3: The deformation gradient is composed of purely isotropic swelling Fs due 
to the pressure of ions and elastic deformation Fe. 


constitutive equations and prescribe thermodynamically consistent mass 
fluxes. 

In the following, we use the a decomposition of the deformation gradient 
into swelling and elastic contributions as seen in Fig. [3.3| with the 


multiplicative split 
F=F.F,, (3.51) 


where Fe describes the elastic response of the material and F, corresponds 


to a volumetric stress free expansion or compression due to swelling. 


3.4.1 Swelling of the Host Material 


The volumetric phenomenon of swelling occurs when species are inserted, 
removed or relocated in a host material. Although we will assume that 
swelling itself is stress free, it can eventually trigger an elastic response 
due to constraints from the boundary or adjacent material. Take for 
example a sponge. If it is clamped and cannot expand but additional 
species are inserted, then either the host or the species are compressed. 


The opposite effect is also possible when the host is already infiltrated 
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by species and is again clamped. If species are removed, then the host 
has to account for the missing material which eventually leads to tension 
forces. 

Here, we assume an isotropic swelling. The deformation gradient of 


swelling then is 
F, =J31 (3.52) 


with the Jacobian Js = det F, as an unknown function of i.e. the 
concentration. With Eq. (2.4) we can deduce from a given material 


volume element dV the spatial volume element dv with 
dv = JeJsdV = JdV, (3.53) 


where J. = det Fe is the Jacobian of elastic deformation. Each of the 
volume elements can be separated into a contribution from the host and 
from the mobile species a. The unswollen and undeformed state is defined 
by 

N 


dV = dVo + Qad NZ, (3.54) 


a=1 


where dN3% = dV ce! denotes the amount of species a in dV and Qa is 
its molar volume, both defined in the in the initial, reference state. The 
swollen and deformed volume is 
N 
dv = dvo +) NadNa, (3.55) 
asl 
where we assume that the deformation of the host volume follows 


dvo = JdVo and that the spatial molar volume of species a is given 
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by Ôa. Further, the amount of species dN, in volume dv is determined 


from the concentration either in material or spatial configuration as 
dNa = Cadv = £ dV. (3.56) 
From Eq. (3.53), Eq. (3.55) and the knowledge that the host material 


deforms only elastically we obtain 


N 
JeJsdV = Jed Vo + Y QadNa. (3.57) 


a=1 


Rearranging Eq. (3.57), using du) = dV — o NactdV from 
Eq. (3.54) and with Eq. (3.56) then yields 


j 1- N Qa ref 
Js a? (Pas, =e -Q a) = Lo aa (3.58) 
1= Y, Mata 


where the second equality is obtained with x = $, An explicit model for 
the deformed molar volume Qa of species a is necessary to evaluate the 
swelling of a material as seen in Eq. (3.58). We consider two assumptions 
for the deformation of a species. 

The first assumption is incompressibility of the inserted species such 
that Ôa = Qa. This gives 
1- aa Ove 


acl 3.59 
1- DI, Qala nn 


J pei aaa _ 
$ = 


and implies an upper bound for the ra with yL 9 Calla < 1. 


The denominator reaches zero for Ca — which leads to an infinite 


> ’ 
swelling (J, — inf). The volume is completely occupied by the species a 


in this case and the host is compressed to zero volume. 
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incompressible incomp. (linearized) compressible 
ref 
1-5 Qa caf N Qa (ca -ch ) N s o ref 
Js 1-3, lala 1 + Yazi 1-Qacref 1 + Pozi B Qe (c2 Ca ) 
AS Js 1 BS 
1-La=1 Maco Gd Rack) 


Table 3.1: The expressions for J; and AS for incompressible, linearized incompressible 
and compressible models of isotropic swelling. 


Developing Eq. (3.59) with a Taylor expansion around ca = Œf yields 
the linearized form 
N Q 
incompressible __ Q _ pref 
Js| =1+ ) I net (ca — cy") (3.60) 


a=1 


which approximates swelling as a linear function of the spatial concentra- 
tion. 

The alternative assumption is compressibility of the inserted species 
with Qu = BJ... The molar volume increases or decreases proportional 
to the mechanical deformation of the host. The factor 85 thereby accounts 


for non-ideality in swelling and could be determined experimentally or 
by atomistic simulations [64]. Eq. (3.58) then yields 


N 
er =1+ Be 5 Qa (a _ eS) (3.61) 
a=1 


and aligns with models in literature [108] [109]. Using \° as defined in 
Tab. [3.1] we compute the derivative of Js with respect to the variable £ as 
0704 


ðJ: e Ss 
gE N ae (3.62) 


Tab. [3.1]sums up the expressions for Js and AS which follow from the three 


model assumptions of swelling. Note that the reference concentration 
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——*— incompressible 
| —©— incompressible (linearized) | 
2 —+— compressible 
S 1.54 
1 s 
0 1000 2000 3000 4000 


c [mol/m?] 


Figure 3.4: The graph illustrates the correlation of swelling and the concentration in 
a clamped system (J = 1) infiltrated by a single species with concentration c with 
cf = 1000 me and Q = 140 7). The incompressible version of Eq. (3.59) diverges 
at > = 7143 mg, No significant difference is observed for the compressible (Eq. (3.61) 
with 8° = 1) and the linearized incompressible case (Eq. (3.60)). 


has to obey > Mac < 1 due to the requirement of Js > 0, for all 
of the three swelling assumptions. Fig. [3.4] illustrates the concentration 
dependencies of the two incompressible and the compressible version. 
The response is very similar around c' for all three assumptions. The 
effect of incompressibility is, however, significant for highly concentrated 
solutions. To reduce complexity, especially due to the factor AS, and 
to be in alignment with similar theories [68] [69], we will consider the 


compressible version of Eq. (3.61) in the upcoming sections. 
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3.4.2 Constitutive Helmholtz Energy 


We assume that the Helmholtz free energy of the electro-chemo-mechanical 
system can be written as the sum of mechanical, electrical and chemical 


contributions as 
~ ~ ~ N ~ 
Ù (F, c, D) = 9M (F,c) +0" (F,D)+ > YSF, &), (3.63) 
a=1 


where we have chosen the Helmholtz free energy of the host at reference 


state to be zero. 


Electrostatics 


For the sake of simplicity, we use an isotropic formulation of electrostatics 
in the spatial domain consistent with Eq. (3.13), which leads to 


pel 
J ` 2eoer 


d-d. (3.64) 


The reference description is then given as 


-E 1 


= zg (5) (DSD) with € =JF7TleFY. (3.65) 


We introduce the material description of the relative electric permittivity 
as €f to keep the notation concise. With Eq. (3.65) we do not include 
dissipative materials [IIO] or anisotropic polarization and assume a 


concentration independent permittivity. 


Mechanics 


In addition to the swelling of Eq. (3.52), we assume a Neo-Hookean 


material law for the description of the elastic response of the host. The 
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mechanical energy, which depends solely on Fe = FF ,'(c?), then 


S a 


becomes 


UF, cà) = RA (Fo) 


i i A (3.66) 
= ¿NH (I; — 3) + ANH (log I3)" — ynu log Is, 
with the Lamé constants ynu, Any and the invariants 
l=tr (ETF) , b= Je = det Fe, (3.67) 


see e.g. Ogden [III]. Note that the Lamé constants can be converted 


from Young’s modulus E or shear modulus G and the Poisson ratio v by 


Ev 2Gv E 
uni e Wer) OO 


ANH = 7 


Hereby, E, G and y have their conventional meaning relating infinitesimal 


strain and change in stress. 


Chemistry 


The chemical contribution to the Helmholtz energy is assumed to be 


7 Ca fala 
ve =cp,+ROf In ( Ta deg. (3.69) 


ref 
Ca 


with the gas constant R and the standard chemical potential u, which is 
independent of concentration and deformation. In concentrated solutions, 
a correction factor, the so called fugacity or activity coefficient fa is 
introduced to incorporate non-idealities in the formulation [112] 113]. In 
dilute solutions, an ideal gas response is obtained with fa = 1. 
Definition is formulated very similarly to classical mixture theo- 
ries in fluids and solids, but there are two features we want to highlight. 


In contrast to conventional usage in solids, where a vacancy motivated 
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undeformed deformed 


DE (ca)lts < PE (ca)lta 


Figure 3.5: Consider an isolated system: The effective concentration cq changes 
from the undeformed (tı) to the deformed state (t2). Therefore, the chemical energy 
increases and gives rise to chemical stresses t© also known as osmotic stress or pressure. 


chemical potential is a function of material concentration c), [63], we 
follow a concept used for fluids and postulate that the chemical energy 
depends on the amount of ions per deformed volume, i.e. on the spatial 
concentration Ca. Therefore, material compression leads to a decrease of 
spatial volume and hence an increase of effective spatial concentration 
(Fig. [??). The importance of this assumption and its implications will be 
illustrated in Example [4.5] The second feature is the integration domain 
of the integral in Eq. (3.69) which represents the entropic contribution to 
the chemical energy. We choose it to be zero at reference concentration 
cif. This will be beneficial for the description of the mechanical response 


in the upcoming section. 


3.4.3 Constitutive Equations 


Following Eq. (3.44), we can deduce the stress tensor, the chemical 
potential and the electric field from Eq. (3.65), Eq. (3.66) and Eq. (3.69). 


The Piola stress tensor then follows as 


P=PM+p®+Pp° (3.70) 
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with the three contributions 


pa =i Ann log Is — qn) Fo” + ynu Fo), (3.71) 
1 
Pe- FT (B@D-5(E-0)1), (3.72) 
= ca əl f. 
P= X` -F "Re 1 ” | de®. 378 
2 a + oa] Ca 20) 


Note that the electric contribution P® is the Maxwell stress tensor as 
derived in [110] or and reduces to an electric body force in the 


material and spatial representations 


Bio = Div (PF) = QF TE > bist = Pa®: (3.74) 


The chemical contribution PŪ originates from the dependency of concen- 
tration on deformation in the chemical energy (Jca = c2). Its equivalent 


Cauchy stress 


is purely hydrostatic. It corresponds to the osmotic pressure 
C 1 Cc 
p = gt (oc), (3.76) 


which is typically introduced in the setting of a membrane permeable 
for a solution but not for a solvent. The jump in concentration across 
the membrane yields diffusion of solvent through it to minimize chemical 


energy. Diffusion through the membrane stops when a critical pressure 


2 McMeeking and Landis formulate the mechanical response with Cauchy stresses 
1 T 
o=+PF*. 
J 
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across the membrane is reached that balances the chemical driving force 
[69]. An ideal solution with fa = 1 and ci = 0 gives the familiar result 
for an ideal gas p = uae R©c. [113]. In contrast to the case of fluids, 
we formulate Eq. (3.69) such that a stress free state with a" (ce!) = 0 
holds for the reference concentration, thereby avoiding further complexity 
with osmotic pressure having to be balanced in the reference configuration 
by a mechanical response. 
Now we evaluate the chemical potential with Eq. (3.44) through 


a Re 
n = See = pd + ROM (seer) 


bea = pl + pM with E o, K 
Ha = ack = Jaa? ap 


(3.77) 


The mechanical part of the chemical potential therefore depends on the 
swelling assumption, see Tab. We continue our analysis with the 
underlying assumption of compressibility which yields AB = 1. The chem- 
ical part € is in alignment with standard chemical potentials with an 
energetic and an entropic contribution [53] [12]. The mechanical contri- 
bution uM due to the elastic part of the pressure p = — str (4PMFT) 
leads to stress assisted diffusion [114]. The factor Je is a consequence 
of large strain considerations. Concentration-dependent properties of 
elasticity and electrostatics would yield additional contributions to the 
chemical potential [64]. 

Finally we evaluate the electric field from Eq. (3.44) and get 


= — (e)! D, (3.78) 
which gives together with Eq. (3.11) the Poisson equation 


—Div (ene; Grad (®)) = pĝ. (3.79) 
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The description in the material configuration of Eq. (3.79) aligns with 


the well known spatial Poisson equation, see Eq. (3.14). 


3.4.4 Multicomponent Mass Flux in Finite Strain 


Thermodynamically consistent transport models have to obey cross- 
relations between driving force/mass flux couples. The Onsager theory 
[LI5}HIT7), the Maxwell-Stefan formalism [102] [103] or the Poisson-Nernst- 
Planck equation are different approaches to satisfy thermodynamic con- 
sistency. In our formulation this is equivalent to postulate, motivated 
by the requirement of the reduced dissipation inequality (3.49), a linear 
ansatz for the mass flux [68] [69] 118] [119] as 


N 
Ja = — Y M2 ¿Grad (wp), (3.80) 
B=1 


where we introduce a mobility matrix Mis» which we require to be 
positive semi-definite to guarantee non-negativity of Eq. (3.49). The 


mobility matrix 


o o 
11 12 
2B = 31 22 oe (3.81) 
aß 
with 
(a. Mi as) >20 Vag,ag E R? (3.82) 


consists of N x N entries where each entry can be again a R?*® tensor to 
represent anisotropic properties. A deeper discussion of the implementa- 
tion of anisotropy within this theory can be found in Section The 
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mobility matrix describes the influence of the driving force Grad (wg) 
on the mass flux Ja. A low mobility is a result of high internal fric- 
tion with respect to the host (on-diagonal terms) and with respect to 
component-component interaction (off-diagonal terms). They are related 
to drag coefficients and binary diffusion coefficients in the Maxwell-Stefan- 
diffusion model as discussed in Section [3.4.6 

We would like Eq. to be consistent with the spatial relation 


N 


jo = — X Mag grad (wa), (3.83) 
B=1 


correlating the mass flux on the deformed geometry with a spatial gradient 
of the electrochemical potential and obtain [P] 


og = IF Map F". (3.84) 


The property Mag is the spatial mobility matrix, a parameter closely 
related to measurable quantities and Eq. (3.84) effectively informs the 
material configuration of changes of transport paths due to deformation. 
Eventually, we solve the material mass flux in Eq. (3.80) with the material 
mobility of Eq. (3.84) which is equivalent to the desired spatial mass 
flux of Eq. (3.83). If the domain is deformed, the material mobility M2, B 
becomes anisotropic due to the conversion of Eq. (3.84), even if isotropic 
material properties are assumed. This anisotropy is solely related to 
deformation of transport paths. The requirement that the material 
mobility matrix has to obey a positive semi-definite structure leads along 
with Eq. (3.84) to the spatial mobility matrix having to be positive 
semi-definite. This property can be checked in the manner of Eq. (3.82) 


3 We deduce from j, - nda = Ja - NdA the relation Ja = JF~'j, [98]. Further, 
the material and spatial gradient operator are linked via Grad (-) = oe grad (-) = 


Ferad (-). 
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material frame spatial frame 
F 


> 
U R 
a Zen 
M? 5 Mag 


Mis 


Figure 3.6: Mis is defined on an intermediate configuration, where the orientation is 
linked to the material frame but the stretches and hence the length of the effective 
transport paths are linked to the spatial frame. 


or by the eigenvalues which are required to be non-negative. Furthermore, 


we can conclude that the determinant of Mag has a non-negative value. 


3.4.5 Modeling Anisotropy in Solid Electrolytes 


In order to include anisotropic material properties as found e.g. in laminar 
structures [37] [57] in our transport framework, we link transport paths 
with the material configuration and use properties like the effective length 
of a transport path from the deformed configuration. For illustration 
purposes, imagine a set of parallel lines on a material which are associated 
with the material points. Their shape and length will change with 
deformation. Now consider an ion traveling through the material without 
being able to cross the lines. The ion then has to follow the material 
structure (the lines are borders) but needs to move the spatial length of 
these lines (deformed geometry). To implement this in our framework, 
we decompose the deformation gradient into F = RU with a right 
stretch tensor U and a rotation tensor R, see Eq. (2.5). This defines an 


intermediate configuration, where the material is stretched with U but 
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not rotated. It serves for the definition of the transport properties, e.g. 
anisotropic mobility matrix Mie: If the structure rotates, the effective 


spatial mobility Mag does so via 

Mas = RM} R". (3.85) 
We get with Eq. (3.84) the material mobility matrix as 

= U "MUT". (3.86) 


The material mobility matrix is ultimately influenced only by stretches. 
The mass flux stated in the material frame takes length changes of the 
transport path into account, the orientation, however, is governed by a 
path associated with the material. 

If an electrolyte is isotropic, we can reduce the formulation by the 


introduction of an isotropic mobility Moe = MagJ, which results in 


-1 
Mas = Map RRT =Magl and M2, =J (FTF) Mas (3.87) 


Therefore, the spatial mobility matrix coincides with the intermediate 
definition. Nevertheless, the material mobility behaves anisotropic due to 
deformation. This anisotropy is required because of the nature of finite 
deformation where the partial differential equations are solved in the 


undeformed configurations. 


3.4.6 Relation to Maxwell-Stefan Theory 


From Eq. (3.83), we can directly deduce the Maxwell-Stefan diffusion 
model [103], which is also the foundation of Newman’s derivation [53] 
120]. For simplicity, we restrict ourselves to isotropic mobilities 
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Mag = Magl. Recalling Eq. (3.83), inserting Eq. (2.15) and multiplying 
both sides with caM_; yields 


N N 
cagrad (we) = Y Bag (va — vo) = Y Kag (va - vB), (3.88) 
a=1 a=0 
with the definition 
-1 
My Mi> 
Kop = —cgM 3 5Ca with M35 = | Ma Moa >> (3.89) 


aß 


where we identify Kag as the drag coefficient as introduced e.g. in [53]. 


The transformations between the two drag coefficient matrices are 


N 
> Iran Kay a= 


Kag = (3.90) 
Kap a#ß 
0 a=ß 
= aay ax p= 
Kag = E E (3.91) 
zu NE, Kys Bfa= 
Kag else 


Note that Kag € RN+!*%*? considers friction with the host but no self 
or background friction. The matrices Kap e RX*N and Mag € RNXxN 
have non-zero diagonal elements. These diagonal elements correspond 
to interactions with the underlying media. All reformulations maintain 


the symmetry of the matrices. The transformation of drag coefficients to 
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mobilities is Mag = CacpKag- The Maxwell-Stefan diffusion coefficients 


Dog are related to the second drag coefficients with 


Kag = Ro 2 


3.92 
A (3.92) 


where cr = DN 


ao Ca 18 the total concentration. Further, a host con- 


centration co is necessary within this description. For dilute solutions 
with no component-component interaction (Mag =0 Va # 8), we ob- 
tain the correlation Maa = po- Zop V8 =1.N. The Maxwell-Stefan 
diffusion coefficients are similar to self-diffusion coefficients obtained by 
experiments [57] [20] or MD simulations [121]. If cross dependencies are 
present, the conversion within the concepts of mobilities and Maxwell 
Stefan coefficients is not as trivial and we refer to Section [3.5.5] for a 


binary electrolyte. 


3.4.7 lonic Current and Mass Fluxes 


The driving force for mass transport in Eq. (3.80) is partially due to 
the chemical potential and partially driven by gradients of the electric 
potential, see Eq. (3.50). We now split the mass flux into a chemical and 


an electrical contribution with 


JS =-—Y¿_¡M?sGrad 
Ja=J%+JE with ll Ya-ı Mas (ne) (3.93) 
JE = -r Grad (6), 
where we have introduced the mass conductivity 
N 
k = FO M2928, (3.94) 
B=1 
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which defines the impact of the electric field on mass transport[?] A 
similar split for the electric current J, = En FzaJa from Eq. (3.10) 


can be obtained by a rearrangement of the summation terms 


JO =- «°Grad 5 
Jı =J +JẸ with a Dont a Grad (Ha) (3.95) 
J? =—k°Grad (8). 


We can identify the mass conductivity of Eq. (3.94) in the chemical part 


of the current density and introduce the ionic conductivity 


N 
kK? = X Dan, (3.96) 

a=1 
which relates analog to Ohm’s law (Eq. (3.15)) the electric part of the 
current density with the electric field. Together with Eq. (3.94) we 
observe that the ionic conductivity K? is directly related to the mobilities. 


Further, we introduce a so-called transference number t, with 


J? =FZ JE => Er’=Fon e E = Fak? (K?) t, 
(3.97) 


as also introduced e.g. in [122] [23] for binary systems. Physically, the 
transference number t, represents the fraction of the overall current 
transported by species a in case of pure migration with no diffusion 
mechanism present (Grad (ua) = 0) [124]. The transference number 
in Eq. (3.97) is therefore a representation of mobilities and thus a 
well defined material property. Note that transference numbers are 
in general anisotropic, their sum is the identity tensor om t =I) 


and, based on our derivation, we find that there is no thermodynamic 


4 In the case of uncharged species—as found e.g. in active particles [63]—the conduc- 
tivities from Eq. (3.93) are zero. Consequently, the electric portion JE drops out 
and only mass transport due to chemo-mechanics remains. 
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restriction demanding positivity of the transference number. In fact, 
negative entries in tí, are possible when one component is more strongly 
affected by the chemical potential of another species, i.e. the component 
under consideration is dragged along by one or more other components. 
Recent experimental work confirms this theoretical result. We 
also refer to Section [3.5.6] which gives further insights on how to test 
thermodynamic consistency of the transport parameters. 

In the following we assume that at least one species is charged (zu # 0). 
In such systems, the conductivity and transference number are measurable 
characteristic parameters and we can rewrite the current density of Eq. 
(3.95) with the last term in Eq. (3.97) as 


N 

1 

Jg=—K > mz ta Grad (Ha) — #°Grad (®). (3.98) 
a=1 2 


Next, we rewrite the electric part of the mass flux of Eq. (3.98) as 


1 

E Cc 

Jo = pz ta (u); (3.99) 
which gives us the possibility of decomposing the species mass flux into 


three components 


N 
1 
Ja = (I-t) JEt? Pa EJ 3.100 
( a) Q a Da ta Etpe q ( ) 


namely, a self diffusion, a drag diffusion and a migration portion, re- 
spectively. Note that Eq. (3.98) and Eq. (3.100) are generic versions of 


Newman type equations for transport in concentrated solutions [53]. 
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Care has to be taken when evaluating the gradient of the chemical 
potential of Eq. (3.77). With cq EI ce /y (m &/(J.Js) due to the 


effect of compression this gives rise to 


Cc 
Grad (pa) = ES On Grad (c2) + €?Q,Grad (pM) + Ele Grad (Je) 


(3.101) 
with the prefactors 
Bol aln Js 
fa = J (1- en) Í 
& = Jep, (3.102) 
o ð C 
Je — sQ M_ Co pd 
Es = Pap 1-3 dc. 


The first term on the the right hand side of Eq. drives classical 
diffusivity, the center term is the mechanical driving force due to pressure 
gradients, whereas the last term specifies a second order effect due to 
volume change. In other words, if the volume changes due to elastic 
deformation, the ions are packed closer and the chemical energy increases. 


If the elastic deformation is inhomogeneous, the third driving force in 
Eq. (3.101) originates. 


3.5 Binary Electrolyte 


Ionic charge in binary electrolyte systems is carried by a positively charged 
cation and a negatively charged anion which we denote by a = {+, —}. 
As we consider Lithium based binary electrolytes, we apply z+ = 1 and 
z_ = —1 henceforth. In solid electrolytes, these ions are embedded in a 


solid host material, whose kinematics are described by the deformation 
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map e; this is fundamentally different from liquids where a third mobile 
component, the solvent, is present. Therefore, we have only two mobile 


species and obtain the binary mobility matrix 


M, Mi 
{Magla, 6 = {+,-}} = ls (3.103) 
M 


with three independent mobilities (which can each be functions of con- 
centration, stress, etc.). The diagonal contains the cation mobility Mı 


and the anion mobility M_. The off-diagonal term M+ in addition to 


the fugacity stems from the non-ideality of concentrated solutions and 
can be associated with the motion of clusters of both cations and anions 
with unlike charge that cause the constituting cations to be influenced 
by the anion potential and vice versa. In the settings of ion-exchange 
polymer-metal composites (IPMC) for sensing and actuation, this 
effect is known as electro-osmosis drag where water is shuttled along with 
hydrogen [126] [127]. Another motivation for the off-diagonal term is a 
counting argument. There are four constitutive quantities, namely the 
chemical potential and the three mobilities. This is the same number of 
parameters as for the standard parametrization of binary ion transport 
[53], see Section [3.5.5] 

In the following, we will present three approaches to compute the 
motion of the mobile species {c+,c_,pg}. The first is a straight forward 
application of the previously derived theory for binary electrolytes. The 
second approach is a reformulation introducing the balance of charge 
and serves as the transition to the third, the electroneutral version that 
constitutes simplifications of the first two. All of these formulations 
describe transport and, hence, the underlying structure of the mechanics 
with the Piola stress tensor P from Eq. (3.70), the material mobilities 
Mis from Eq. (3.84) and swelling from Eq. (3.61), is unaffected and 


remains the same in all versions. 
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3.5.1 Full Formulation 


For the full formulation in terms of (py, cf.,c&, Dj, we state the balance 
of linear and angular momentum, add mass conservation for positive and 


negative ions derived from Eq. (3.83) and include Gauss’ law as 


& + Div (Ji) =0 
J, = -M$ Grad (w+) — Mí Grad (w_) 
& +Div(J_)=0 


(3.104) 
J_ = -M$ Grad (w;) — M? Grad (w_) 
Div(P)=0, PF” = FPT 
Div (D) = pĝ, D = -e0e,Grad (®) 
where the electrochemical potentials are given by 
w4 = u4 + F0 and w =p — FO. (3.105) 


The equations in formulation (3.104) constitute the system of field equa- 
tions necessary for the solution for the fields {4+, c$., c2, ®}. The ini- 
tial boundary value problem requires initial conditions for cf ,c? and 


boundary conditions (BC) for ı,c%,c” and ®. Typical electrochemical 


boundary values (BV) are those for blocking electrodes given by 
-N:Jılos, = -N -J-log, =0, Flas, = ©, (3.106) 
where © is a prescribed, position dependent electrical potential on the 


boundary. Also typical are galvanostatic BC specified by 


Japp 


-N s J4ləBo = “PF? -N i J_|os, = 0, (3.107) 


63 


3 Theory of Ion Transport in a Binary Solid Electrolyte 


where a boundary current per undeformed area J*PP is applied and only 
cations leave the system. The boundary current is usually given by 
reaction kinetics such as the Butler-Volmer equation [128]. Note that 
for global electroneutrality, it is necessary that f. Bo J’PP = 0, where 
Bo denotes the boundary of the electrolyte. To within an additive 
constant the electric potential is fully described due to its coupling to 


ion transport. 


3.5.2 Towards Balance of Charge 
To reformulate Eq. (3.104), we recall the cation mass flux from Eq. (3.100) 
and use its binary form as 
t 
Ty =t- J +1,.J_+ za 


= — (t? MS, + t? M2) Grad (u+) (3.108) 


te 
— (EM? +t? M2.) Grad (u-) + Ia, 


where we identify an ambipolar mobility matrix as mentioned e.g. in 
119| as 


Mamb = EM} + MS =tM? +42 MS 
2 o o o 2 o—1 (3.109) 
= P? (M}M? - M3?) 6°77. 


Note that in contrast to the treatment in [53], all concentration depen- 
dencies are hidden in the mobilities. Further, we exchange the anion 
concentration c? with the field variable total charge density py via 


C 4 and obtain from Eq. (3.98) a statement for the charge flux 


o 


J,= -5 [t Grad (u4) — t? Grad (u_)] — k*Grad(B). (3.110) 
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From Eq. (3.104) we deduce with 


+O “oO 
Pq = Fac, 4 


H Fz ¿e 
(3.111) 


= Div (Fz} J4) + Div (Fz_J_) = Div (J4), 


H+ 


the electric potential ® = = — % and some rearrangements the second 


set of transport equations 


& + Div (J4) 
J4 


pe + Div (J4) 


=0 

© t+ 
= —M emp Grad (u+ T p-) + Fu 
=0 


= E (3 -MGrad(u+n))| (8112) 


=0, PF? = FPT 


where the current density is driven by the electrochemical potential of 
the cation, see Eq. (3.105). With respect to the BC, we exchange the 
anion flux in Eq. (3.107) with a boundary current with 


-N -Ja = J? (3.113) 


consistent with the cation flux of Eq. (3.107). 


3.5.3 Electroneutral Description 


Ion polarization occurs only in the double layer close to interfaces and 


results in strong electric potential gradients. The width of this domain 


can be approximated b 


y the Debye length [104] and is usually of the order 
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of nanometers. Away from this region, electrostatic forces are dominant, 
drawing unlike charges together, and therefore electroneutrality prevails. 
The relationship 

SR (3.114) 
is an adequate approximation for simplifying the transport equations and 
the individual concentrations c and c? can be both replaced by a salt 


o 


concentration cf = c°. Consequently, the charge density becomes zero 


pg = 0 everywhere and Poisson's law decouples from mass transport [116]. 
Besides the replacement of c$. and c? by c?, it is beneficial to replace ® 


by the electrochemical potential of the cations 
We = H+ + Fz,0, (3.115) 


because it is directly accessible in probe measurements [53]. Consequently, 
we reduce the set of field variables {y;,c§,c2,®} to 
{y:,c°,w+}. Conservation of charge, as given in the form of an ex- 
tended Ohm’s law in Eq. (3.112), then serves as the governing equation 
for the electrochemical potential. 

Besides eliminating a field variable and thus saving computational effort, 
another advantage of the electroneutrality assumption is the reduction in 
the number of required parameters. This is significant especially since 


individual ion properties are very difficult to measure. With the molar 
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volume of the salt Q4 


L = Q} + Q, the chemical potential of the salt 


ug := uf + pE and considering Eq. (3.101), we obtain 
& + Div (J) =0 
t? 
J, = -D°Grad (c) + Y 
— EPQ+G°Grad (pM) — ¿t G°Grad (Je) 
Div (J,) =0 (3.116) 
KÊ 
J, = D¿Grad (c”)- 7, Grad (w+) 
+ €PNy G Grad (p™) + &£G7Grad (Je) 
Div(P)=0, PF” = FPT 


with contributions of diffusion, migration, stress-assisted diffusion and 


volumetric effects. The prefactors are given by 


uf (c) 
o eny? Æ 
D & amb dc > 
o apo _ 1 (Ou£lc) 
G” = amb — = ( dc 
& Oz (c) 
D, — F (I t$.) dc °, 
1 


(3.117) 


with the Fickian diffusion coefficient D° formulated in the material setting. 


- to be of the 


Choosing the chemical potential ¡us 


mean salt fugacity fy 


obtain 


_ ln fs 


( 


form of Eq. (3.77) with 


(f-)(f+) as introduced in Eq. (3.69), we 


(1 


Olnc 


(3.118) 


2), 
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where the expression in brackets is usually called the thermodynamic 
factor. 

The set thereby constitutes an extended version of the classical 
Newman theory [53] [129] with which it coincides for the case of p, = X 


(no deformation) and 24 = 0 (no swelling). Under these conditions the 


last two terms on the right-hand side of the mass and current flux in 


Eq. drop out. 


The osmotic stress in Eq. (3.75) becomes 


ex 1 E 
p° =2RO (1 p? nt) de, (3.119) 


cref 


and thus is solely depended on the thermodynamic factor. 
Other than the parameters associated with swelling in Eq. (3.61) due 


to the salt molar volume, i.e. Q4 and the reference salt concentration 


cf no further properties than the standard electrochemical parameters 


{D°,«°,t§,14+ om le) are required to couple mechanics and ion 


transport in this formulation. 


3.5.4 Steady State Current in a Binary Electrolyte 


Charging and discharging a battery involves cation shuttling through 
the solid electrolyte. During this process, a concentration gradient of 
the salt builds up, see discussion in Section [4.4] When the steady state 
is reached, the concentration gradient has the highest slope eventually 
determining a steady state conductivity. The steady state conductivity 
property is important to judge the performance of an electrolyte and is 
also easy to access in measurements. 

Since the anions are restricted to the electrolyte domain and are, unlike 
the cations, not inserted and extracted during operation, their dynamics 
determine the steady state. In other words, when the driving forces 


for anions are in equilibrium and their mass flux vanishes, i.e. J_ = 0, 
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then the system is in steady state. From the mass fluxes in formulation 
(3.104), where we eliminate w_ in the expression for J}, we obtain with 
(Ja = FJ, an expression of a steady state Ohm’s law as 


(3% = - (193% Grad (=) (3.120) 


with the steady state conductivity 


(r°) = F? (MS, — M3 (M?) M2). (3.121) 


Expression (3.121) shows the importance of off-diagonal mobilities. The 
steady state conductivity would only be governed by the cation mobility if 


MÍ, is assumed to be zero. It is worth noting that effects from mechanics 


and the chemical potential do not enter the calculation explicitly. In the 
case of an isotropic system (Mag = Mag), we can deduce requirements 
on the entries of the binary mobility matrix from Eq. (3.103) from a 
positive steady state conductivity. To do so, we write the steady state 


conductivity in the spatial form as 


(Ki) = iF (K?) FT = - (M,M_—M2) I. (3.122) 
The expression in brackets is the determinant of the spatial mobility 
matrix. Since both the determinant and the steady state conductivity 
are non-negative, we conclude M_ > 0. Due to the structure of the 
determinant, one can further conclude that My has to be non-negative, 
too. 

Alternatively, on can follow Balsara and Newman [130] and use the 
formulation to derive the steady state conductivity. It is reached 


when 


t_ 


J- = -MampCrad (44 + H-) = 7 


JSS Ło. (3.123) 
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and gives a correlation between the steady state current J “i and the 


chemical potential. Rearrangement yields 


1 Ke 
(3,5 = ae (e 7, Grad (ur + 1) E (3.124) 


where we have introduced the so called Newman number 


iso ,.iso iso+iso? 
_ er)? _ (M- Me) o Oz (19) 
Fi det Map det Mas F2 piso Oc . a 


Ne 


Here we have used Eq. for MS, and assumed isotropic transport 
properties. Eq. (3.125) is in alignment with to work of Balsara and 
Newman [130]. The Newman number is non-negative due to the ther- 
modynamic constraint det Mag > 0. The right hand side of Eq. (3.124) 
corresponds to the diffusion part of the electric current in Eq. (3.112). 


Consequently, we obtain 


Ko 


== (3,12 
1+ Ne (84120) 


(IL = - (193% Grad (=) with (K9) 
which is consistent with the result of Eq. (3.120). 


3.5.5 Conversion of Transport Parameters 


Binary electrochemical systems are usually characterized by macroscopic 
measurements yielding isotropic transport parameter {rl tis°, Diso 1 + 
2l fel) We now want to correlate these parameters to the mobilities 


and drag coefficients introduced in Section and Section as 
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these give an idea of the interaction of the anions and cations in the solid 


electrolyte. With Eq. (3.96), Eq. (3.97) and Eq. (3.117) we obtain 


piso ue My = M+ 
+  M,+M_ -2M4 (3.127) 
M;M_-M? 2RO& a 
pe = + ote ROEL 1 + 0 In f+ (c) 
M,+M_-2My c Olnc 


The spatial properties transform into the material description with 


=1 ; 3 =l 
Parr) Baw, DIr D 
(3.128) 


Taking the inverse of Eq. (3.127), we deduce mobilities from measured 


parameters as 


2 


Dc . tise 
M, = i 150 Y y= {+ =} 
c Oln f+ 2l fale) F2’ i 
Die a oe i 
M+ = c aln f+ (c) i F2? ` 
2ROE (1+ 2,0) 


The first part in Eq. describes the influence of diffusion, the 
second equation reveals the diffusion and migration contributions to the 
mobility. The ratio of both for the anion mobility M_ is given by the 
Newman number, see Eq. (3.125). 


Further, if we assume f+ = f+ = f-, we can compute from the fugacity 


the thermodynamic factor by integration as 


i “ Əl fe(c) 1 
fr lc ) = exp ( : Oae) R (3.130) 
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Eq. (3.129) and Eq. (3.130) therefore yield the transformation rules to 
convert transport properties used in formulation (3.116) to the properties 


used in Eq. (3.104). Using Eq. (3.90) with the inverse of Eq. (3.89) gives 
the correlation between binary drag and diffusion coefficients for a binary 


system as 


M+ M+ — 1 Ko + K K 
M+ M_ BA (Ko+Ko- + Kt (Kor T Ko_)) K+ Ko+ + K+ 


(3.131) 


The on-diagonal mobilities Maa stem from a complex interplay between 
friction of host as well as the adjacent component. If the drag coefficients 
are positive then the mobilities are also positive. Combining Eq. (3.127), 
Eq. (3.131) and Eq. (3.92) provides us with 


iso _ F? (Ko+ + Ko_) 
c? (Ko+Ko- + Ks (Ko+ + Ko-)) 
Ko- Do+ 
he = = 3.132 
+ Kor + Ko Dor + Do ( ) 
io E ere Doy 2o- cre 


=! D 
amb — Koy + Ko- ~ ROC Doy + 2o- ROco 


where we have linked the isotropic ambipolar mobility M'S°, from 


Eq. (3.109) and the thermodynamic diffusivity Y = oe un [53] [120]. 


It is worth noting that only the conductivity depends on ie component- 


component friction K+. Further, the transference number can only be 


negative if one of the drag coefficients Kg_ or Kg, is negative. 

It is important to note that measurements of electrochemical transport 
parameters do not usually take the mechanical response into account. As 
seen in the derivations, conductivity and transference number are not 
influenced by mechanics. The steady state conductivity from Eq. (3.122) 


is another measurable parameter without direct impact of mechanics. 
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The diffusion coefficient is, however, usually identified by experiments 
where quite large concentration gradients occur. In such measurements, 
we expect a strong influence of mechanical driving forces which assist the 
diffusion. The measured relaxation time can be then only associated to 
an effective diffusion coefficient. An alternative approach to gain diffusion 
coefficients without contributions of stress driven diffusions might be 
atomistic simulations [131]. 

Finally, we want to point out that a third component would introduce 
three new transport parameters due to the non-diagonal entries. In 
general, having N species requires N(N+1)/2+1 independent measurable 
parameters to account for all possible interactions between ions and the 


chemical potential. 


3.5.6 Thermodynamically Consistent Parameters 


It is important to use thermodynamically consistent parameters in order 
not to violate the fundamental second law of thermodynamics. Obviously, 
this has to hold for both the formulation with mobilities (Eq. (3.104)) 
and the parameters of the concentrated solution theory (Eq. (3.116)). 
For the former, we have concluded that the mobility matrix Mag has 
to be positive semi-definite, see Eq. (3.81). The eigenvalues, therefore, 
can be checked for validation. In binary isotropic electrolytes, an easier 
assessment is possible. From Eq. (3.109), and (3.118) we deduce 
that the determinant is positive for positive conductivity, diffusivity and 


thermodynamic factor. In isotropic systems, this reads 


— piso „iso c f On f(c) = 
det (Mag) = Di°k A (1. ae ) . (3.133) 


Note that a positive determinant does not necessarily yield a positive 
definite matrix. However, we know that the diagonal elements are positive 
due to positive steady state conductivity in Eq. (3.122). If the symmetric 
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mobility matrix has a strongly diagonally dominant structure, it has a 
positive definite form. This feature is given with 0 < t§° < 1 as seen in 


Eq. (3.127). Therefore, we can conclude that 


be OS)”, oops 
| dln f(c) (3.134) 
438° 1 1 TN 

O<t <1, Doda Jire 


is sufficient to fulfill the second law of thermodynamics. If not, the 
eigenvalues of Mag have to be computed and checked for non-negativity. 
Therefore, the transference number tio cannot just independently take 


any real value as proposed by Monroe and Delacourt [102]. 


3.6 Summary 


In this chapter, we have presented a rigorous formulation of ion transport 
in a solid electrolyte. Using the concept of finite strain, the Coleman- 
Noll procedure and the Helmholtz energy, we derived in a consistent 
manner the driving forces for species transport, the mechanical response 
of the host material and their coupling mechanisms. We thereby chose 
the transport description in the deformed geometry but solve it —as 
common in continuum mechanics— in a reference description. The 
influence of mechanics turned out to be broad. Most notable is the 
consideration of ion size by means of swelling and eventually stress 
driven diffusion. Further, we considered the effect of the change of 
geometry on mass flux and in the chemical potential. The former was 
also discussed for an anisotropic material and the latter gives rise to 
an osmotic pressure and is closely related to the thermodynamic factor. 
Using electroneutrality assumptions for a binary electrolyte, we have 


derived an extended version of concentrated solution theory [53] taking 
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diverse mechanical coupling mechanisms into account. We identified the 
connections to the Maxwell-Stefan theory and presented the correlation 
of mobilities, drag coefficients and transport parameters by means of 
conductivity, diffusivity, transference number and thermodynamic factor 
from concentrated solution theory. Further, the relation to the steady 
state current was shown and thermodynamic consistency of transport 
parameters was discussed. The off diagonal terms of the mobility matrix 


and hence the anion-cation interaction were found to be significant. 
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In this chapter, we complement the theoretical description of the coupled 
electro-chemo-mechanical modeling of a binary solid electrolyte with 
numerical studies as shown in Ganser et al. [101]. We utilize simulations 
to highlight the coupling mechanisms between electrochemistry and 
mechanics in solid electrolytes and some of their implications on battery 
performance. Whenever possible and meaningful, we utilize full sets of 
measured, concentration-dependent transport properties. Three studies 
with polymer electrolytes and thus mobile anions and cations are shown. 
In particular, we examine the influence of the changes of the mechanical 
material properties on the response of the electrolyte, while keeping the 


electrochemical properties fixed. 


4.1 Simulation Environment 


We utilize the finite element solver COMSOL Multiphysics® [132]. Pois- 
son equation (from Eq. (3.104)), conservation of mass (from Eq. (3.104) 
and Eq. (3.116), respectively) and balance of charge (from Eq. (3.116) 
were implemented as partial differential equations in the "general form" 
and discretized with quadratic Lagrange elements. Since no theory de- 
velopment was done in the mechanical part of the transport theory in 
Chapter[3] we were able to use the COMSOL modules "structural mechan- 


ics” and “nonlinear structural materials” to solve the deformation field 
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with cubic Lagrange elements. All equation were solved in a fully cou- 
pled manner and all field variables involved were normalized to optimize 


accuracy and convergence. 


4.2 Material Data 


The most common polymer electrolyte is poly(ethylene oxide) (PEO) 
which shows sufficient transport and chemical stability. However, PEO 
creeps rapidly at elevated temperatures and its mechanical response then 
comes closer to that of a liquid than to that of a solid. Composite systems 
such as the block copolymers PS-PEO are under investigation to combine 
the transport property of PEO with enhanced mechanical properties 
of polystyrene (PS). However, besides improving mechanical properties, 
such compositing strategies generally reduce the volume fraction of the 
conducting phase and thus lead to reduced transport properties. This 
means that improving, e.g., mechanical stiffness generally negatively 
impacts electrochemical transport. In the following studies where the 
influence of mechanics is investigated, we keep the electrochemical param- 
eters constant and thus do not account for these competing effects as we 
are interested in highlighting the effect of mechanics on electrochemistry. 


This is easier to study if electrochemical properties remain unchanged. 


4.2.1 Electrochemical Transport Parameters 


Only very few complete electrochemical characterizations are reported for 
polymer electrolytes. We utilize a set of concentration-dependent param- 
eters for sodium salt (NaCF3SO3) in a PEO polymer electrolyte and 
a recently published work [56] with lithium salt (LiTFSi). Furthermore, 
we take data from Timachova et al. for the PS-PEO polymer with the 


lithium salt LiTFSi. In all cases, we assume a homogeneous material and 
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Figure 4.1: Transport properties data of PEO with sodium salt [55], with lithium 
salt and PS-PEO with lithium salt [57]: (a) ionic conductivity, (b) diffusivity, (c) 
transference number and (d) thermodynamic factor. 


define the concentration with respect to the overall PEO and PS-PEO 
volume, respectively. 

The concentration-dependent transport parameters are shown in 
Fig. In general, we can observe that the conductivity of PS-PEO 
is almost one order of magnitude lower compared to that of the two 
PEO electrolytes. The diffusivity of the sodium system varies almost 
one order of magnitude depending on the concentration. In contrast 
to that, PEO with lithium salt has an almost constant diffusivity of 
Di 21x 10-1 =. which is much higher than the PS-PEO system. 
The transference number as shown in Fig. is strongly negative for 
high concentrations of sodium salt while the other two lithium systems 
show very similar response, where only ira of PEO with LiTFSi is 
slightly negative for higher salt concentrations. The strong deviation of 


the sodium system compared to the lithium systems might be either a 
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Figure 4.2: The graphs in the first row show the converted mobilities (a), drag 
coefficients (b) and chemical stresses (c) from PEO polymer electrolyte with sodium 
salt [55]. The graphs in the second row show the converted mobilities (d), drag 
coefficients (e) and osmotic pressure (f) converted from PS-PEO block-copolymer 
electrolyte [57]. Fig. (f) shows also the osmotic pressure from PEO with lithium salt 
196]. 


consequence of the different material or simply due to a different evalua- 
tion method: Both Pesko et al. [56] and Timachova et al. [57] utilize the 
recently introduced concept of the Newman number, see Section [3.5.4] to 
measure transport properties. The thermodynamic factor is similar for 
the lithium salt systems and increases almost linearly with concentration. 
The sodium system on the other hand shows a very small thermodynamic 
factor. 

We now take the measured transport properties of Fig. [4-T]and convert 
them into mobilities and drag coefficients. We do this for PEO with 
NaCF3503 and for PS-PEO with LiTFSi salt. The mobilities of the 
former are shown in Fig. and were computed with Eq. (3.129). 
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It can be observed that anions are slightly more mobile than cations 
throughout the concentration range. Strong cross interactions between 
the ions are present due to the high off-diagonal mobility M+. This is 
confirmed by the drag coefficient K+ in Fig. [4.2b obtained with Eq. 
(3.90). The negative drag coefficient K_o is a consequence of the negative 


transference number. Although the negative sign is non-intuitive, one can 
show by an eigenvalue analysis that the system is still thermodynamically 
consistent, see Section [3.5.6] 

The resulting mobilities of the PS-PEO/LiTFSi system are depicted 
in Fig. and show a low cation mobility in the concentration range of 
interest. The mobile anion has a small drag coefficient (Fig. |4.2b). The 
negativity of the off-diagonal mobility goes along with a negative drag 
coefficient between anion and cation. Note that this negative mobility 
is not related to a negative transference number. Again, these data are 


admissible from a thermodynamic point of view, see Eq. (3.134). 


4.2.2 Mechanical Properties 


Measurements of mechanical properties of polymer electrolytes are chal- 
lenging due to several reasons. The electrolyte layers used in batteries 
are very thin, the operating temperature range is usually elevated and an 
inert atmosphere is required. The few available characterization access 
mostly the elastic properties, although mechanically irreversible and time- 
dependent behavior is strongly expected. In the following, we will leave 
the Young’s moduli as design parameters and disregard inelastic effects 
due to a lack of information. The orders of magnitudes are as follows: 
PEO provides a stiffness of roughly E = 10 MPa at room temperature 
[133] [34]. However, its ability to resist deformation at elevated tempera- 
tures decreases rapidly [48]. The mechanical properties for PS-PEO are 
dominated by the PS-phase, leading to roughly E = 150 MPa [45] [48] at 


operation temperatures of 80°C. We assume a Poisson ratio of v = 0.4 
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for all electrolytes as found for e.g. polystyrene [135] and an isotropic 
response throughout the work. 

The reference value of the molar volume is Q4 = 140 m for LiTFSi 
[136] and Q4 = 80 ~ for NaCF3SO3 [137]. Note that the ions are 
assumed to be compressible in the swelling model used in this work, 
see Section The reference salt concentration for which swelling is 
assumed to be zero is set to be cf = 1000 me for PEO (both salts) and 
cf — 1500 me for PS-PEO. Fig. and f shows the osmotic pressure 


of Eq. (3.75) as a function of salt concentration. With values up to 2 MPa 


for PEO with sodium salt it is small but not negligible. The chemical 
stress in a PS-PEO electrolyte is an order of magnitude higher compared 
to PEO with sodium salt and in the same order as PEO with lithium 
salt. The similarities of the lithium salt systems come either from the 
same salt or can be traced back to the experimental assessment of the 


thermodynamic factor. 


4.3 Double Layer in a Solid Electrolyte 


In the first example we want to show the capabilities of the full model as 
given by Eq. (3.104) that treats anions and cations separately and can 
thus be applied to situations where electroneutrality cannot be assumed 
everywhere and where the ions behave differently such as in a double 
layer. This effect is characterized by local anion/cation polarization and 
is most pronounced in the case of blocking electrodes. Related studies are 
found for liquid electrolyte systems [104] [138], solid electrolytes with fixed 
anions and references cited therein. Double layer properties can be 


obtained experimentally e.g. via impedance spectroscopy measurements 


[140]. 
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Figure 4.3: (a): Setup of the boundary value problem for blocking electrodes solved 
with Eq. (3.104). (b) illustrates the drop of electric potential over the whole domain 
and (c) zooms into the double layer regions where the electric potential decreases and 
the chemical potential increases. In systems with a high elastic modulus, mechanics 
has an influence on the result. 


4.3.1 Material 


We will restrict our investigation to the impact of mechanics on the 
ion distributions in a solid binary electrolyte. Note that concentration- 
independent material parameters are used in this study since the avail- 
able parametrization as discussed in the previous section resolves the 
concentration dependency only with respect to the salt and does not 
distinguish between anion and cation concentrations. With xi = 0.01 5, 
Dis =5x 10-1? =, tise = 0.3 a system similar to PS-PEO is utilized. 
We assume an ideal mixture such that f+ = f- = 1, see Eq. (3.69). 
These parameters transform with Eq. to My = 1.1 x 10-12 mel 


mJs? 


M_ = 1.5 x 10-12 mE and M4 = 7.6 x 10-13 mel, showing that the 


mjs’ 


anion is more mobile than the cation. The electric permittivity is set 
to er = 10 [141], the molar volume of the salt is Q4 = 140 cm? [136]. 


mol 


Although ion sizes for some non-aqueous solvents are reported [142], there 
is a lack of such information for solid state electrolytes and therefore 
we assume Q_ = 100,. The reference concentration is assumed to be 


r 


exe = 1000 mol, As the influence of mechanics is the focus of this work, 


we vary Young’s modulus in the range 10 MPa to 1000 MPa. 
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4.3.2 Boundary Value Problem 


We consider a blocking electrode (Fig. 14.3h) of width L = 0.5 ym with 
no-flux conditions for anions and cations on the boundaries. We assume 
that the right boundary is electrically grounded (®, = 0) and on the left 
boundary subject to ®; = 10mV. The system is mechanically fixed on 
both sides. 


4.3.3 Results 


A stationary solution exhibiting a double layer is reached within 1 ms. 
Fig. |4.3b shows a constant electric potential for most of the domain. The 
potential drops occur very locally at the two electrodes within the double 
layer (Fig. [4.3F) which extends to a width of 0.5 nm, comparable to the 
results in [139] and the Debye length of Ap = 0.1 anf!] 

Fig. shows the related concentration profiles of anions and cations 
within this double layer. Away from the electrode, the concentrations are 
basically equal which confirms electroneutrality for the inner part of the 
electrolyte. In the region close to the negative electrode the concentration 
of anions decreases whereas the cations accumulate. The concentration of 
cations is such that their chemical potential gradient balances the electric 
field (Fig. [1-3p). Comparison of the results with different Young's moduli 
reveals that the anion concentration depends strongly on the mechanical 
stiffness: The higher the material stiffness, the more pronounced is the 
non-symmetric response of the ions (Fig. [4.4h). The reason for this is that 
the large size of the anions leads to swelling that cannot be countered by 
depletion of the smaller cations. Consequently, a pressure gradient builds 
up and increases the mechanical contribution to the chemical potential 
(Fig. 4.3). Although both ions are affected by this driving force, the 


1 The Debye length represents the length scale of the diffuse layer and is given by 


Ap = y/ 22 [T38]. 
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Figure 4.4: The figures show the concentration (a) and the pressure profiles (b) in the 
double layer for various material stiffnesses. The higher the elastic modulus, the more 
pronounced is the pressure gradient which acts as a mechanical driving force and 
yields a lower concentration of the anions. The impact on the smaller cations is less. 


influence on bigger ions outweighs that on smaller ones and results in 
a considerably lower increase of anion concentration at the electrode 
for high Young’s moduli. The cation concentration, however, increases 
slightly as the modulus is increased. This is a result of the lower drop of 
the anion concentration. The size of the double layer region stays nearly 


the same when the modulus is varied. 


4.3.4 Discussion 


In summary we see that the suggested full model of Eq. (3.104) is well 
suited to describe double layers and can also predict the influence of the 
mechanical properties of the electrolyte on the related electrochemical 
phenomena. Note that of course mechanical effects on the double layer 
are more pronounced in solid electrolytes than in liquid electrolytes, 
where the solvent is mobile and can relax the constraints in liquids. 
Finally we point out that the electroneutral transport model from 
Eq. (3.116) cannot capture the potential drop associated with blocking 
electrodes, since only a full model with dedicated Poisson equation and 
separate concentration fields for anions and cation such as given by 
Eq. can accurately describe the situation, specifically at the 
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boundaries. However, in contrast to the case of blocking electrodes, 
charge separation is negligible when flux boundary conditions are used 
such as for galvanostatic charge and discharge (see Example |4.4). In 
this situation, the discontinuities in electric potential adjacent to the 
electrode are features of the redox reactions occurring at the interfaces 
between the electrolyte and the electrodes. As such they are accounted 
for by models like Butler-Volmer reaction kinetics (see Chapter[5) or Tafel 
equations. As a result, the electroneutral model serves as an appropriate 
approximation when flux BC are invoked. The blocking electrode BC 
is therefore the condition that most strongly requires an approach that 


accounts for local charge separation. 


4.4 Galvanostatic Charging 


When a battery is charged, cations pass from the positive electrode (cath- 
ode) to the negative electrode (anode) through the electrolyte. When this 
process occurs at a given current density it is referred to as galvanostatic 
charging. In a binary electrolyte, both anions and cations act as carriers 
of charge, but only one kind of ion (usually the cation) enters and leaves 
the system. The other (usually the anion) accumulates against what 
are effectively blocking electrodes. The cation is locally attracted to the 
polarizing anions due to strong electrostatic forces (leading to effective 
local electroneutrality) and eventually a salt concentration gradient builds 
up. These profiles dictate the performance of a battery and are thus one 
of the classical tasks for electrochemical simulations to predict. Locally 
low concentrations are not favorable due to their influence in slowing 
reaction kinetics on the interfaces, while locally high concentrations can 
lead to phase separation, e.g. depletion of lithium. The concentration 
gradient is usually governed by the diffusion coefficient, the transference 


number and the thermodynamic factor. In this study, we emphasize 
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Figure 4.5: Setup of the boundary value problem of galvanostatic charging for (a) the 
full model, Eq. (3.104) and (b) the reduced model, Eq. (3.116). (c): A concentration 
gradient builds up during charging for E = 10 MPa. Both, model (3.104) and model 


(3.116) lead to the same result. 


that mechanical properties can also be important in controlling the salt 
concentration gradient and thus the performance of the electrolyte. We 
utilize isotropic parameters of the PEO polymer with NaCF3503 salt 
at 85°C operating temperature, see Section We will vary 
the Young’s modulus in the range of 10 MPa to 500 MPa. For the full 
description as given in Eq. (3.104), we utilize e, = 10, f} = f- and 
Q- = 1094. 


4.4.1 Boundary Value Problem 


We use both the full model without electroneutrality from Eq. (3.104) 
and the extended Newman model from Eq. (3.116) to study the complex 
interplay of electrostatics, chemistry and mechanics. The charging process 
can be treated as a 1D problem where we assume a separator thickness 
of L = 30pm. We apply a current of J?PP = 0.5 n5 and set © = 0V 
for the full and w+ = 0 V for the electroneutral model on the right, see 
Fig. Again, we fix the left and the right boundary mechanically. 
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Figure 4.6: (a): Steady state concentration gradient which occurs during charging for 
three different Young’s moduli. (b): Mechanical pressure p™, osmotic pressure pe 
and total pressure p within the electrolyte. 


4.4.2 Results - General Remarks 


Focusing first on the electrochemical response, we carry out a simulation 
with fixed Young’s modulus of E = 10MPa. From the results we see 
that a concentration gradient builds up due to migration of anions in 
the first 100s, see Fig. [4-5p. A stationary state is reached when the 
chemical driving forces on the anions equal the electric driving forces and 
hence J_ = 0. Comparison of the extended Newman model and the full 
formulation yield the same concentration and electrochemical potential 
profile within numerical accuracy. In the full formulation, only a small 
double layer develops with differences in anion and cation concentration of 
< 0.01 ma, Note that spatial and material quantities are almost identical 


in this example because mechanical deformation is small. 


4.4.3 Results from the Full Model 


To gain further understanding on the coupling mechanism, it is instructive 
to investigate the driving forces for the ions grad (wa) within the solid 


electrolyte and their chemical, mechanical and electrical components. 


Awa 
L > 


The first order approximation of the gradient, e.g. grad (wa) = 
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enables us to access the driving force by measuring Aw, = wH8ht — wett P] 
For E = 10MPa, we measure Aw, = 13.2mV for the cations and 
Aw_ = —8.2mV and hence a relatively high potential drop for the anions 
although the mass flux 7_ vanishes in steady state. This can be confirmed 


by using 7_ = 0 in Eq. (3.104) which gives 


Mz SS 
rad (w_) = => ; 4.1 
grad (ot (4.1) 
The cation mass flux || j=] = J®PP/F is constant for a 1D-problem and 
consequently, even though the anion flux vanishes in steady state (j_ = 0) 


the gradient of the anionic electrochemical potential can only be zero for 
M+ = 0. 

To understand the influence of mechanics, we repeat the simulations 
with E = 100 MPa and E = 300 MPa. We observe that concentration 


gradients are less pronounced with increasing mechanical stiffness E, see 


Fig. [4.6h. In general, the concentration polarization has two effects on 
mechanics: (1) The osmotic pressure p© increases sightly for an increase 
of salt concentration and (2) the swelling mechanism leads to a pressure 
gradient (Fig. [4.6b). The mechanical pressure pM thereby shows different 
characteristics for different moduli. It has a reversed slope with respect 
to the concentration polarization for E = 10 MPa. This is a consequence 
of the osmotic pressure p? which couples with the mechanical stress 
via Div (PM + P°) = 0. For high moduli such as Æ = 300MPa, the 
picture changes such that the sign of the pressure gradient follows the sign 


of the concentration gradient because the swelling mechanism becomes 


dominant. Local volume change is observed to be small with J = 10.025. 
Fig. shows the potential drops with respect to the mechanical 


moduli E. The electrochemical potential for the cation Aw, is composed 


2 We use the spatial description to motivate the findings and avoid therefore the formal 
material anisotropy as discussed in Eq. (3.84). 
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Figure 4.7: (a): Potential drops across the electrolyte of the individual contributions 
for various Young’s moduli computed with Eq. (8.104). (b): The cation mass flux J+ 
for different stiffnesses at X = 3L is composed of diffusion (blue), migration (red), 
stress driven diffusion (yellow) and a volume change contribution (purple), see Eq. 


En. 


of AD, Aus and Au. The chemical contributions Auf for both ions 


decrease with increasing stiffness. Note that the increased stiffness also 
leads to a less pronounced salt concentration polarization, see Fig. [4.6h. 
The mechanical contribution AM and Ay is in tandem with the me- 
chanical pressure as shown in Fig. |4.6b. However, due to the small size of 
the cation, Ape and therefore its contribution to the electrochemical po- 
tential is small. The effect on the anion is substantially more pronounced, 
see Fig. [1-7p. 

There is a small decrease of Aw, which comes from the concentration 
dependencies of the mobilities. Integrating the left equation in Eq. (4.1) 
we find that 


M_ 
P, M¿M- — M4 


Awy = l grad (w4) dv = j} dv, (4.2) 
Pi 


where the mobilities are concentration dependent. The integral can be 
interpreted as a macroscopic steady state resistance. Using the concen- 
tration profile from Fig. we obtain a smaller value from the integral 


and therefore a smaller Aw, for an increase of E. Further, we see a small 
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increase of A® for increasing E in Fig. [4-7p. To explain this behavior, we 
consider the case of constant mobilities and therefore get from Eq. (3.83) 


the estimation 


j} ~ grad (w+) = Grad (u) + Grad (u) +FGrad(®) (4.3) 
~ ~_—’ 


const. const. ~O 


The gradient in electric potential has to compensate the gradient of the 


chemical potential. 


4.4.4 Results from the Extended Newman 
Formalism 


Now, we move our focus to the solution using the extended Newman 
formalism as given in Eq. (3.116) and probe the contribution on the mass 


fluxes in the middle of the separator (X = 1L). The cation mass flux 


=3 

of —5.2 mo, is composed of diffusion, migration, stress driven diffusion 
and a contribution from the effect of volume change, see Fig. [1-7p. The 
fractions due to migration and the effect of volume change stay almost 
constant for the range of investigated moduli. The stress driven contri- 
bution gains importance for higher material stiffness and, therefore, the 
diffusion contribution is less pronounced. The latter can be seen from 
the smaller gradient of concentration in Fig. [4.6h. Note that probing 
at a fixed location with varying concentration introduces an additional 
contribution due to the concentration-dependence of the transport proper- 
ties. Positivity of the migration originates from the negative transference 
number (Fig. ). 


4.4.5 Discussion 


This example illustrates several aspects of the complex mechanisms 
occurring during galvanostatic charging. Both models, Eq. (3.104) and 


91 


4 Understanding Ion Transport in a Binary Solid Electrolyte 


Eq. (3.116), align and no pronounced double layer effects are observed. 
The parametrization shows strong interactions between anions and 
cations by means of cross mobilities M+. They have a direct influence on 
the steady state current of Eq. (3.120) and potential drop of Eq. (4.1) 
and yield Grad (w_) # 0 although J_ = 0. Young’s modulus has a great 


impact on the concentration profile. The decrease of chemical potential 


drop with increasing Young’s modulus goes together with an increase 
of electric potential drop. The mechanical contribution pa is found 
to be small. Overall, the electrochemical potential drop Au is barely 
influenced by a change of Young's modulus. 

The mechanical driving forces in Eq. due to swelling show a 
pronounced influence on ion transport. Especially for stiffer systems, 
the influence of stress driven diffusion is on the order of the diffusion 
mechanism and must be accounted for. Although we see a strong impact 
on the concentration distribution in the simulation, the absolute value 
of mechanical stress is small and therefore not in the regime of inelastic 
deformation. For soft materials, osmotic pressure leads to small strains 
and stresses which induce (small) transport contributions in the direction 
opposing the concentration increase. On the other hand, for stiffer mate- 
rials, the importance of osmotic pressure is negligible for the parameter 
set investigated. Transport effects due to volume change are small in all 


cases. 


4.5 Mechanically Deformed Electrolyte 


For the last example in this chapter, we study a two dimensional problem 
where the solid electrolyte undergoes an inhomogeneous external defor- 
mation. This can be caused for example by inhomogeneous deposition on 
a metal interface on one side of the electrolyte layer, i.e. the growth of 


an intrusion penetrating the separator. The change of effective transport 
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paths is an obvious outcome, but also a stress field in the host will 
build up which couples with ion transport. With this example, we want 
to demonstrate the importance of the large deformation formalism for 
correctly predicting ion transport in a deformed solid electrolyte and we 
want to show that external deformation has a strong influence on the 


concentration distribution within a solid electrolyte. 


4.5.1 Setup 


We use the the full model as given in Eq. (3.104) to solve this problem. 
The formulation includes various aspects of coupling between mechanical 
deformation and electrochemistry. An analysis of the different contribu- 
tions in the fully coupled system is difficult since all effects interplay with 
each other. Therefore, in addition to the full model, we consider two 
simplifications to isolate specific effects and to investigate the importance 


of the individual coupling mechanisms: 


Setting A BS =0 (No swelling) O(c) (No spatial diffusion) 
Setting B 8°=0 (No swelling) nO (ce) (Spatial diffusion) 

Setting C ß°=1 (With swelling)  pC(c) (Spatial diffusion) 

Both Settings A and B exclude the swelling mechanism and therefore 
consider no contribution of stress driven diffusion. In addition to this 
simplification, Setting A excludes the contribution of the change of volume 
to the chemical energy with uS (c°) such that e.g. no osmotic pressure is 
present. Therefore, Eq. (3.116) simplifies with ¿£S = 1 and €? = €% = 0 
and only the change of transport path due to deformation is taken 
into account. On the other end of the spectrum, Setting C includes all 
features. For this calculation we use the electrochemical data measured 
for a PS-PEO block-copolymer with LiTFSi [57], see Section [4.2.1] and 
E = 300 MPa. 
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Figure 4.8: Boundary conditions for Example Mechanical deformation with 
displacement BC is applied at Phase 3. P1 and P2 mark the probe points for Fig. 


4.5.2 Boundary Value Problem 


We consider a rectangular domain of size 50 ym x 30pm with periodic 
BC left and right, see Fig. and apply plain strain conditions. The 
top edge is fixed in space and grounded (® = 0). The bottom edge will 
be deformed into a Gaussian-like shape with ugc = [X, APw(X/RP)], 


where we utilize the Wendland function 
wer) =(1—[Irll)*0+4 Ilr) Wr <1 (4.4) 


with amplitude AP and radius RP to define the protrusion with compact 
support of size 2RP. The amplitude is chosen to be AP = 5 pm and the 
radius is RP = 20m. Further, we impose a current density J®PP on the 
top and bottom edges. To highlight the results of our theory, we will 


successively apply varying boundary conditions in six phases: 
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Figure 4.9: Spatial concentration c and electrochemical potential w+ measured at the 
tip (P1) and the flank (P2). 
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Phase 1 serves as a comparison datum for a current flow on the deformed 


geometry in Phase 5. The relaxation Phase 2 is used to obtain a steady 


state subsequent to a transient phase. In Phase 3, we apply a deformation 


much faster than the diffusion time to highlight the effect of mechanics. 


Since we obtain via the next relaxation in Phase 4 a steady state once 


more, the features in Phase 5 and Phase 6 are independent of the rate 


of deformation as long as a purely elastic response is considered. The 


outcome would coincide with a slowly penetrating protrusion. 
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Figure 4.10: Color contour plots showing the salt concentration and the electrochemical 
potential in the electrolyte immediately after deformation (end of Phase 3). The 
spatial concentration shows a peak close to the protrusion tip. The electrochemical 


potential shows no gradient for Setting A, small effect for Setting B and a strong one 
for Setting C. 
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Figure 4.11: Color contour plots showing the salt concentration and mechanical 
pressure in the deformed electrolyte for the relaxed state in equilibrium (end of 
Phase 4). The concentration in Setting A does not relax, Setting B shows complete 


relaxation and in Setting C the ions are pushed away from high pressure regions. 


a | SettingA b | Setting B C | Setting C 


x10°m -2 - x10°m 


Figure 4.12: Color contour plots showing the salt concentration and the electrochemical 
potential in the electrolyte during galvanostatic charging (Phase 5). Setting A still 
shows a peak at the concentration tip, an almost linear concentration polarization 


occurs in Setting B and Setting C shows, even during concentration polarization, a 
depletion of salt at the tip. 
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4.5.3 Results 


Fig. [4.9] shows the evolution of electrochemical potential w and spatial 
concentration c at two probes attached to the boundary of the material. 
One probe is associated with the tip of the protrusion (Xp, = [0,0]) and 
the other is on the flank (X py = [-15 pm, 0]), see Fig. |4.8 

During application of current in Phase 1, a concentration gradient 
builds up. Therefore, the concentration at probes Pl and P2 increases in 
all settings, see Fig. Since there is no deformation, the results for tip 
probe P1 and the flank probe P2 coincide. The concentration gradient 
is not as pronounced if stress-driven diffusion is taken into account and, 
therefore, Setting C has the least increase of salt concentration. The 
electrochemical potential is linked to the salt concentration and hence has 
the same characteristics. Afterwards, in Phase 2, the salt distribution 
relaxes to the reference concentration c' for all settings. 

In Phase 3, the volume above the tip (P1) is compressed with Je = 0.8, 
where swelling effects are small (J, ~ 1) due to the fast deformation 
process. Therefore, the same amount of ions is located in a smaller 
(compressed) volume. Fig. [4.10] shows contour plots of the salt con- 
centration and electrochemical potential immediately after deformation. 
If no swelling is considered, the salt concentration increases to a value 
of c & wee = 1.25c"*. Indeed, this happens in Setting A (Fig. (4.9), 
but the electrochemical potential stays constant because of its material 
description (Fig. |4.10p). For the other two settings, this value is not 
reached. The diffusion mechanism due to increased chemical potential 
counteracts the concentration build-up from the beginning of the defor- 
mation. The peak of salt concentration at Pl in Fig. |4.9p increases 
only up to c ~ 1.08f. The mechanical pressure locally reaches 50 MPa 
and therefore, in Setting C, stress driven diffusion plays a major role, 
increases w+ locally and therefore grad (w+) dramatically, and squeezes 


ions away from the tip immediately. The snapshot, Fig. |4.10¢, taken 
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directly after deformation, shows depletion in the region of high pressure. 
An opposite effect can be seen in P2. This region experiences a small 
tension and thus an increase of volume (Je ~ 1.02). Therefore, the ions 
occupy more space and the concentration decreases. 

During relaxation after deformation (Phase 4, Fig. |4.11), no change 
of salt concentration occurs in Setting A since the chemical potential 
defined with respect to the material concentration is "not aware" of 
the deformation. Its value stays constant as seen in Fig. and 
consequently the gradient, which drives transport, is zero. Setting B, 
in contrast, equilibrates the concentration profile to a value slightly 
above cf. This is a consequence of the no-flux-boundary conditions and 
conservation of mass in the deformed computational domain. Since ® 
is prescribed on the top boundary, the change of 50 me also affects the 
electrochemical potential which increases to 5mV. The time constant to 
reach steady state is comparable to Phase 2. The response in Setting C 
shows again a gradient in salt concentration throughout the electrolyte. 
The salt concentration decreases at the tip with c = 0.87c"* (Fig. |4.11Ẹ) 
and therefore goes even below the reference concentration. The high 
pressure in this region makes the tip an energetically less favorable 
place for the ions. On the other hand, the region of negative pressure 
experiences a small increase of concentration. 

In Phase 5, we solve for galvanostatic charging on the deformed ge- 
ometry (Fig. [4.12). All of the three settings account for the deformed 
transport path associated with Eq. (3.84), i.e. the electrochemical poten- 
tial at the tip (P1) needed for transport is 6 to 9mV less than to that at 
P2 (Fig. 4.9). A concentration polarization with less salt accumulation on 
the tip (Pl) and more ion accumulation on the flank (P2) originates and, 
roughly speaking, is superposed to the steady state of Fig. There- 
fore, the very high concentration in Setting A is observed. Setting B -on 
the other hand- shows a moderate increase of salt concentration at P1, 


an effect of the shorter transport path (Fig. |4.12b). The gradient of 
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concentration in the electrolyte is almost homogeneous and therefore, 
one finds the highest concentration in the region furthest away from 
the counter electrode, i.e. in P2. Setting C shows a similar increase of 
concentration at both probes, Pl and P2, but the squeezing effect is still 
dominant and the concentration difference between Pl and P2 remains 
at 300 me. Nevertheless, due to the additional contribution of pressure 
in the electrochemical potential, the difference of w} between P1 and P2 


is smaller compared to the other two settings. 


4.5.4 Discussion 


These examples lay out the complex coupling mechanisms for ion transport 
in an inhomogeneously deformed solid electrolyte and illustrate the various 
contributions to the coupling in our model. Each setting (the simplified 
Settings A and B as well as the full Setting C) with their different 
assumptions shows not only qualitative, but also quantitative differences 
in terms of the mechanical and electrochemical response to a given 
penetration. The missing relaxation mechanism after deformation in 
Setting A is not compatible with the understanding of mobile ions within 
an electrolyte. Similarly, Setting B neglects the property of ion volume 
and thus omit the energetic unfavorability of salt in high pressure regions. 
Only Setting C provides a consistent coupling including considerations 
of the change of transport path, appropriate relaxation and stress driven 
diffusion. Therefore, we can predict a concentration and electrochemical 


potential profile along a deformed separator. 


4.6 Summary 


Three examples illustrated the impact of mechanics on ion transport 


by using both a transport model distinguishing cations and anions with 
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full considerations of electrostatics and one that only describes the salt 
concentration by assuming electroneutrality. Only the first is capable 
of resolving the double layer at a blocking electrode. In the case of 
galvanostatic charging, the electroneutrality condition holds throughout 
the domain and both models align. Maxwell stresses therefore are neg- 
ligible. Also the osmotic pressure was found to be small for relatively 
stiff materials. Using electrochemical transport parameters from PEO 
based electrolytes and extending them by mechanical stiffness gave us 
the possibility of investigating the impact of Young’s modulus on double 
layers (in blocking electrodes) as well as the salt concentration profile and 
electrochemical potential (galvanostatic charging). The salt concentration 
profile is strongly affected by the stiffness of the host material although 
only small pressure gradients occur. The electrochemical potential, how- 
ever, alters only little. Finally we have studied the response of a deformed 
electrolyte. Only the comprehensive mechanical considerations of the 
fully coupled model enabled us to obtain reasonable simulation results. 
External deformation thereby influences the electrochemical transport 
significantly. The concentration profile along the interface varies strongly 
and suggests that a depletion of salt occurs in the region of highest 
intrusion. 

With a constant boundary current we have excluded interface effects 
in these examples to analyze solely the bulk transport. Reaction kinetics 
at the electrode-electrolyte interface and their coupling with the bulk 
transport need to be described thoroughly for a comprehensive model of 


charge and discharge and will be discussed in the upcoming chapter. 
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Interface kinetics are a subject of non-equilibrium thermodynamics and 
describe the speed of a reduction / oxidation process as given e.g. in the 
context of reaction (1.2). Its generalized form is the monovalent reactions, 


i.e. 
M = M +e, (5.1) 


where an atom M reacts to a cation M* and an electron e`. The reduction- 
oxidation processes are influenced by electric fields and, potentially, 
configurational entropy and mechanics. In analogy to the transport theory 
of Chapter [3] we approach the interface kinetics on a macroscopic level 
since the micro-structure of the interface is either barely understood, as 
in the case of the composition of the so-called solid-electrolyte-interphase 
(SEI) layer, or too difficult to treat computationally, e.g. requiring the 
atomistic structure to be resolved. 

This chapter follows Ganser et al. [128] and discusses in the first 
section the state of the art of modeling reaction kinetics. We then use 
transition state theory from a purely energetic point of view and derive 
the generic form of the Butler-Volmer equation. This is then specialized 
to various combinations of different types of electrode and electrolyte, 
putting the presented theory in the context with commonly used formulas. 
A special emphasis will thereby be placed on the impact of mechanical 
stress and mechanics as they play varying roles depending on the specific 


combination of electrode and electrolyte considered. The mechanisms 
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behind deposition and stripping of a metal electrode will be discussed in 
detail. In the course of our formulation of the electro-chemo-mechanical 
reaction kinetics at the interface we suggest a new form of the Butler- 


Volmer equation for describing the influence of mechanics. 


5.1 State of the Art 


The Nernst equation and the Butler-Volmer equation are two widely 
used formulae for describing interface electrochemistry. The Nernst 
equation —an equilibrium description— correlates the difference between 
the energetic states of the adjacent materials, usually described by their 
activities, and the difference in their electric potentials [143]. The Butler- 
Volmer equation [70] on the other hand constitutes a non-equilibrium 
description that links the current flux at the interface to the electric 
potential drop across it via Arrhenius kinetics. If chosen consistently, the 
Butler-Volmer equation should thus yield the same result in conditions 
of equilibrium as stated by the Nernst equation. The Butler-Volmer 
equation is therefore the more general formulation. 

Different approaches have been used to derive Butler-Volmer reaction 
kinetics. A prominent one is transition state theory, which is also our 
choice in this work. Our derivation is, however, slightly different to the 
derivations documented in literature so far. The basic idea of e.g. Bockris 
et al. [143], Hamann and Vielstich [144] and Newman and Thomas- 
Alyea is as follows: The current density across the interface between 
an electrode and an electrolyte is governed by three quantities, namely 
concentration-independent reaction rate constants (k*'*, kY) determining 
the frequency of activation, the concentrations of the species present at 
the interface (cm, c+) specifying the amount of species which can react, 
and a contribution due to the activation energy, adjusted by the drop 


in electric potential across both sides ot the interface 6°" — 6°’, where 
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ptr PIY denote the electric potentials at the surface of the electrode 
and the surface of the electrolyte, respectively. Following this approach, 


the net current density then reads 


JBN kW cM exp ( 


Alu + (1 - B)P (8% — 0%) 


RO 
(5.2) 


vere _ BF (Be = oly) 
RO f 


— k% c} exp | 


where jBV is the current density across the interface, defined to be 
positive for ions moving from the electrode to the electrolyte, Ad 
is the activation energy barrier for the forward reaction, eer is the 
activation barrier for the reverse reaction (both defined at reference 
state), F is the Faraday constant, 8 is the symmetry factor, R is the 
universal gas constant and O is the absolute temperature. Several authors 
[71] 145] [146] use Eq. as a starting point for deriving an extended 
Butler-Volmer equation (e.g. to incorporate mechanics) - but each with 
a slightly different model and outcome. In addition, most of the derived 
theories are formulated in a way that makes their application to relevant 
problems very cumbersome. For example, Bockris et al. [143] refer 
to surface concentrations (i.e. moles per unit area) in the prefactors 
of Eq. (5.2), whereas usually volume concentrations are utilized in its 
application (i.e. moles per unit volume). Such volume concentrations are 
commonly utilized in applications of the Butler-Volmer equations because 
they are directly accessible from the theory of volumetric bodies[!] 

Furthermore, only the dilute version of the Nernst equation can be derived 


from Eq. (5.2), as activity coefficients do not appear in its formulation. 


1 To avoid issues with units, some authors simply define the rate constant as a reaction 
velocity, see e.g. Christensen and Newman or Hamann and Vielstich [144]. Other 
authors set bulk and surface concentration equal or leave inconsistencies in 
their derivation [119]. 
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Bazant and coworkers resolve this issue by extending the reference Gibbs 
energy by a non-ideality contribution, which is denoted as an excess 
chemical potential [112] [148] [149]. The comprehensive work of Bazant 
[112] describes various interfaces under non-ideal conditions with a special 
emphasis on phase transition materials such as LiFeP O, in the framework 
of non-equilibrium thermodynamics. 

Other aspects of the formulation of electrochemical kinetics at interfaces 
and their relationships to the Butler-Volmer equation have been addressed 
and shall be outlined in the following. As an alternative approach to 
transition state theory, one can use the law of mass action to 
derive the Butler-Volmer equation, as seen in Latz and Zausch [151]. 
A rigorous asymptotic analysis was deployed by Dreyer et al. [152] to 
motivate the jump of electric potential in the interface region. Yang 
[147] extends the Nernst equation to the non-equilibrium condition and 
considers (surface) stresses in the overpotential via a Taylor expansion. 
Lai and Ciucci [122| propose an exponential relation between the current 
and the electrochemical potential drop across the interface. However, this 
drop is not modeled as an energy landscape with an activation energy to 
be overcome. We note that all the derivations mentioned above lead to 
the same overall form of the Butler-Volmer equation. However, they yield 
slightly different prefactors, especially when mechanical contributions are 
of interest. 

A somewhat more sophisticated formulation is the electron transfer 
theory by Marcus and Hush, which assumes that the charge transfer 
kinetics are dominated by a reorganization of the solutes in the adjacent 
layers due to a change of charge. Well defined experiments on surface- 
bound redox couples showed that the electron transfer theory is superior 
to the Butler-Volmer equation, though both theories align for small 
current densities [153]. Further extensions are possible including effects 
such as electron tunneling in metal electrodes [153] [154], N-shaped free 


energy landscapes in intercalation materials or proton-electron 
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interactions to increase accuracy. A quantum mechanical approach 
for the description of reaction kinetics provided by the Marcus theory is 
also used e.g. in [112] 1574159]. Finally, we mention that the diffusive 
layer at the interface can be explicitly resolved through a Poisson-Nernst- 
Planck formalism. In this way only the remaining Stern layer needs to be 
modeled by a Butler-Volmer equation [65] [160]. However, we will restrict 
our analysis to the Butler-Volmer equation as it is commonly used and 
parameterized in the field of energy storage, especially in contexts where 


influences of mechanics are important such as solid state batteries. 


5.2 General Framework 


Following Bockris et al. [143], we utilize transition state theory to describe 
the reaction kinetics on an electrode-electrolyte interface. We restrict 
our considerations to a monovalent reaction, see Eq. (5.1). Intercalation, 
de-intercalation, deposition and stripping of lithium serve as concrete 
examples. As seen in Fig. the material state changes from the 
electrode phase on the left hand side (in this case a metal electrode) 
to the electrolyte phase at the right hand side. This process follows a 
specific energy landscape along reaction coordinates. Note that these 
coordinates are a theoretical construct and are not necessarily related to 
spatial dimensions. Proceeding along the reaction coordinate may thereby 
describe various processes such as adsorption, electron transfer, transport 
through an SEI layer, passivation etc. The energetic states of the reactants 
in their reduced form G®* (in the electrode) and in their oxidized form 
(in the electrolyte) G*Y are specified by Gibbs energies which consist 
of enthalpic, entropic and mechanical contributions. At this point, we 
assume the absence of an electric potential and postpone consideration 


of electrostatics. In the following, we use the notion of partial molar 
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Figure 5.1: (a): Illustration of a metal electrode-electrolyte interface with a mechanical 
contribution due to electrolyte pressure and normal stress at the interface. The energy 
barrier G? is associated with the most likely deposition/stripping reaction $, involving 
the first atomic layer of the electrode. Reactions into hidden layers, e.g. reaction in 
the second layer f are not considered. (b): Energy landscape for a reaction involving 
the first layer at the interface between an electrode and electrolyte without an electric 
contribution (black) and with an electric contribution (blue). The double headed 


arrows indicate that the energetic level of electrode and electrolyte can change e.g. by 
ref 


deposition YIP, swelling Q4p°!Y or configurational entropy RO In (c+ /cte 


Gibbs energy instead of chemical potentials as it is commonly used in 


the context of reaction kinetics. 


5.2.1 Gibbs Energies of Electrode and Electrolyte 


The most prominent representatives of electrolytes are binary ion con- 
ducting systems as discussed in Chapter [3] We assume an additive split 


of the Gibbs energy of the atoms in the electrolyte, in energy per mole, 


Gy = Y" Gr,  I*%% = (ref, conf, ni, mech, V,...}, (5.3) 
kerey 
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where the different contributions 7°% such as reference value, configura- 
tional entropy, non-idealities and mechanics are dependent on the chosen 
electrolyte. Note that electrical effects are excluded. 

Recalling the chemical potential of a solid electrolyte of Eq. (3.77), and 
noting that for pure substances the partial molar Gibbs energy equals 


the chemical potential, we define the different contributions 


Cors = ROln(c+ /c), (5.4) 
Gy = ROMAN (5.5) 
a= a, (5.6) 


ely 


in addition to the reference value of the Gibbs energy per mole Gref 


defined at standard pressure and at a reference value of concentration of 


the cation. If the concentration cy deviates from c, one has to take the 


conf 
holds for dilute solutions. Eq. (5.5) describes non-idealities, as usually 
present in concentrated solutions, see Section The mechanical 


contribution of Eq. (5.6) in a binary solid electrolyte is a consequence 


configurational entropy G“ „into account. The contribution of Eq. 


of swelling and shrinkage of the electrolyte due to local change of salt 
concentration. Note that we describe the reaction kinetics in the spatial 
configuration and for simplicity we assume volume change in Eq. (3.77) 
with JeA® = 1 in this chapter. 

Electrodes in battery applications are often aggregates of active parti- 
cles where the ions are intercalated and de-intercalated in a host material, 
or metal electrodes where the ions are deposited and stripped layer by 
layer. We will consider four representative situations — a metal electrode 
with mechanical effects neglected, a metal electrode with mechanical 


effects taken into account, an active intercalation particle in the limit of 
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dilute concentrations of stored cations without mechanics and an inter- 
calation particle with a high concentration of stored cations and with 
mechanics taken into consideration. 

Similar to the electrolyte, we define an additive split of the Gibbs energy 
of the atoms in the electrode, without considering electrical effects, in 


energy per mole, 


Get = 5 Ge. I°" = {ref, conf, mech, ...}, (5.7) 
keletr 


where the different contributions in /°" such as a reference value, config- 
urational entropy and mechanics are dependent on the kind of electrode. 

The contribution G¢ is a reference value of the Gibbs energy per 
mole. If we assume that the electrode is at standard temperature and 


pressure, then G“%r is the standard Gibbs energy per mole. The reference 
etr 
ref 


electrode (M) and G<t:]AP for an active particle (AP). It is also related to 


a reference concentration c!, in moles per volume, of intercalated atoms. 
M > > 


values are specific for a system, i.e. we denote G“tE|M for a lithium metal 


The configurational entropy is important if the atom concentration cM 


deviates from ci: 


Get. = ROln(cu/a). (5.8) 
For higher atom concentration cm in the active storage particle, one has 
to take the maximum capacity of the storage particle in moles per unit 
volume into account. Each spot is either occupied by an intercalated atom 
or a vacancy. The latter has a concentration cy and the concentration of 
total sites is given by cp" = cm + cv. The reaction is then 


M = Mt +e +V. (5.9) 


108 


5.2 General Framework 


Since the vacancies appear in the reaction as oxidation products, their 


Gibbs energy 
GSP = RO In (eR — em) /AF), (5.10) 


is added to the electrolyte energy G% of Eq. (5.3). A mechanical 
contribution can be found in both active storage particle and metal 
electrode. The former is 


Get i = up, (5.11) 


mech 


where Qm is the partial molar volume of the intercalated atoms in the 


r 


storage particle, and p°" is the pressure in the storage particle minus 


standard pressure [63]. For a metal electrode, we have 
Guar = ee, (5.12) 


where the contribution W@P to the Gibbs energy due to deposition and 
stripping will be discussed in Section [5.3] 


5.2.2 Reaction Kinetics 


Now consider the reaction and the path it takes, still in a setting where 
the electric potential is zero everywhere so that electrostatics need not be 
considered. The reaction itself is a complex and nonlinear phenomenon 
with different multidimensional reaction paths possible. We choose the 
most likely path and associate a reaction coordinate between O and 1 
to the process, as shown in Fig. |5.1| The reaction coordinate is not 
necessarily associated with a spatial length but, obviously, 0 corresponds 
to conditions in the electrode and an energy per mole for the cation equal 
to G° and 1 to conditions in the electrolyte and an energy per mole for 


the cation equal to G°. Somewhere between 0 and 1 on the reaction 
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coordinate axis there is an energetic barrier G* which is the highest point 
in the energy landscape as seen in Fig [5.1P?] The higher the barrier, the 
less likely is the successful completion of the reaction for each attempt. 

Within transition state theory it is difficult to link a clear material state 
to the reaction coordinate and consequently it is not trivial to model the 
energetic state of the barrier. In analogy to the Bronsted-Evans-Polanyi 
principle (used e.g. in Huang et al. [161] and Bazant [112]) one can 
assume that the barrier has a relationship with the states of the electrode 
and electrolyte such that the barrier energy can be fully dependent, 
partially dependent, or completely independent of the energies of the 
cation in the electrode and the electrolyte. We therefore introduce the 


barrier energy, in energy per mole of the cation, as 


GaGie+ SO vay + SY age (5.13) 


ke [ely keletr 


where Gr is independent of the reduced (electrode) and oxidized (elec- 
trolyte) states and the prefactors 6%" and ad account for the relationship 
between the barrier energy and that of the cation in the electrode and the 
electrolyte. Eq. (5.13) can be regarded as an extension of Bazant [112], 
who applies the Bronsted-Evans-Polanyi principle solely to the reference 
energies. The prefactors have, as yet, an undetermined character, but 
possible choices for them will be discussed and utilized below. In general, 
the barrier can be expressed in terms of a transition state activity coeffi- 
cient and it is explicitly dependent on the chemical or mechanical state of 
the activated complex, see e.g. Bazant [112]. However, these properties 
of the activated complex are not easily accessible. In the absence of 
electric potential, Gł — G°* is the energy barrier height standing in the 


way of a successful anodic reaction (oxidation). Similarly, in the absence 


2 The energy barrier corresponds to the state of an activated complex. 
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of electric potential, Gł — G%Y is the energy barrier height standing in 
the way of a successful cathodic reaction (reduction). 

We now consider the energy landscape in the presence of an electric 
potential. It changes if we consider charged species and electric potentials 
PtT on the electrode and P% on the electrolyteß] The Gibbs energies 
for the electrode and electrolyte for monovalent reactions, including the 


effects of electric potentials, are then, in energy per mole of cations, 
. 1 
GE =G* + Fe, GS" = GY + Fe, (5.14) 


where the second term is equivalent to the electrochemical potential of 
Eq. (3.115). Since the bulk properties of the interface regions (e.g. an 
SEI layer) are unknown, classical bulk theories for electrostatics like the 
Gauss law cannot be used to resolve the exact spatial distribution of the 
electric potential. Nevertheless, we know that the electric potential has 
to change continuously from the electrode to the electrolyte. Therefore, 
we make the assumption that the electric potential varies linearly along 
the reaction coordinates. Note that this assumption is tantamount to 
use of the value of the potential as a reaction coordinate. The energetic 


state of the barrier, due to electric potential, then becomes 
Gh = GË + F (8% + 8 (0 — 99), (5.15) 


where 8, the symmetry factor, expresses the position of the peak value 
of energy in the energy landscape, i.e. the position of GĦ on the reaction 
coordinate, see Fig. Marcus-Hush theory indicates that 6 = 0.5 
and aligns with several experimental findings on symmetric reactions 


[153]. However, more recent studies have revealed cases of asymmetric 


3 Most authors we found use only the electric potential drop instead of the absolute 
values. We will see that only the drop is relevant and the absolute values do not 
affect the reaction rates. 
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heterogeneous electron transfer processes which do not allow an accurate 
parameterization in terms of Marcus-Hush theory unless one assumes 
unequal force constants for the solvent reorganization in reduced and 
oxidized state [162] [163]. We will proceed with the symmetry factor as a 
parameter. 

If an ion is "sitting" on the electrode and attempts to jump into the 
electrolyte (anodic reaction) or is "sitting" in the electrolyte and attempts 
to jump onto the electrode (cathodic reaction), it has to overcome the 


activation energies 
AGs = Gt, — GF, AGs = G4 - GY (5.16) 


for the forward (anodic) and backward (cathodic) reactions, respectively[] 


Use of Eq. (5.13), and yields 
AGs =i + > (ae — DG + Y ra 


keletr keley 
JE (8 PER 1)F (pe _ po) , 
AGs = Gi, + > a + Y (rc 


ref 
keletr kerey 
L BF (om BR po) . 


(5.17) 


Having established the activation energy barriers, we can now compute 


the rates of the reactions taking place at the interface. For this purpose, 


4 Bockris et al. use the opposite definition of the arrows. The right arrow in 
their publication describes electronation which is a cathodic reaction. 
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we use the Eyring theory [113] where the rates are determined from the 


activation energy through 


(5.18) 


where 77 it is the success frequency for the forward (anodic) reaction, kp 
is Boltzmann’s constant, h is Planck’s constant and is the success fre- 
quency for the reverse (cathodic) reaction] Note that the configurational 
entropy contribution in ACs and AG‘ yields a concentration dependency 
of F and F. Excluding the configurational entropy, one can identify the 


remainder of Eq. (5.18) with classical rate constants. The Boltzmann 


AGs 


RO 
universal gas constant R gives the probability of the ion successfully 


factor, for the forward reaction exp (- ) with temperature O and 


accomplishing the jump, from the electrode to the electrolyte, with an 
equivalent Boltzmann factor for the reverse reaction. The coefficient EL 
is the attempt frequency for such jumps, see. e.g. Bockris et al. [143]. 
Full consistency of Eq. (5.18) with Eyring theory requires the exclusion 
of the dependency of Gs on the configurational entropy of the reactants 
by setting 
cont = Scone = OV” = 0. (5.19) 
This choice is based on the understanding that Gibbs energies of reactants 
naturally depend on their concentrations, while a transition state only 
exists at a single point in time and, in this instance, is not an ensemble 


of states. We emphasize that values for Be a and 6°, other than 


5 The transmission coefficient occurring in the Eyring theory is assumed to be one. 
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zero would lead to invalid kinetic rate laws even for simple chemical 
reactions independent of electrical potential. See also Savéant [156] for a 
critical discussion of previous studies that used invalid assumptions in 
the context of electron transfer processes. 

To obtain an interface current from the success frequencies in Eq. (5.18), 
we introduce an area density c*, in moles per unit area, of active surface 
sites. This property is independent of the actual amount of cations in the 
electrolyte and atoms in the electrode. In the case of a metal electrode 
this might be associated to the amount of atoms per area at the interface. 
For a sufficiently smooth surface the density c® is directly related to the 
lattice constant and therefore a material constant. Implicit in this choice 
is the concept that only the outer layer (see 1, in Fig. [5.1) is able to 
react with rate 7; the reaction in hidden layers (e.g. t5) are unlikely 
to occur because their energy barrier is much higher than GÈ, and so 
their reaction rates will be negligible. Therefore, for a metal electrode, 
as shown in Fig. the anodic reaction (stripping) then gives rise to a 
current density 

j =FS?. (5.20) 
We now argue that the same area density, c5, controls the reverse (cathodic 
deposition) reaction. Consider an attempted jump of an ion from the 
electrolyte to the electrode. If the target site of that jump is already 
occupied by an atom (e.g. of an inner layer), there exists an additional 
energy barrier for completion of the jump in the form of repulsion from 
that site. Thus the success rate for such a jump will be negligible. 
Successful jumps are therefore only possible to sites on the surface whose 
area density on the metal electrode is given by c?. Therefore, the current 


density associated to the reverse (cathodic deposition) reaction is 


GS = FOF. (5.21) 
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A similar argument holds for non-metallic active storage particles but 
with the area density cS associated with the intercalation sites at the 
interface. Note that the entropic energy contribution and not cS causes the 
forward or backward current density to vanish if the volume concentration 
approaches zero. For example, Eq. (5.8) and Eq. (5.4), respectively, yield 


0 ee > see F0 EN 
c4>0 => GY +-—o => AGs > o > 70 


In the case of a metal electrode, the anodic current density cannot vanish 
since a reaction partner is always available on the surface. Surface 
passivation, however, can decrease the amount of active surface sites c° 
and limit the overall current density. But such aspects of the behavior 


are outside the scope of this work. 


Combining Eq. (5.20) and with Eq. (5.16) and yields the 


anodic and cathodic current densities 


7 = doo exp ( Cros + Exam (R = UGE + Vier eo) 


RO 


(1 _ B)F (oer _ pel) 
exp RO , 
Eo GH + Dreri + ers pa 7 DES 
J = Joo EXP 
RO 
BF (De u ory) 
exp RO 5 


(5.23) 
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where the prefactor is given by 


kp® 
joo = y PE (5.24) 


We note that Eq. and Eq. differ in that independent reaction 
rates and cation concentrations appear as coefficients of the anodic and 
cathodic current densities in Eq. (5.2), whereas they are identical in Eq. 
(5.23). The concept that the forward and reverse current densities are not 
explicitly a function of electrode and electrolyte cation concentrations is 
also applied by Gutman [150] who derives the current densities via the 
kinetic law of mass action. 

Before evaluating the effective current densities, it is beneficial to 
consider a steady state where forward and reverse current densities are 
equal, and thus the net current density is zero. Under these conditions, 
the electric potential difference equals the equilibrium or open circuit 
potential V¿* = (@°* — pely) “1. The individual anodic and cathodic 
current densities prevailing at this equilibrium potential difference are 


then called the exchange current densities and read 


Jo =Joo exp 


— eq 
E DIVE 


: 3 1 el 
=Joo exp ( Clos BA kert One GR + peter (& 7 na) 


Chos + J kert (ojo 7 1)Gy* + Dike rely ee) 
RO 


RO 


(5.25) 
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We deduce from the equality of the two lines in Eq. (5.25) a general 


version of the Nernst equation between the electrode and electrolyte as 


FV =G- G", (5.26) 


in agreement with Bockris et al. [143] who derive the Nernst equation 
from consideration of thermodynamic equilibrium. 

We emphasize that the result connecting Eq. and is 
a consequence of the equality of the prefactor joo in Eq. for 
forward and reverse reactions. If the anodic and cathodic terms in Eq. 
had coefficients that differ from each other and that are not both 
equal to joo, the prefactors would enter Eq. (5.26) as exponents and 
a conventional Nernst equation would be unattainable. Note that the 
equilibrium potential in Eq. (5.26) is independent of the barrier height, 
the symmetry factor and the coefficients dg and ge f 

Using Eq. (5.26) we rewrite the exchange current density from Eq. 
(5.25) as 


_ gi de En getr get: 
Jo wen ( Sa) exp (== (B= ORIG ) 


RO 


Feri- pn 
RO 


(5.27) 


It is not surprising that the exchange current density depends on the 
energetic states of electrode and electrolyte. The consequence of this 
will be discussed in the upcoming section. Note that the open circuit 
potential (Eq. (5.26)) and the exchange current density (Eq. (5.27)) 
are both determined by the same energetic states of the electrode and 


electrolyte and thus are not thermodynamically independent [112]. 
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electrode electrolyte 


Ei > F < œ applied overpotential 


Ger (VET +n) 


etr eq 
G +FVs 


0 pl reaction coor. 


Figure 5.2: The energy landscape at open circuit conditions (blue) attains an equal 
level at electrode and electrolyte. In the case of an anodic reaction (7 > 0), the 
equilibrium of forward and reverse reactions are violated and a positive net interface 
current arises. 


> ,+ 
Now we consider the non-equilibrium case where j # j. The net current 
density in the direction from the electrode on the left side of the interface 


to the electrolyte on the right is defined as 


PY a= jj. (5.28) 


and is therefore anodic if positive and cathodic if negative. With the help 
of Eq. (5.23) and Eq. (5.27), we obtain the well known Butler-Volmer 


equation 


= exp (a B) 25") exp ( 675) (5.29) 


where we have introduced the surface overpotential 


n= 0% _ pv _ yea (5.30) 
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as a function of the open circuit potential V¿* given in Eq. (5.26). Note 
that only the exchange current density jo in Eq. (5.29) depends on the 
barrier, and the term in the brackets does not. If the overpotential is 
positive, the equilibrium given by the Nernst equation (5.26) is violated, 
the Gibbs energy of the electrode increases and consequently, the an- 
odic current density j dominates the cathodic current density j (see 
Fig. [5.2). The net current density is then positive. If the overpotential 


is smaller than zero, a negative current density prevails. 


5.3 Mechanical Energy of Deposition 


In the previous section we presented a generic concept to link the energetic 
state of the electrode and electrolyte to the interface reactions which are 
usually accessible from bulk theories. Active storage particles [63] and 
solid binary electrolytes [IOI] have shown a strong coupling mechanism 
between mechanical stresses and electrochemical transport due to a 
swelling mechanism. Each ion inserted or removed displaces surrounding 
material, i.e. the host, due to its molar volume. Consideration of these 
concepts for solid electrolytes and for intercalation storage particles yield 
the mechanical effects accounted for in Eq. or Eq. 6.11). We note 
that experimental verification for the influence of mechanical stresses 
on electrochemical response has been obtained, for example, for silicon 
active particles [164] and energy harvesting applications [76]. The work 
of Pannikkat and Raj on the correlation of mechanical stress and 
equilibrium potential of a platinum-zirconium interface also proves that 
mechanics influences reaction kinetics. 

In contrast, deposition on and stripping from metal electrodes is quali- 
tatively different from intercalation, de-intercalation and solid solution 
in electrolytes. The stress free size of the electrode grows/shrinks atom 


layer by atom layer at the oxidation/reduction process and thus cannot 
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Electrolyte 


Ah 


Electrode 


Figure 5.3: A sketch of a solid state half cell loaded by a dead load F. The weight is 
lifted by Ah during deposition. 


be treated as a volumetric swelling process. This change of geometry 
is likely to cause deformation and therefore contributes to the overall 
energetic state of the metal electrode. In this section we consider the 
mechanics of deposition on and stripping from metal electrodes and its 


implications for electrochemistry. 


5.3.1 Macroscopic Motivation 


Imagine a half-cell consisting of an ion conducting solid electrolyte and a 
metal electrode that rest on a rigid support. The electrolyte is connected 
to a second electrode which serves as an ion supply. Furthermore, the 
electrolyte is assumed to have very good transport properties such that a 
homogeneous deposition on the metal electrode takes place when a current 
is applied. A weight is placed on the solid electrolyte and consequently 
a force F acts on the electrolyte. As the material is deposited on the 
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electrode with a current density jPY, the thickness of the metal electrode 


increases by 
Ah = a (5.31) 


for each time increment At, where Qy is the molar volume (defined in 
the undeformed state). As a consequence, the weight is lifted by this 
amount and therefore, the potential energy of the weight rises by 


Q 
AW, =FAh = P JEY At. (5.32) 


Obviously, the potential energy decreases in the case of stripping. To 
capture the full mechanical response, we also have to take into account 
the mechanical deformation within the newly deposited material. In 
the case of a linear elastic response, we can use a spring model. The 
stiffness is thereby k = EME A/h with Young’s modulus of the electrode 
EME, electrode thickness h and area A. The material will be compressed 


following Hooke’s law 
F = ku > u=F/k, (5.33) 


by the displacement u. The new height of the metal is then h=h-uand 
the generated energy due to deposition is composed of potential energy 


and elastic energy 


AW, = FAh+ ku? = F(Ah — 5u) = FAR(1 — u ). (5.34) 


In cases where Er is small (i.e. if forces are small and stiffness high), 
only the mechanical energy due to the displacement as in Eq. (5.32) is 


relevant. 
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Now move attention to the impact on the electrochemistry at the cell 
level. To store the electrochemical energy AE by depositing a certain 
amount of ions, one has to apply the voltage U. The energy supply in the 
half cell is then ¡BY AU At and splits to electrochemical and mechanical 


energy 
AE + AW, = jPY AU At. (5.35) 


Consequently, to store electrochemical energy AE at current density JPY 
requires 
AB=(U-FM jPY AAt. (5.36) 
AF 
Hence, in order to store the same electrochemical energy in the cell, the 


applied voltage must be higher if the electrolyte is loaded by a weight. 


5.3.2 Deposition on lon Scale 


In the previous section we have motivated the fact that the mechanics 
of deposition and stripping, have an effect on electrochemistry. We now 
want to transfer our macroscopic insights to resolve the local energy 
contribution due to mechanics that is required for transition state theory. 
Assume a planar interface, where one layer of atoms is deposited on a 
surface with area A. There are Nm moles of atoms with volume Lu 
necessary to create this layer on the metal electrode and the interface 


then moves by 


Nußm 


Ah = A 


(5.37) 
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We also assume a mechanical force F in the normal direction acting on the 
interface with area A. The energetic contribution due to the deposition 


of Nu moles of atoms is then 


NuÚm 


Wry = AR = 2% 


F. (5.38) 


The molar energy generated by displacing the interface is 


Wy, Ou 
yer = v= F = -Nuo 5.39 
Nu A M n ( ) 
where the normal compressive stress o&" = n? on in the normal direction 


n to the interface and o is the Cauchy stress tensor evaluated at the 
metal electrode interface. Note that the amount of atoms drops out in 
Eq. (5.39) and the concept of discrete atom layers changed to a continuous 
setting. As obvious from Eq. (5.34), we can neglect the elastic strain 
energy contributions for Er <1. 

The result of Eq. (5.39) is comparable with, for example, the findings 
of Monroe and Newman [71] and Ma et al. [145] if making a few strong 


assumptions. For example, we can assume we are dealing with a solid 
that behaves like a fluid with 05% = —p**, We then obtain 


PP = Qyp*, (5.40) 


where the mechanical contribution is identical to that arising from the 


process of swelling during intercalation in an active particle [63]. 


5.4 Examples 


The Butler-Volmer equation derived in Eq. (5.29) is independent of the 


specific (interface) material system. However, the corresponding exchange 
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current density and overpotential need to be adjusted accordingly. We 
will do this for a selection of representative interface systems in the 


following. In all upcoming equations the definition 


$ Gi + (1-8 E ely Yu B- seh CH Z 
le = joo exp ref ( ref ne F E) fl 


(5.41) 


is used, where joo is defined by Eq. (5.24) . The quantity Hee corresponds 
to a reference state and hence is independent of changes in volume 
concentrations of electrode and electrolyte and mechanical state. The 
index Y represents the class of electrolyte and Z the class of electrode. 
We further define 


AG reel = =G 


AG vos] Y 


etr netr |Z 
ref =Ó ef Er ef | 


_ gel Gey (5.42) 


Cos ref Gref 


as the activation energy barriers in the forward and reverse reaction, 
respectively, at reference state and 
Y Goy Y Gett Z 
AG yet|z = = Gl = Gref| (5.43) 


as the difference of standard Gibbs energy between electrolyte and elec- 
trode. 


5.4.1 Dilute Electrolyte | Storage Particle (ideal) 


With above framework at hand, we first investigate the reaction kinetics 
between two ideal systems, namely an ideal active particle (Z=AP), 
using contribution of Eq. (5.8) and a dilute binary ion conducting 
electrolyte solution (Y=BIC), using contribution of Eq. (5.4) where 


in both materials non-idealities and mechanics are neglected. Thus, we 
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have /°" = (ref, conf) and 1%Y = {ref, conf}. Making use of Eq. (5.23) 
together with Eq. (5.28), we write the current density as 


qe =k em exp 


AE + (1 — BF (oe = =) 


RO 
a (5.44) 
= kre, exp ( Gal BR (0% 0%) 
T RO 
with the two rate constants 
a etr (Yetr el el 
petr — J00 exp Feat Ce + eet, GE 
ae RO , 
(5.45) 


ref RO 


3 etr (Yetr ely nely 
kely — I exp ( ere Ch + peroo ôk Gk ) 
+ 


With that, we have determined the exact form of the rate constants in 
Eq. from the introduction. Linear reaction theory implies concen- 
tration independent rate constants which is in line with the finding that 
the barrier is not affected by the configurational entropy of the adjacent 
regions as discussed in connection with Eq. (5.19). Therefore, we use 
Eq. (5.19) in Eq. (5.45) and henceforth throughout this work. 
Eventually, we obtain for the Butler-Volmer equation the open 
circuit potential with respect to the electric potential ® via Eq. (5.26) as 


eq_ 1 BIC cy /ep 
iS T AGzelap + RO In Se E (5.46) 
M 


The exchange current density from Eq. (5.27) computes to 


refyBIC | C+ ee CM il 
Jo = Joo lap ( +) ( =) . (5.47) 
c+ CM 
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The Butler-Volmer equation together with Eq. (5.46) and Eq. (5.47) 
reassembles the form as derived in, for example, Newman and Thomas- 
Alyea [53] and Bockris et al. [143]. Although we arrive at the same 


outcome, we point out that we made the specific assumptions of a dilute 


electrolyte and an ideal active particle. Therefore, if non-ideal systems are 
of interest one has to adjust the formulation of the open circuit potential 


and the exchange current density accordingly. 


5.4.2 Dilute Electrolyte | Metal Electrode 
(w/o mechanics) 


Another idealized electrochemical system is an electrolyte at dilute con- 
centration of cations, see Eq. (5.4), together with a pure metal elec- 
trode (Z=ME), again without considering mechanics. Thus we have 
I°'y = {ref, conf} and 1%" = {ref}. We obtain from Eq. and 
with Eq. (5.19) the open circuit potential with respect to the electric 


potential and exchange current density as 


1 
Ved = = | AG rer + ROIn E JI , (5.48) 
F ee 
ref BIC [ C+ en 
Jo = Joo IME (=) . (5.49) 
+ 


Thus we obtain the commonly used form of the Butler-Volmer equation 
for metal electrodes [60]. 


5.4.3 Concentrated Electrolyte | Storage Particle 


If we consider a high concentration of atoms intercalated in active stor- 
age particles, we have to take the Gibbs energy of the vacancies into 
account, see Eq. (5.10). This energy contributes to the electrolyte Gibbs 
energy, as given in Eq. (5.4) and Eq. (5.5). Together with the Gibbs 


126 


5.4 Examples 


energy of the active particle, using Eq. (5.8) and (5.11) and taking me- 
chanical effects into account we then have 1% = {ref, conf, ni, mech, V}, 
re = (ref, conf, mech} and can derive the open circuit potential with 


respect to the electric potential as 


yea = 1 AG BIC ROI Cr J+ (cM = cm) Q ely Q etr 
6 p refl AP + n ref fref + Slip = MP 
F er Jy CM 


(5.50) 


Accordingly, applying Eq. (5.27), the exchange current density is given 
by 


l= 1—p-6°'" 1- 
j ES) e B E o a 
0 = a M M 
00 AP cre ref = z 


— Setr —_ R _ gely 
exp (£ rech) Omp) exp (2 B nech) ae) ; 


RO 
(5.51) 


where we have already applied Eq. (5.19). The dependence of the ex- 
change current density on the fugacity f+ of the electrolyte and intercala- 
tion concentration in the active storage particles aligns with the currently 
available literature on the subject in regard to modeling if 5°” = 0 is 
assumed The decline of the exchange current density 
for fully charged (cu = cx) and fully discharged (cu = 0) particles as 
predicted by Eq. has been shown to occur for LiFePO, [167]. 
Next, we turn our attention to mechanical effects. The mechanical 
contribution for an active storage particle in Eq. (5.50) aligns with results 
of Stein et al. [168]. The exchange current density depends on both the 


pressure in the active particle p“" and the pressure in the electrolyte p*. 


6 Note that and treat the configurational entropy of the vacancies within 
the activated state and thus differ from the treatment in Eq. 69. The symmetry 
coefficient in the fourth expression of Eq. (5.51) then drops out. 
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These pressures can be significantly different in magnitude. In liquid 
electrolyte systems the pressure in the fluid is usually small. On the 
other hand, active storage particles can experience high mechanical stress 
magnitudes during the intercalation process [62] [63]. Thus the difference 


in pressure between the electrode and the electrolyte must be taken into 
account, see Eq. (5.50). 


5.4.4 Concentrated Electrolyte | Metal Electrode 


The study of pure metal electrode systems which are, for example, present 
in solid state batteries requires, comprehensive consideration of electro- 
chemistry and mechanics. Therefore, for a metal electrode subject to 
mechanical stress we use the energy as defined in Eq. (5.39) and consider 
a binary ion conducting electrolyte (e.g. a polymer) with utilization 


of the energy contributions of Eq. (5.4), (5.5) and (5.6). Hence we 


have I°% = {ref, conf, ni, mech} and /°® = (ref, mech}. The result in 


Eq. (5.26) then gives, together with Eq. (5.39), 


a 1 is 
Vs" = = [AGS + RO In (a) + Qp% + Oma + (5.52) 
+ whe 


Making use of Eq. (5.19), the exchange current density follows 


(1-8) (1-8-87) 
. _;refiBIC [ C+ f+ 
Jo =Joo ME ref ref 


e E A (5.53) 
— setr 1-8- 
exp (Fe oor") exp ( en) ur) 


and is consequently a function of the entropic state of the electrolyte and 


the mechanical states of both the electrolyte and the electrode. 
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To discuss the influence of deviatoric stresses, we split the normal stress 
into a hydrostatic and a deviatoric portion with 


etr etr ely 


On =Ona—P with 079 =n: o!:n (5.54) 


where øf is the deviatoric stress tensor. If the molar volumes of the 
cation in the salt and the metal electrode are such that Q} = My, then 
Eq. (5.52) simplifies to 


1 
Vet = — | AG ¡e HS + RO In oe Opol (5.55) 
F O i 


The open circuit potential V$“ is therefore independent of mechanical 
contributions if o°*, is small. This follows because G3” and G$" increase 
and decrease in the same magnitude if p°” > 0 and p°” < 0, respectively. 
The difference of the energies remains constant. Note that, although 
the mechanical stresses drop out in Eq. (5.55) if a = 0, a gradient in 
mechanical stress induces a change of concentration due to stress driven 
diffusion and thus implicitly influences the open circuit potential. The 


exchange current density reads 


ely 
. ref f C+ (=A) f+ Gob Far} (1 = caen — Ser a Op 
Jo = Joo ref ref exp RO 


cy + 
êz —(8 T IE ca 
P RO ' 


etr 


For small 0,4: 


the mechanical contribution in Eq. (5.56) is independent 


of the symmetry factor 8. Depending on the choice of prefactors, the 


etr ) 
’ 


mechanical contribution in Eq. (5.56) may vanish (9%, +6, = 


; 1 
decrease the exchange current density as the pressure increases (5,77 4, + 


etr ; a sd : 
ôS p > 1) or increase the exchange current density as the pressure in 
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ely 


oY a + en < 1). The latter is a consequence 


the electrolyte increases (ô 


of a decrease in activation energies. 


5.5 Application within Concentrated 
Solution Theory 


Recall that the Butler-Volmer equation obtained in the previous section, 
Eq. (5.29), is derived via a difference in electric potentials. In the context 
of concentrated solution theory [53], the field variable in the electrolyte 


is, however, the electrochemical potential 
wy = FOY + Gey, (5.57) 


see Eq. (3.115), and is therefore not directly compatable with the 
definition of the overpotential of Eq. (5.30). With the help of Eq. 
(5.57) and the Nernst equation (5.26), we rewrite the overpotential from 
Eq. (5.30) as 


w4 — YA (5.58) 


where we have introduced the open circuit potential with respect to w+ 


as 
Vet = ge (5.59) 
wı = F . . 


This overpotential is independent of the energetic state of the electrolyte 
GY, i.e. its cation concentration or mechanical condition. Note that the 
exchange current density (e.g. Eq. (5.27)) is not a function of the electric 
potential, thus remains unchanged and, as before, is a function of G%®. 

Open circuit potentials can be accessed by measurement in an equi- 


librium state, i.e. in the absence of concentration gradients [169]. Both 
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definitions of the open circuit potential, V7 and Vz", differ only by the 
additive constant (eae! in the equilibrium and stress free state. Concen- 
trated solution theory as discussed in Chapter [3] considers implicitly the 
electrolyte energy state in the field variable w} and thus only requires 
knowledge of the electrode in the parametrization of the open circuit 
potential V$. 

At this point, let us briefly recall selected examples of Section [5.4] with 
respect to a formulation with the electrochemical potential as electrolyte 
field variable. Recall that the exchange current densities (Eq. (5.49), 
Eq. (5.51) and Eq. (5.53)) remain unchanged. The metal electrode in 
combination with a dilute electrolyte (Section [5.4.2) yields 


1 
V4 = =a) (5.60) 


Wy ref 


In experiments it is beneficial to select one material as datum and measure 


only potential differences with respect to this material. This is achieved in 


etr | ME 
ref 


for the metal electrode then vanishes. This concept is in agreement with, 


our formulation by setting G = 0. The open circuit potential Ma 
for example, Newman and Thomas-Alyea who use the metal electrode 
as the reference potential. We emphasize that the specific choice of zero 
open circuit potential is different to the notion of standard Gibbs energy 
and further does not take mechanical effects into account. The system of 
Section [5.4.4] with the electrode energy contributions 1%" = (ref, mech} 
(see Eq. and Eq. (5.39)) yields the open circuit potential with 


respect to the electrochemical potential as 
Vict = = (GSM + Oyo) (5.61) 


with the stress dependency arising due to normal stresses. Eq. (5.59) 
gives for they system of Section that describes the combination of 
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concentrated electrolyte and active particles with mechanics (Eq. (5.8), 
Eq. (5.10) and Eq. (5.11)), the open circuit potential 
1 
Ved === Ca + ROln (=e) + ou") : (5.62) 
M 


w+ F er 


The exchange current density is given by Eq. (5.51). Eq. provides 
a clear correlation of VS} and cm based on the electrode model with 
energy contributions of Eq. (5.8), Eq. (5.10) and Eq. (5.11). However in 
battery application, direct measurements of Vo, as a function of cm (state 
of charge) are feasible. Therefore one usually relies on the measured open 
circuit potential instead of using e.g. Eq. (5.62) although the theoretical 
formulation is able to resolve the mechanical contribution. We note 
that Bohn et al. [63] and Christensen and Newman [61] have considered 
deviations from ideal conditions in storage particles by inclusion, inter 


alia, of an excess Gibbs energy for the intercalated cations. 


5.6 Influence of Mechanics 


We have defined a generic form of the energy barrier with respect to the 
energy states of the adjacent electrode and electrolyte in Eq. (5.13). The 
way that the parameters d¢"* and q are chosen has a strong impact on 
the exchange current density. Recall that the overpotential (Eq. (5.30) 
is independent of the barrier and is therefore not affected by go and 
de. We have assumed that the energy barrier is independent of the 


configurational entropy of the adjacent layers. In the following, we focus 


etr 


on the mechanical contribution and discuss four assumptions for ôS% p 


and 6°, as given in Table[5.1] and relate them to literature. 


mec 


Before doing so, we note that Go as introduced in Eq. (5.13), might 


also be stress dependent. In general, it is difficult to determine the 
t 


mechanical contribution to the activated state G/,, due to its unknown 
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een ae implication on the barrier similar work 

0 0 barrier is independent [112] [145] 

1-aM aM linear drop of Gibbs energies [7I] 

B 1—8 equivalent to electric potential [71] [146] [161] [170] 
0 1 only electrolyte affects barrier [71] [72] 


Table 5.1: Different assumptions for the influence of mechanics in the electrode and 
electrolyte energies on the energy barrier. 


location on the reaction coordinate. Therefore we will not consider the 
stress dependency of > in the following. 


A barrier which is independent of the mechanical state of the electrolyte 


ely = getr 


and electrode is given by rach = Omec 


„ = 0. Mechanical stress influences 
the exchange current density in this case, but not due to the Bronsted- 
Evans-Polanyi idea, where the barrier height itself is a function of stress. 
Here, mechanical stress stems from the stress-dependency of the open 
circuit potential (see Eq. (5.62) and (5.61)). Ma et al. [145] follow the 
concept of a fixed barrier, with both 6, and 5°”, equal to zero, but 
come to a slightly different Butler-Volmer formulation than we do since 
they define the open circuit potential to be zero and thus independent of 
stress. It also aligns with Bazant [112] where the Bronsted-Evans-Polanyi 
principle is applied solely to the reference energies of oxidized and reduced 
states. 

A second assumption is to propose a linear drop of the mechanical 
part of the Gibbs energy between the electrode and electrolyte with 
Se, = 1- aM and 8%, = aM where we have introduced a second 


symmetry factor a”, again restricted to the interval between 0 and 1. 
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The exchange current density of Eq. (5.27) as a function of electrode and 


electrolyte energies then reads 


a aM_ = etr _ [ely 
nme SEN (1- BG -G ) có) 


RO RO 


and, for the specific case of a metal electrode and a solid electrolyte (see 
Eq. (5.53)), takes the form 


amla O alae, 
Jo =Joo M (5) ( +) (5.64) 
+ + 
M _ 1- -Q etr _ Q ely 
exp (€ Sa L a +E ) (5.65) 


This structure with symmetry factors is similar to that of Monroe and 
Newman [71], but with only mechanical effects taken into account[”] 
Motivated by the idea that mechanics should solely affect the equi- 
librium potential, Monroe and Newman assume in a second step 
that aM = 1 — 8 (which translates in our formulation to 5° B and 


mech = 
Oe =1- 8). Another motivation for this choice is the idea that the 
difference in the mechanical part of the Gibbs energy along the reaction 
coordinates parallels the difference in electric potential, see Eq. (5.17). 
This approach is utilized e.g. by Lu et al. [146], Huang et al. [161] and 


Joos [170]. 


A second simplification of the derivation used by Monroe and Newman 
[71] is to set a™ = 1, leading to 5°, = 0 and 9%, = 1. This assumption 
is, for example, used in [72] and several subsequent publications studying 


the interface stability of a metal electrodes in solid state batteries. Here, 


+ 
T Monroe and Newman introduces the term Au” as the change of the mechanical 
energy drop at the interface between undeformed and deformed state. This is 


equivalent to GSH h — Ela in our formulation. 
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only the energetic state of the electrolyte influences the energy barrier 
and implies a non-symmetric influence of electrode and electrolyte state 
on the reaction. 

Since all assumptions in Table [5.1] are thermodynamically consistent 
from a theoretical point of view, only experimental validation will enable 
comprehensive understanding of reaction kinetics and clarify which as- 
sumption is most appropriate regarding the energy barrier. Until better 
knowledge is available, we suggest the use of the simplest model with 
oly, = oe, = 0 because no further assumptions, such as the same 
linear drop of energies between electrode and electrolyte, is required. 
Furthermore, with a theory for the thermodynamics of the adjacent com- 
ponents in hand, see Chapter [3] it seems to us more appropriate to treat 


mechanical effects in the same manner as entropic ones. In Eq. (5.19) 


u 1 A 
we have seen that 00% p = 6, = 0 and thus we conclude the mechanical 
1 
parameter to be °? p = Oey, = 


5.7 Interface Properties 


We now want to deduce the interface properties of a polymer electrolyte 
adjacent to a metal electrode. Within the formulation of the extended con- 
centrated solution theory with electrochemical potential as field variable, 
we utilize the equilibrium potential from Eq. (5.61), choose lithium metal 
as reference and thus set G£tt][M = 0 under standard stress conditions. 
Regarding the exchange current density, the interface resistivity Rpy can 
be experimentally measured by, for example, impedance spectroscopy. 
The correlation to the exchange current density is given by linearization 


of the Butler-Volmer equation (5.29) yielding 


F 
-BV . 
ee E 5.66 
j Jone” (5.66) 
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where the symmetry factor $ drops out. With the knowledge of Ohm’s 
law y = RgyjPV we obtain with 


=a (5.67) 


the correlation of exchange current density and the interface resistance. 
Experimental data for an interfaces between lithium metal and a polymer 
electrode, for example, showed an interface resistance in the range of 
6.7Q cm? and 27Qcm? [171]. Under the assumption that the measure- 
ments were done in a stress free environment, Eq. then yields an 
exchange current densities between 1.1 a5 and 4.5 més. These measure- 
ments, however, do not take concentration dependencies into account. 
To resolve that, we recall Eq. (5.53), use the electroneutrality condition 


c = Cy and choose a symmetric reaction with 3 = 0.5. We then obtain 


Q eo Q ely 
jo = KPY c25 exp (- ne ) exp ( aH ) ‘ (5.68) 


0.5 
where we have introduced KBV = gie BIG ( put . Eq. 1.68] is 


similar to but extended by mechanical effects. Neglecting the 


mechanical contribution and assuming an ideal solution with c = 1000 me 


yields for the before mentioned interface resistances values for KPY 


between 0.356 Aa and 1.437 JzA For the upcoming investigation, 


we will use with KBV = 0.589——— an intermediate value which is also 


Vm mol 
used in 


8 Note that the exchange current density in miss the Faraday constant. 
Therefore, one has to multiply the rate constant by the Faraday constant to achieve 
the values presented in our work. 
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5.8 Summary 


This chapter dealt with a consistent derivation of the Butler-Volmer 
equation including mechanics and together with the transport theory of 
Chapter |3} it completes the electro-chemo-mechanical model for a solid 
state electrolyte adjacent to a metal electrode. 

We have used the Gibbs energies of cations in electrode and electrolyte 
bulk and an energy barrier based on the Bronsted-Evans-Polanyi idea 
to define activation energies for reaction kinetics. The reaction rates are 
then related to a surface concentration of cations at the interface, which is 
the same for both the forward and reverse reactions. This yields a generic 
Nernst equation and resolves discrepancies between surface and volume 
concentrations that are present in the literature. It also suggests that 
interface kinetics depend strongly on the thermodynamic states of the 
adjacent bulk materials, i.e. the coupling of open circuit potential, the 
exchange current density and the chemical potential of active particles. 
We have embedded the Butler-Volmer equation in concentrated solution 
theory, where the electrochemical potential is used instead of the electric 
potential as the field variable in the electrolyte. 

We have additionally discussed the energetic contribution of the mechan- 
ics of deposition and stripping at a metal electrode, which is conceptually 
different from that of intercalation and de-intercalation processes. The 
normal stress, not the hydrostatic pressure, is shown to be the variable of 
interest for the thermodynamics of deposition and stripping, with an im- 
pact on both the open circuit potential and the exchange current density. 
We have elaborated the link between our very generic formulation and 
other modified Butler-Volmer equations discussed in literature, with an 
emphasis on mechanical stress. Although not definitive, we suggest inde- 
pendence of the energy barrier from the Gibbs energy in the electrolyte 
and electrode. This is a simple assumption but, with mechanics taken 


into account, constitutes a new form of the Butler-Volmer equation. 
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We now recapitulate the thoroughly derived reaction kinetics for a metal 
electrode adjacent to a solid binary electrolyte which we are going to 
apply in the morphological interface stability analysis in the upcoming 


chapters. The Butler-Volmer equation for this system reads 


F'n F'n 
BV _ 
q 40 op (55) Sap ( ans) | 


n = oOo — > (wi + Quo"), (5.69) 


where we have repeated the Butler-Volmer equation (5.29), the overpoten- 
tial of Eq. (5.58) with the mechanical contribution of a metal electrode 


(Eq. (5.61)) and the exchange current density of Eq.(5.68). 
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A Quasi Static Analysis 


After carefully deriving an electro-chemo-mechanical battery model in 
the previous chapters, we now present our first assessment of morpho- 
logical stability of a metal electrode in a solid state battery. Recall that 
the mechanisms by which lithium dendriting and intrusion into solid 
electrolytes may be suppressed are not fully understood. For example, 
an open question concerning a polymer electrolyte such as PS-PEO is 
whether there is a composition that provides reasonable transport prop- 
erties and, at the same time, achieves sufficient mechanical stiffness and 
strength to suppress intrusions. Other polymer systems struggle with 
the same trade-off. It is widely believed that a stiffer electrolyte is more 
likely to keep the lithium metal under control, but a full understanding is 
still missing. Mechanical properties alone cannot explain morphological 
stability because polymers under consideration have Young’s moduli that 
are an order of magnitude lower than that of lithium. Therefore, it can 
be inferred that coupling between electrochemistry and mechanics in the 
battery plays a key role in any possible suppression mechanism. 

To understand possible suppression mechanisms, we utilize the rigorous 
formulation of coupled electro-chemo-mechanics for bulk ion transport as 
introduced in Chapter [3] and for reaction kinetics as derived in Chapter [5] 
Our approach allows us to utilize a full set of measured electrochemical 
parameters and to solve for the ion transport in the deformed geometry 
caused by lithium protrusion. The model is used to numerically analyze 


the interface current density at a given deformation state. If the current 
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Figure 6.1: The comprehensive simulation model in this study considers electro-chemo- 
mechanical interaction in the electrolyte bulk as well as at the interface in order to 
determine morphological interface stability of a metal electrode. 


density is higher at the tip of the protrusion than at an undisturbed region, 
we expect an unstable behavior and growth of the protrusion. Otherwise, 
we expect stable behavior, see Fig. Note that the basic definition 
of stability used here is somewhat similar to that used by Monroe and 
Newman [72] even though the further considerations and results differ 
considerably. With the stability criterion at hand, we ultimately aim for 


a morphological stability map. 


6.1 The Boundary Value Problem 


A protrusion of a metal electrode into a solid electrolyte stems either 
from an initial geometry (surface roughness), from inhomogeneous electro- 
chemical deposition or a combination of both. In this study, we prescribe 
a pre-existing protrusion and study the response of the interface kinetics 
without specifically identifying the source of deformation. For such a 
protrusion, we can conclude that if the current density at the tip is higher 
than a reference value on an undisturbed region, then the tip will grow 


faster than the surrounding, indicating an unstable system. 
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We consider a similar setup to Section [4.5] with a two dimensional geome- 
try of the solid electrolyte at plain strain, see Fig. (6.2h. The electrolyte 
has a thickness of L°'Y = 30 pm [33] and the width of the computational 
domain is chosen to be W = 100pm. While the top side of the separator 
is assumed to be held fixed, symmetry boundary conditions are applied on 
the sides. The bottom boundary is deformed with a Gaussian shaped pro- 
trusion. The displacement is upc = |X, APw(X/RP)], where we utilize 
the Wendland function of Eq. (4.4). We choose an amplitude AP = 4 um 
and radius RP = 15pm for the deformation. The metal electrode is 
assumed to be very stiff compared to the electrolyte and hence only the 
separator undergoes mechanical deformation. Furthermore, we know 
that lithium has a much higher electric conductivity compared to the 
ionic conductivity of the electrolyte and therefore we set the electrode 
electric potential as 6°’ = 0 along the bottom interface of the separator. 
Other than in Section [4.5] where we have assigned a fixed current density 
for the top and bottom interface, we now apply a constant current of 
JePP = 0.5 my only at the top of the solid electrolyte, see Fig. (6.2). and 
use the Butler-Volmer equation to determine the interface current density 
on the bottom electrolyte-electrode interface. The current density at 
the electrolyte-electrode interface is therefore the consequence of coupled 
bulk and interface transport described by the extended concentrated 
solution theory, see Eq. (3.116|) and reaction kinetics, see Eq. (5.69). 
All upcoming results correspond to the steady state solution. 

Polymer solid electrolytes are the main scope of this investigation. 
We hence utilize two recently published sets of parametrization, one 
for PEO [56], and one of the composite material PS-PEO [57], each 
with the salt LiTFSi. The latter host material has a higher Young’s 
modulus, dependent on the ratio of PS to PEO phase. The polystyrene 
phase (PS) enhances stiffness but also degrades the electrochemical 
transport properties since less ion conducting PEO material is available. 


Fig. [4.1] shows the concentration-dependent transport parameters. The 
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Figure 6.2: (a) Boundary value problem of the deformed solid electrolyte. (b) Color 
contour plots showing the salt concentration and the mechanical part of the hydrostatic 
pressure in the case of a PS-PEO electrolyte with E = 100 MPa. 


operating temperature is set to O = 353K. The interface conductivity 


is KBV = 0.589 2. We assume a Poisson ratio of v = 0.4 and will 


vary the Young’s modulus E. 


As in Chapter [4] the simulations are carried out with the finite element 
tool COMSOL Multiphysics® 5.3a. All field variables are carefully scaled 
to improve conditioning of the system matrix and thus gain appropriate 
numerical convergence. The non-linear equations of mass transport, 
balance of charge, mechanics and the non-linear boundary conditions 
by means of the Butler-Volmer equation are solved in a fully coupled, 
direct scheme. The outcome shows insensitivity to the width of the 
computational domain and to mesh resolution for the fine meshes used 


in the numerical analysis. 


6.2 Effect of Deformation on Interface 
Current Density 


Fig. shows the steady state solution for the deformed electrolyte 
during galvanostatic charging in terms of the salt concentration and the hy- 


drostatic pressure for one exemplary system (PS-PEO and E = 100 MPa). 
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Figure 6.3: The graphs show the response of the interface subjected to lithium 
protrusion for selected Young’s moduli: (a) Current density jPV; (b) Exchange 
current density jo. There is a qualitative change of the current density for increasing 
Young’s moduli suggesting morphological stability. 


The electro-chemo-mechanical coupling in the bulk is the superposition 
of two mechanisms. First there is a vertical salt-concentration gradient 
along the height of the separator due to the nature of the binary salt, see 
also the discussion in Section [4.4] Recall that the electric field acts on 
both the negatively charged anions and the positively charged cations. 
However, anions are restricted to the electrolyte and cannot leave or 
enter it. They then accumulate at the electrode where the cation enter 
the electrolyte and become dilute at the other such that, at steady state, 
diffusion equalizes the electrically driven transport. 

The second effect is mainly driven by the external deformation. As 
discussed in Chapter |4| a stress field builds up with a pressure peak 
close to the tip. The hydrostatic stress causes a depletion of the salt 
concentration due to the stress driven diffusion mechanism. The perimeter 
shows the opposite response. We observe tension on the shoulders due to 
the fact that the relatively stiff separator wants to follow the deformation 
of the tip, but is constraint by the given displacement of the boundary 
condition since we do not allow delamination. The increase in tension 
leads to tensile stresses as also predicted by Harry et al. [31] and yields 


an increase of the salt concentration. The findings show that the ion 
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Figure 6.4: The graphs show the response of the interface under protrusion for selected 
Young's moduli: (a) mechanical contribution of hydrostatic pressure, (b) the surface 
stress, (c) salt concentration, (d) open circuit potential, (e) electrochemical potential 
in the electrolyte, (f) overpotential. 


transport in the electrolyte is also strongly coupled with interface kinetics, 
especially in the case of inhomogeneous deformation. 

We now consider PS-PEO parameter values for the transport properties 
and three Young's moduli, namely E = (100,300, 600}MPa. Fig. [6.3h 
shows the current density at the interface for the three stiffnesses at 
steady state. For a relatively small elastic modulus (E = 100 MPa), a 
two times higher current density can be observed at the protrusion peak 
compared to the flatter region. This changes when the Young’s modulus 


is increased leading to a current density at the protrusion tip that is 
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almost equal to (E = 300 MPa) and smaller than (E = 600 MPa) the 
current density in the flatter, undisturbed region. With E = 100 MPa the 
simulation shows a faster deposition process at the tip which therefore 
means an unstable surface. The opposite occurs for E = 600 MPa. The 
qualitative change of current density for different Young’s moduli is an 
indication of morphological stability and we will rely on this idea for the 
upcoming analysis. Before proceeding, we wish to discuss the origin of the 


change of current density distribution when the Young’s modulus is varied. 


For the analysis of the current density at the electrolyte-electrode 


interface we recall the Butler-Volmer equation from Eq. (5.69) 


PY = jo ep (a 0751) exp ( 575) (6.1) 


as a function of the exchange current density jo and the overpotential 


n= —Lhwy Va with Y? HQmosa, The interface current density is 
therefore dependent on the electro-chemo-mechanical state of the cation 
in the solid electrolyte and the stress state of the electrode. In the 
following, we systematically investigate the different contributions. 

We first focus on the bulk transport in the solid electrolyte. The 
protrusion leads to a positive hydrostatic pressure pM at the protrusion 
tip as seen in Fig. [6-4h. The hydrostatic stress increases for higher 
mechanical stiffness. We also observe a state of tension at the perimeter 
of the protrusion. Fig. [6.45 shows the spatial salt concentration at the 
interface for the three investigated Young’s moduli. There is a depletion 
of salt concentration at the protrusion peak due to the effect of stress 
driven diffusion as it is also observed in Fig. [6.2p. Recall that the driving 
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force of the electrolyte for the reaction at the interface is the cationic 


electrochemical potential 


wy= u? + RO Mn (=) + ROln (4) Jeo + PFS, (62) 


which has been derived in Eq. (3.115) and (3.77). This field variable 
is plotted in Fig. and stems from a superposition of mechanical 
(Fig. [6.4h) and entropic effects (Fig. 6-4). We see that the mechanical 
contribution dominates. For a change of Young’s modulus, we observe 
that w./F increases from 20mV (for E = 100 MPa) and up to 45mV 
(for E = 600 MPa) at the protrusion tip compared to the flatter region 
with a constant value of roughly 7mV. 

We now investigate mechanical influences on the open circuit potential 
Vo. and exchange current density which enters the reaction kinetics. 
Fig. shows the normal stress at the interface for the three different 
Young’s moduli. We observe a compression at the protrusion peak which 
is more pronounced for higher stiffnesses, but also some tension at the 
perimeter. The form is similar to the hydrostatic pressure in the solid 
electrolyte (Fig. |6.4h), but the absolute values are a bit higher. The 
open circuit potential V$? scales linearly with the normal stress due to 
Eq. (5.52), which yields the values of V§{ at the interface in Fig. [6.4H. 
While the open circuit potential is close to zero at the undeformed regions, 
we see values of V{9 from —7mV (for E = 100 MPa) up to —46mV (for 
E = 600 MPa). Positive values of V9 are present at the perimeter. 

The overpotential 7 at the interface for different Young’s moduli as 
seen in Fig. is computed using Eq. from the electrochemical 
potential (Fig. [6.4fl) and the open circuit potential (Fig. [6.4p). Recall 
that the electric potential of the electrode is assumed to be constant 
along the interface. The surface overpotential as shown in Fig. is 
strongly negative (anodic reaction) at the tip for E = 100MPa and 
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positive (cathodic reaction) for E = 600 MPa. The intermediate modulus 
E = 300 MPa shows some oscillation but is overall fairly uniform over 
the protrusion surface. 

The second variable that is necessary to compute the current density 
and that is influenced by mechanics is the exchange current density, see 
Eq. (5.69). The mechanical stresses at the tip of the protrusion cause an 
increase in jo as seen in Fig. [6.3p. The exchange current density at the tip 
is a factor of 1.5 higher than that in the flatter region for Æ = 600 MPa. 
In contrast to the overpotential, which the normal stress and electrolyte 
bulk stress enter with opposite signs, the mechanical influences add up 
together and therefore increase jo as the Young’s modulus increases. Note 
that there is a small dip for E = 600 MPa due to a strong depletion of 
the salt concentration. 

Ultimately, the distribution of the current density in Fig. broadly 
follows that of the overpotential y in Fig. |6.4f. The increase of the 
exchange current density in Fig. |6.3p with the increase of stress is 
especially pronounced for materials that are quite stiff. The current 

mA 


density at the tip is at —1 en: approximately twice as high in magnitude 


as the applied current density for E = 100MPa and much less for 
E = 300 MPa. For E = 600 MPa, it even changes from deposition to 
depletion due to high mechanical stresses. The exchange current density 
has only a small impact on the overall current density in our example. 
The current density for materials with higher Young’s modulus decreases 
at the protrusion tip and thus aligns with experimental findings of Harry 
et al. Bij. Another aspect observed by the experimental study is an 
increase of current density on the perimeter of the protrusion. Our 
simulations also show this effect for a high Young’s modulus material 
(Fig. [6.3h). 

To summarize the analysis, we see that a deformation of the solid 
electrolyte triggers stabilizing and destabilizing effects for the growth of 


a protrusion. Three contributions are competing. A positive hydrostatic 
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pressure p™ increases the electrochemical potential of the electrolyte, 
thus promotes a faster deposition and is therefore destabilizing. The 
effect of stress driven diffusion yields a depletion of salt concentration 
which lowers the electrochemical potential of the electrolyte as well as the 
exchange current density and thus decreases the current density at the 
protrusion tip (stabilizing). The third effect is the open circuit potential, 
which eventually lowers the current density at the peak of the protrusion 
for higher normal stresses and is also stabilizing. Depending on the trans- 
port parameter and mechanical properties, either the stabilizing or the 
destabilizing effects dominate. Therefore, we will develop a morphological 


stability map dependent on these factors in the next section. 


6.3 Morphological Stability Map 


The findings of the previous section, which showed a qualitative change of 
current density along the deformed interface due to a variation of Young’s 
modulus of the electrolyte, suggest a morphological stability criterion for 
solid polymer electrolytes. The idea is that there is a critical Young’s 
modulus, for which the current density at the peak of the protrusion 
becomes smaller than that on the undeformed interface. Our aim is to 
find this critical modulus as a function of material properties, geometry, 


interface conductivity and applied current density. 


6.3.1 Definition of Morphological Stability 


To avoid influence of the small oscillation in results between the peak of 


the protrusion (see e.g. Fig. |6.4h), we use a weighted average current 
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Figure 6.5: The current density evaluated at the side, the tip and with the smoothing 
scheme of Eq. for (a) PEO and (b) PS-PEO. The crossing point between 
evaluation at the side (dashed line) and smoothed evaluation scheme (straight line) 
determines the critical Young's modulus and is in the same region compared to the 
tip evaluation (dotted line). 


density over the protrusion instead of a single value and compare it to 


the current density at the side of the computational domain. This yields 


wy2 
de T JPY (X) w (X/R?) de | ia >0 unstable 
Dem (X/RP) da <0 stable, 
(6.3) 


where the Wendland function w is defined by Eq. (4.4) as for the 
protrusion itself. The contribution of the current density adjacent to the 
peak of the protrusion is therefore smoothed out, and the current density 
in the regions where it oscillates is given less weight. 

Our intent is to find the smallest Young’s modulus where Eq. (6.3) 
changes from a stable to an unstable interface. We denote this value as 
the critical Young’s modulus Et and we will see that this value strongly 
depends on the properties of both the bulk material and the interface. 

The concept of the stability criterion of Eq. (6.3) is thereby similar 


to the one used by Monroe and Newman [72] who compare current 
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density at the peak of the protrusion to the current density at the lowest 
point in the valleys. Fig. shows the difference of the two evaluation 
schemes for PEO and PS-PEO for the setup herein. The current densities 
are evaluated at the side ¡BY (0.5W), at the tip ¡BY (0) and with the 
smoothing scheme of Eq. (6.3) as a function of Young’s modulus. The 
values from the side remain almost constant over the Young’s moduli 
range and thus serve as a reference in the stability criterion (6.3). The 
current density at the tip crosses the reference in a similar range of the 
Young’s modulus compared to the values obtained by the smoothing 
scheme. This is 14MPa (smoothed) and 20MPa (tip) for PEO and 
230 MPa (smoothed) and 310 MPa (tip) for PS-PEO. Since protrusion 
growth is usually not a local effect but a process along the perturbed 
interface (see Fig. [1.3p) we find the smoothed scheme in stability criterion 
(6.3) more appropriate than the tip evaluation. Nevertheless, we want to 
emphasize that this is not conservative since the tip evaluation indicates 


up to 40% higher requirements to the electrolyte stiffness. 


6.3.2 Influence of the Protrusion Height 


First, we investigate the influence of the amplitude of the protrusion on 
the critical Young’s modulus. The stability-instability boundaries for 
PEO and PS-PEO are shown in Fig. as a function of the protrusion 
amplitude AP. This boundary is the locus of points for which the term in 
parenthesis in Eq. is equal to zero. Above the lines, the interfaces 
are predicted to be stable for the respective materials and below to be 
unstable in the sense that the height of the protrusion will grow. We 
varied the amplitude between 0.1 um and 4pm. The results in Fig. 
show that the critical Young’s moduli for stability are dependent on the 
amplitude of the protrusion. However, the critical Young’s modulus does 
not change much for both PEO, where it varies from 13 to 28 MPa, and 
PS-PEO, where it ranges from 246 to 330 MPa. Although the results of 
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Figure 6.6: (a) The critical Young’s modulus separating stability from instability de- 
creases only slightly as the amplitude of the protrusion is increased. (b) Morphological 
stability map for protrusion amplitude AP = 4pm. The critical Young’s modulus 
(blue line) is plotted as a function of a measure of how good the transport properties 
of a material are. Above the blue line there is a stable electrolyte in the sense that 
the protrusion is predicted not to grow. Below the blue line the electrolyte is unstable 
and the protrusion will grow. 


the critical Young’s modulus should be understood more qualitatively, the 
order of magnitude corresponds quite well to experimental findings. The 
polymer electrolyte PS-PEO with a Young’s modulus of approximately 
150 MPa [48] showed stable behavior with lithium metal anodes. On the 
other hand, the soft polymer electrolyte PEO that is usually operated 
above the glass transition temperature has a high likelihood for protrusion 


growth. 


6.3.3 Variation between PEO and PS-PEO 


We have already seen that the PEO material with relatively good trans- 
port properties (see Fig. requires a much lower stiffness to ensure 
stability compared to the PS-PEO material, see Fig. [6.6h. To broaden 
our insights and to more generally consider the trade-off between interface 


stability, elastic modulus and lithium transport, we interpolate between 
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the sets of transport properties for PEO and PS-PEO by use of the pa- 
rameter €, that is equal to 0 for PS-PEO and 1 for PEO. This parameter 
can be thought of as charactering a family of composites with different 
amounts of PS present compared to PEO. With such a parameter we 
interpolate between the electrochemical properties of PEO and PS-PEO 
by use of 


Y = EYPEO + (1 — £) rPS-PEO, ¿€ e [0,1], (6.4) 


where € is the interpolation parameter and Y € {K, D, t4, Bln falo), en 


are the transport parameter. We carry out a series of stability simulations 
using the values of transport parameters produced by various values of 
€ in Eq. (6.4). We assume y = 0.4 for both materials throughout this 
analysis. 

The results are shown in Fig. |6.6þ where the critical elastic modulus 
is plotted against the interpolation parameter £. The blue line indicates 
the boundary between stability and instability for the protrusion growth, 
see Eq. (6.3). As shown in Fig. [6.6p, the stability regime is above the 
blue line and the unstable regime is below it. We note that there is a 
monotonic decay in the value of the critical Young's modulus required 
for stability of the protrusion as one ranges through composites from 
PS-PEO to PEO. 


6.3.4 Variation of Selected Parameter for PS-PEO 


We now study the correlation of the critical Young's modulus with respect 
to selected parameters. We use PS-PEO as reference material system 
and modify only one parameter per each analysis. Fig. shows the 
correlation of conductivity and critical Young's modulus Et. We see 
that a higher conductivity is beneficial for stability. Especially a decrease 


of the ionic conductivity has a strong impact, e.g. if it is two orders 
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Figure 6.7: Parameter study of the critical Young’s moduli with PS-PEO as reference 
material. The following parameters are varied separately while keeping all others 
fixed: (a) ionic conductivity; (b) diffusivity; (c) transference number; (d) interface 
conductivity; (e) separator thickness; (f) applied current. 


of magnitude lower than PS-PEO then a Young’s modulus of up to 
E = 1200 MPa is required. If the diffusivity is increased, the value of 
E“it drops significantly, see Fig. [6.7p. A strong benefit with respect 
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Figure 6.8: The graphs show the critical Young’s moduli for stability of the protrusion 
as a function of € from Eq. (6.4), a measure of how good the transport properties 
of a material are. (a) Results are shown for three charging rates indicating that 
the protrusion is less stable at higher charging rate. (b) Results are shown for two 
non-uniform and a uniform value of KPV. As can be seen, non-uniformity causes the 
protrusion to be less stable because it enhances the interface conductivity at the peak. 


to morphological stability is already obtained for the relatively small 
increase of diffusivity by a factor of five. Smaller diffusivities lead to a 
complete depletion of salt and thus a termination of the simulation. Fig. 
6.7¢ shows the impact of the transference number. Note that we set the 
transference number constant in this study. For tf" = 0 and t° = 1 
we obtain Et = 400 MPa and 21 MPa, respectively. The intermediate 
values scale almost linearly. Next, we study the impact of the interface 
conductivity which enters the simulation model in Eq. (5.69). Fig. |6.74 
indicates that degrading the interface conductivity is harmful for stability 
since a higher E is required. However, there is a saturation limit, which 
is in the case of PS-PEO at about 240 MPa. The relation of Æt to the 
thickness of the separator is presented in Fig. [6.7p. We see that a very 
thin separator is beneficial for the morphological stability|"] Values higher 
than LY = 50 um, however, are equally stable with Et ~ 300 MPa. 


1 Note that we assume a stiff cathode. However, for very thin separators, a deformation 
of the cathode is more likely and thus should resolve to increase accuracy. 
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Next, we study the influence of current density. Experimental findings 
and theoretical model formulations indicate that a higher charging rate 


supports the growth of an intrusion [32]. The result of a corresponding 


mA 
cm?’ 


analysis is shown in Fig. [6.8p. Three sets of calculations, with 0.1 
0.2 ma, and 0.5 25, are carried out for composites with a range of 
transport parameters as characterized in Eq. (6.4). Fig. indicates 
that higher charging rates demand higher Young's moduli. Furthermore, 
we see in Fig. a linear correlation between applied current and critical 
Young's modulus for PS-PEO. This outcome holds also for PEO and the 
intermediate values, see Fig. [6.8h. 


6.3.5 Heterogeneous Interface Conductivity 


One explanation for the source of protrusions is a heterogeneous interface 
conductivity, e.g. triggered by the presence of an inhomogeneous SEI 
or defects in the metal electrode [50], see the discussion in Chapter 
These inhomogeneities are assumed to decrease activation energies and 
therefore to increase the exchange current density locally. We model 
this phenomenon by applying a factor to the standard exchange current 
density and use the Gaussian like function w(r) defined in Eq. as 


the perturbation. The modified interface parametrization then reads 
KPY — (1+ \w(X/R?)) KPY, (6.5) 


where X defines the amplitude of heterogeneity and X is the horizontal 
coordinate of the interface. Therefore, e.g. for A = 1, the region of 
highest protrusion has an interface conductivity twice that prevailing in 
the valleys on either side of the protrusion. 

Simulations are carried out at an applied charging rate of 0.5 ma for 
A=0 (uniform), A = 1 and A = 2 and the results are shown in Fig. [6.8p. 


The results are plotted versus € and show that the non-uniform interface 
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conductivity causes the protrusion to be less stable and more likely to 
grow. This occurs because the non-uniform interface conductivity has a 
higher value at the peak of the protrusion and therefore more lithium 
is deposited there than on the shoulders adjacent. PEO demands an 
almost ten times higher stiffness for stability of the protrusion when 
the non-uniformity in the interface conductivity is present, whereas the 
necessary increase of Et for PS-PEO is not as pronounced. For PS-PEO 
the critical value of Young’s modulus that ensures protrusion stability 
increases only by ~ 20% to roughly 300 MPa for A = 2. Fig. shows 
also, that the requirement on the stiffness is overall similar for A = 1 and 
A= 2. 


6.3.6 Discussion 


The morphological stability analysis of the deformed solid electrolyte 
with the fully coupled electro-chemo-mechanical model indicates a strong 
correlation of the material stiffness and local deposition current density 
while protrusion growth. Although we assume the worst case scenario, 
an infinitely stiff metal electrode, a growth of the protrusion can be 
suppressed if the electrode material has a certain mechanical stiffness. 
This is contrary to analyses such as that of Monroe and Newman [72], 
who find a relationship between material properties of the metal electrode 
and solid electrolyte. It further contradicts the simple purely mechanical 
consideration that the electrolyte has to provide a certain stiffness such 
that the lithium protrusion is suppressed by "pushing it down". 

By using a full continuum description of transport and interface, we are 
able to show that the electrolyte transport properties correlate strongly 
with the stability of the interface. The better the transport properties are, 
the lower are the requirements on a "stable" separator’s Young’s modulus. 
The stability map from Fig. indicates that even a mechanically soft 


material such as PEO has almost the capabilities to suppress protrusion 
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growth. Material modifications like cross-linking or using a 
gel electrolyte with a double polymer network (E = 44 MPa [42]) are 
possible paths to successfully suppress protrusion growth. Of course, these 
materials cannot be directly associated to a specific € in the interpolation 
scheme of Eq. (6.4), but they are most certainly closer to PEO and thus 
do not require as high stiffness as PS-PEO does. The critical Young’s 
modulus for PS-PEO is in the same order of magnitude compared to its 
measured mechanical properties and therefore aligns with experimental 
findings of stable deposition. If the glass transition temperature of PS is 
exceeded, then PS-PEO looses its stiffness and therefore its capabilities 
to suppress the growth of protrusions [172] . 

Furthermore, we varied several parameter of PS-PEO as shown in 
Fig. |6.7| The outcome confirms the correlation of transport properties 
and stability: If the transport parameters of PS-PEO are degenerated, 
then the requirements on the separator are higher. We see an almost 
linear correlation between the applied current density J*PP and E**, 
Note that with the values of J®PP in our simulations we are below the 
limiting current density of approximately 2.7 mo for PS-PEO as predicted 
e.g. by [30] 32] [77] who correlate dendrite growth with a full depletion 
of salt at the SEI. A pronounced likelihood of a complete depletion of 
salt at the SEI is given by a lower ionic conductivity, diffusivity and 
transference number. Fig. [6.7p-c agrees with this understanding as it 
shows that a higher stiffness is required when the transport parameters 
are degraded. A higher transference number showed to be beneficial for 
stability as already suggested by Tikekar et al. [89]. The interface 
conductivity also has a strong influence on stability, especially for small 
interface conductivities, which aligns with the observations of Liu et al. 
[173]. A thinner separator is also beneficial for stability. The trend of 
higher requirements on the separator stiffness for higher current densities, 
lower transport properties and larger separator width align with the idea 
of a depletion triggered dendrite growth [77]. 
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If we solely consider the analysis of an interface with homogeneous 
transport properties, we can conclude that increasing the stiffness by 
sacrificing good transport properties is not necessarily helpful for the 
morphological stability. Obviously, this is good news in fulfilling the 
demands of high charging rate performances for batteries since a material 
design that optimizes transport properties can also be beneficial for 
morphological stability. 

Furthermore, it is worth to consider the inhomogeneities of the interface 
with a local increase of interface conductivity. We then see that the 
mechanically weak PEO is far away from the stable region. This is 
different for PS-PEO where the impact of locally increased interface 
conductivity is relatively small. The finding highlights the importance 
of homogeneous interfaces, especially for mechanically weak separator 
materials. We have further seen an almost linear correlation between 
charging rate and critical Young’s modulus no matter which material. 
This is a problematic feature of polymer electrolyte battery systems 


especially for fast charging applications. 


6.4 Summary 


This chapter discussed a quasi-static analysis of morphological stability 
of a metal electrode in contact to a binary solid polymer electrolyte 
based on a rigorous electro-chemo-mechanical description of the bulk 
and reaction kinetics. We utilized measured transport parameters of the 
polymer electrolytes PEO and PS-PEO. Analysis of the current density 
at a deformed interface showed that the normal stress dependency of the 
equilibrium potential, the mechanical influence onto the electrochemical 
potential in the electrolyte bulk as well as the value of the exchange 
current density are equally relevant. The current density at the protrusion 


tip showed a clear dependency on the Young’s modulus of the electrolyte. 
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This effect is present even though the metal electrode is assumed infinitely 
stiff - which is in agreement with the fact that the Young’s modulus 
of a polymer electrolyte is at least one order of magnitude lower than 
pure lithium. The simulation is able to quantitatively reproduce the 
experimental findings of Harry et al. [31], i.e. the increase of current 
density at the protrusion perimeter. 

The analysis showed a correlation between required stiffness and trans- 
port properties. PEO demands only a fraction of the critical Young’s 
modulus compared to PS-PEO electrolyte, yielding in a plausible range al- 
though some drastic model assumptions were made (i.e. elastic behavior). 
The study suggests that it is important to find the right trade-off between 
the electrochemical and mechanical properties of the solid electrolyte. 
It is not necessarily better to increase the mechanical stiffness if this 
increases is accompanied by a decrease of transport properties. 

We further find that the charging rates and local modifications of in- 
terface conductivities alter the morphological stability plot. The required 
stiffness scales linearly with the charging rate and therefore suggests 
upper bounds for the charging rate for a given Young’s modulus for 
fast charging applications. Local defects experimentally showed a strong 
impact on the likelihood of protrusion growth. Again, this could be ratio- 
nalized by the presented calculations. Especially electrolytes with good 
transport properties such as pure PEO require a much higher mechanical 
stiffness if interface inhomogeneities are present. 

The presented methodology is therefore an attractive tool to predict 
the "sweet spot" of a solid polymer electrolyte which has sufficient 
transport properties by obeying at the same time the mechanical stiffness 


requirements for morphological stable lithium/polymer interfaces. 


159 


7 The Metal Electrode 


Metal electrodes, in particular lithium metal electrodes, and their usage 
within batteries are an interesting modeling subject for several reasons. 
The first and most obvious property of metal electrodes is their ability 
to conduct electrons. Secondly and most importantly, the oxidation 
and reduction process of cations (i.e. lithium) at the interface to the 
electrolyte leads to a deposition and stripping of material. The mass and 
therefore the stress free volume of the electrode thus changes due to an 
applied current. The third modeling subject is the mechanical response 
of lithium itself. 

This chapter focuses on the kinematics of material deposition and strip- 
ping. This includes the volume description, the interface compatibility 
and the growth rate. We will also briefly discuss the material properties 


of lithium metal such as elasticity and electric conductivity. 


7.1 Kinematics of Material Deposition 


The stress free growth of lithium during deposition (charge) and shrinkage 
during stripping (discharge) is associated to the current density at the 
metal electrode-solid electrolyte interface. Assuming no significant surface 
migration of ions [174], one can deduce the boundary velocity [| wder || in 
analogy to Ferrese and Newman at a specific point by counting ions 


passing the interface. We obtain 


Ni -de 
eter = Mer, (7) 
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Figure 7.1: Illustration of the two mappings necessary to describe deposition or 


stripping and deformation: The geometry frame describes the metal electrode at initial 
etr 


state. Deposition is modeled with the mapping ¢f* which yields the undeformed 


material frame. With the mappings of the electrolyte po and the electrode y$"", we 


obtain the deformed state, where the interface shapes S; of electrode and electrolyte 
coincide. 


where jP is the current density at the interface i.e. obtained by the 
Butler-Volmer equation (5.69), Qui is the molar volume of lithium atoms 
in the lattice structure and F is the Faraday constant. 

Eq. (7.1) is usually utilized in liquid electrolyte systems [58] [79], 
where the liquid solvent can be easily displaced during inhomogeneous 
deposition. In contact with a solid electrolyte, however, the newly de- 
posited material on the metal electrode is able to generate mechanical 
stresses in both electrode and electrolyte, which will lead to local defor- 
mations. Therefore, we describe the dynamics of deposition and stripping 
at the metal electrode together with the mechanical coupling to the solid 


electrolyte, and we do so within the concept of finite strains. 


7.1.1 Volume Description 


In Chapter [2] we have introduced the idea of finite strain continuum me- 
chanics that distinguishes between a material and a spatial configuration 


due to deformation. Recall that the process of deformation underlies 
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the conservation of mass. Deposition and stripping, however, change the 
mass of a given, undeformed (material) electrode domain B£t" and are 
therefore dependent on time t € T', see Fig. [7.1] To depict the temporal 
dependency, we introduce a third configuration Bet independent of time 
and which we refer to as geometry configuration. This domain might 
be associated to the initial shape of the metal electrode. Each point 
X, € Bs can be uniquely mapped to X € Bft" by 


de Bs x T — Bet 
ne Pe (7.2) 
(Xat) > X= Kt). 


The geometry displacement field d = X, — X describes the geometry 
change in the metal electrode due to deposition or removal of mass. In 
addition, we characterize the mechanical deformation by the deformation 


map 
ptt = . (7.3) 


which is equivalent to Eq. (2.1). The displacement field u = X — g is 
determined by the balance equations, the constitutive material model 
and the boundary condition. 

Here we have to turn to the numerical treatment to fully describe the 
methodology used. Describing only a deformation process, the associated 
discretization (grid points) in the interior of the domain is usually associ- 
ated to the material body description and therefore are displaced in a 
physically meaningful way by the deformation map $t". This is different 
for the deposition and stripping process ¢é', where the material in the 
interior and the associated discretization points are not physically altered. 
The newly generated material is therefore not adequately resolved by the 


discretization and, even worse, mesh entities can disappear if the removal 
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of material is greater than the mesh size. To avoid this problem, we 


equally redistribute the mesh in the interior domain y$ with a Laplace 


smoothing [132] in every time step 


Divx, (Gradx, (X)) =0, (7.4) 


where the operators Divx, (-) = He), and Grad xy (-) = wets 


defined with respect to the geometric coordinates X,. Note that the 


are 


choice of the Laplace operator is not physically motivated and therefore 
not the only option. Since it showed reasonable results for our application 
we will use Eq. henceforth] 

Solving the Laplace smoothing of Eq. (7.4) requires appropriate bound- 
ary conditions. The geometric boundary displacement d is zero for 
surfaces which are not affected by material deposition and stripping. The 
other boundaries are determined by the growth rate and will be discussed 


in the upcoming section. 


7.1.2 Interface Compatibility 


It is obvious that electrochemical reactions require interfacial contact 
between metal electrode and solid electrolyte. If an interface region 
delaminates, the reduction and oxidation reactions are most likely 
blocked and no further deposition and stripping of material can take 
place. It was shown that polymer electrolytes are able to maintain 
adhesion to the lithium electrode [48]. This finding is supported by X-ray 
microtomography images [50] which do not show delamination effects of 
PS-PEO electrolyte and lithium metal at a protrusion while depositing. 
The same holds for PEO [30]. What is not known from these image 


1 Other smoothing methods provided by COMSOL Multiphysics® are computationally 
more expensive and have not shown an improvement of the smoothing quality 
Multiphysics [132]. 
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techniques is whether a tangential slip along the interface occurs during 
deposition and stripping or not. Measurements of friction coefficients 
have not been conducted yet, and therefore, we will describe two extreme 
situations, one where the interfaces are perfectly glued together and the 
second where tangential slip is not constrained. 

The kinematics of a perfectly aligned interface are as follows. We 


denote the interface of the metal electrode and solid electrolyte as 
Si = 91(BSY) N p(B"), (7.5) 


given in the spatial configuration (Fig. [7-1right). Recall Section [3.1] for 
the kinematics of the electrolyte such as the deformation map pS”. The 
interface S; is either computed from the electrode S; = pf" (Cf (Sgt) 
or from the electrolyte side S; = y?” (a): The shape of the interface 
is therefore determined by three contributions: the deformation in the 
electrolyte pi”, the deformation in the electrode $*™" as well as the 
geometry change Cet". 

Note that the interfaces of electrode and electrolyte do not align in the 


material configuration (Fig. [7.1] center): 


(pi) (St) FA). (7.6) 


We therefore have to apply special treatment for the growth of material 
since the material boundary conditions defined on the electrolyte (such as 
the current density) have to be converted to obtain consistent boundary 
conditions for the metal electrode. 

To obtain the perfect contact given by Eq. (7.5), we have to couple 
the mechanical response of the electrode-electrolyte interface. Let’s 
assume that the interfaces align in the initial state, i.e. S¿% = Ses The 
assumption of perfect compatibility is then fulfilled if and only if the 


boundary displacement of each interface point of the electrolyte u%Y is 
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equal to the combined geometry change d° and the deformation of the 


etr 


metal electrode u“. To couple both, we demand pointwise 


Ul sev = dl sg + User, (7.7) 


which is also used by Tikekar et al. [89]. Note that the compatability 
condition of Eq. does not allow any tangential slip or in other words, 
the boundaries are glued together. 

Alternatively, one obtains a transition condition with perfect contact 
together with tangential slip by imposing equal boundary velocities in 
the spatial normal direction to the interface. The boundary condition 


then reads for each point on the interface S; 


(vn) ef =—(v-n) Sr (7.8) 


with the electrolyte velocity JE = ap (X), the electrode velocity 


vš = E (eps (CE (X,))) and spatial normal vector n that points in 


opposite direction for the electrode and electrolyte. 


7.1.3 Growth Rate 


We now derive the relation of current density and interface velocity taking 
the geometric changes of the interface during deposition and stripping 
into account and therefore generalize Eq. (7.1). The current density ¡BW 
is computed with respect to the deformed geometry with area element da 
and normal direction n whereas the deposition is defined on the material 


configuration with area element dA and normal direction N. 
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We make use of the fundamental law of conservation of mass to resolve 
the growth rate of the interface between electrolyte and electrode. It is 
given by 


d gor 
—mri+ | Mi 
Ge dee 


da =0 (7.9) 


and states that the rate of change of (lithium) metal electrode mass my; 
is equivalent to the material with molar mass Mg; deposited or stripped 
due to the current density j*°?. Note that the current density is defined 
on the deformed interface as it specifies the active area for electrochemical 
reactions. The mass of the electrode mz; can be written in terms of the 


density, yielding 


d d / 
—ML = > 7 pLidv (7.10) 
dt dt Pr (Betr) 


with the mass density pri. In order to exchange the integral and the 
temporal derivative, we switch from the spatial via the material to 
the geometry configuration and thus use the time-independent integral 
domain B§"". With the help of Reynold’s theorem (Eq. (2.14)), we take 


the moving boundary due to deposition into account and obtain 


d d o d , f a 
LS dV = < o;dV eter. NdA, (7.11 
dt E dt I... PLid Je dt PLi g Set PL? a , i ) 


Na, pm) 
=0 
where the first term on the right hand side drops out due the balance of 
: d 
mass of an undeformed medium, see Eq. (3.6) ‚and vP = 2X, | sere 
is the unknown interface velocity for deposition and stripping defined 
in the material configuration. Note that this velocity is different to the 
one in Eq. (7.8) because it takes only deposition but no deformation 
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into account. Using Eq. (2.8) we obtain the spatial representation of the 


boundary integral 


d 
—myi = / priv’? . F’ nda (7.12) 
dt S: 

Inserting Eq. (7.12) into Eq. (7.9) and postulating that material deposi- 


tion and stripping are local effects gives 


Myj@er 


= — ppv? . FTn. 7.13 
F PLiY n ( ) 


We are now interested in the normal velocity of the boundary in the 
material configuration. We thus multiply Eq. (7.13) with the material 
normal N and get with N: Fn = F: (n & N) the normal velocity as 


Mui jap _ _ Mui J 


dep . N = 
j A FA; F:(n@N 


gee (7.14) 
) 
where the scaling factor between JA and j is taken from Eq. (2.10). 
Note that all properties in Eq. (7.14) such as the deformation gradient 


F are associated with the metal electrode. 


7.2 Material Properties of Lithium Metal 


In the past, mechanical response of pure lithium has received little 
attention due to its very reactive nature and limited fields of application. 
Studies on anti-proton sources [11] [2] provided first measurements on 
mechanical stiffness with E = 1.9—4GPa and non-elastic behavior 
such as the mechanical strength with about 1 MPa [12]. Due to the rise of 
solid state batteries with lithium metal anodes, the material properties of 
lithium got into the focus of research. Microindentation [15], compression 
experiments [14] and nanoindentation [175.177] showed higher yield 
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strength with up to 105 MPa, size effects and anisotropy due to the grain 
structure of lithium metal. Due to its low melting point of 180.6°C and 
for an application at elevated temperatures (e.g. 80°C) such as necessary 
for polymer battery systems [50], there is also a high likelihood that the 
material starts to creep. This might influence long term applications 
in battery systems. Recently, Narayan and Anand [178] published an 
elastic-viscoplastic material model which describes the complex response 
of lithium. 

In this work, however, we put our focus on the electrochemical response 
with a simple mechanical description. Therefore, we apply the same level 
of complexity as we have used for the solid electrolyte in Chapter [3] and 
assume purely elastic behavior. We also note that by omitting non-elastic 
effects, we model a material which is more prone against deformation. 
This assumption is therefore conservative with respect to a morphological 


interface stability analysis. 


7.2.1 Material Model 


Despite the knowledge of non-elastic effects of lithium as stated before, 
we reduce complexity by the assumption that the full deformation as de- 
scribed by the deformation gradient FF is purely elastic. This assumption 
is especially true for small strains. Similar to the solid electrolyte, we 
assume an isotropic hyperelastic Neo-Hookean material model [111], but 
without considering swelling effects (see Eq. (3.66)). The strain energy 
is then given by 


~ 1 1 
Na (PF) = 3 NH (I, — 3) + zANH (log 13)” — ynu log Is, (7.15) 
with the Lamé constants ynu, Any and the invariants 


h=tr (FTF) y k=J=detF. (7.16) 


169 


7 The Metal Electrode 


The elastic material model helps us to understand the basic electro-chemo- 
mechanical response and is in analogy to the works of e.g. [151 [71] [86]. 
This approach of using solely an elastic model also circumvents questions 
whether newly deposited material "inherits" some history of the pre- 
existing material such as strain hardening, etc. Considering in-elasticity 
also requires a mapping of history variables such as plastic strain or 
hardening parameters to the geometric configuration and not to the 


material configuration as it is usually done. 


7.2.2 Ohm’s law 


Electronic conduction in a metal electrode is governed by the balance of 
charge where the current density is deduced by Ohm’s law. It states a 
linear relationship between electric current density Jq and voltage drop. 
Recall Eq. for the differential form 


div(j,)=0 with j, =—K°grad(®), (7.17) 


etr is equal to the 


where the conductivity of the metal electrode is k 
inverse of the resistivity. Lithium metal, for example, has a conductivity 
of Kt" = 10.6 x 10° 31 [179] which is almost ten orders of magnitude 
higher compared to the one of a polymer electrolyte (see Section [4.2.1) 
since the small electrons can more easily be transported compared to the 
bulky ions of the salt. Hence, small gradients of the electric potential 
already yield a very high electric current. Since the length of the electron 
path between current collector and electrolyte is small, roughly = 100 pm, 
the drop in electric potential in the electrode is expected to be small. 
For the consistent implementation of Ohm’s law on a deformed geome- 


try, we again use the formulation of finite strains and solve the partial 
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differential equation in the material configuration] The electric field in 


the metal electrode then is computed by 
Div(J,)=0 with J,=-JIFIR"FTGrad(®), (7.18) 


where the electric conductivity is transformed in analogy to Eq. (3.84) 
into the material configuration. Appropriate boundary conditions such 


as Cauchy or Neumann condition, respectively, are 
Plage = 6, — —N-Jgloge = Jar. (7.19) 


The boundary conditions are defined in the undeformed geometry and 
we apply Eq. (2.10) to convert the current density of deposition J4°P to 


j*P to take the deformation of the interface length into account. 


7.3 Summary 


This chapter dealt with the description of the metal electrode in a solid 
state battery. Two mechanisms lead to a change of geometry. One process 
is related to deposition and stripping, which adds and removes material 
from the electrode. Mechanical deformation of the electrode is the second 
process that has to be modeled. We derived a consistent formulation 
within the finite strain framework by utilizing in addition to the commonly 
used displacement field (deformation) a geometric displacement field 
(for deposition and stripping). We postulated interface compatibility 
between electrode and electrolyte motivated by experimental findings. 
Furthermore, we presented material properties of lithium metal concerning 


elastic response and its capability to conduct electrons. 


2 We use the capabilities of COMSOL Multiphysics® 5.3a which automatically maps 
the material to geometric configurations correctly if no history variables are involved. 
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This chapter deals with a second approach to describe the morphological 
stability of a solid electrolyte - metal electrode interface. In addition 
to the scope of the first study in Chapter [6] we now model the metal 
electrode and the growth of material during charging explicitly. 

The study utilizes the electro-chemo-mechanical model of ion transport 
in a solid polymer electrolyte (Chapter [3), the extended Butler-Volmer 
reaction kinetics by means of mechanics (Chapter [5) and describes explic- 
itly the lithium electrode and its growth and shrinkage due to deposition 
and stripping (Chapter [7). 

One of the few explicit treatments of lithium metal growth in battery 
applications on the continuum scale was presented by Ferrese et al. 
[58] who describe the ion transport by using the classical concentrated 
solution theory [53] but disregard any mechanical influences. The follow- 
up works [59] also treat mechanics, including non-elastic effects, 
describing the redistribution of lithium within the metal foil due to 
inhomogeneous deposition. However, the mechanisms of swelling and 
stress driven diffusion and the effect of reduction of transport paths due 
to deformation are not considered. 

In our analysis, we first investigate trigger mechanisms in order to un- 
derstand possible reasons for inhomogeneous material deposition leading 
to damaging of the solid electrolyte. The calculations provide insight 


into the mechanisms dominating interface evolution during deposition 
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and are capable of predicting protrusion branching. The idea of protru- 
sion branching thereby stems from the current density distribution of 
the analysis of Section [6] which suggest elevated current density at the 
perimeter of a protrusion. We then put our focus on the morphological 
interface stability for PEO and PS-PEO electrolytes and compare the 
outcome with the results of Chapter 6] We will see that inhomogeneity 
of interface conductivity is a likely reason to cause protrusion growth and 
therefore we develop an analytical network model to depict a correlation 
of interface and bulk properties on morphological stability. We conclude 


this chapter with a recap of coupling mechanisms in protrusion growth. 


8.1 The Boundary Value Problem 


In this analysis, we model the negative metal electrode and the solid 
electrolyte and substitute the cathode by a current density and mass 
flux at the top side of the electrolyte, see Fig. Again, we assume 
an isothermal system at 80°C in two dimensions at plane strain. The 
electrode and the electrolyte are 20pm and 30pm thick, respectively. 
The width of the computation domain is set to 100 qm. 

A linear elastic material model is assumed for the lithium metal elec- 
trode, see Eq. (7.15). At elevated temperature, lithium has a Young’s 
modulus of approx 4 GPa [12] and the Poisson ratio v = 0.34. Its elec- 
tronic conductivity is k** = 10.6 SI [179]. The coupled system of ion 
transport and mechanics in the solid electrolyte is described by the 
extended concentration solution theory of Eq. (3.116) considering mi- 
gration, diffusion and stress-assisted diffusion of the ions. It takes also 
hyperelasticity and swelling of the host material into account. We again 
utilize the electrochemical parametrizations of PEO [56] and PS-PEO [57] 
with LiTFSi salt introduced in Section [4.2.1]and the interface properties 
discussed in Section i.e. a symmetric reaction with the exchange 
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Figure 8.1: The boundary value problem of the second stability analysis deals with a 
metal electrode at the bottom and the electrolyte on the top. 


current density given by Eq. with KBV = 0.589 Ss: The 
Poisson ratio is assumed to be v = 0.4 for both PEO and PS-PEO and 
the Young’s modulus of PEO and PS-PEO will be varied to gain more 
understanding of the system behavior. 

The boundary and transition conditions of the metal electrode and solid 
electrolyte are summarized in Fig. Symmetry boundary conditions 
are applied on the left and right side. Therefore we enforce continuity of 
the electrochemical properties {w+,c°} and zero horizontal displacement 
for {497 , pe", Cet}. It is assumed that the lithium metal is attached to 
a highly conductive and stiff current collector at the bottom interface 
and thus we have ® = 0 (electrically grounded) and ug = 0 (no 
deformation). The lithium electrode is attached to the solid electrolyte. 
Between the metal electrode and electrolyte, we assume compatibility in 
normal direction and allow tangential slip according to Eq. (7.8). The 
interface current density is computed by the Butler-Volmer equation 
(5.69) as a function of {wı,cı,®,0°"}, which determines the growth of 
the lithium metal interface, see Eq. (7.14), the electrical current at the 
electrode, and the current and mass flux of the electrolyte. The cathode 
is not explicitly modeled and thus we relate the mass and current flux 


to the applied current density J?PP. We further mimic a second metal 
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electrode at the top where lithium is homogeneously stripped with a 
boundary displacement proportional to J?PP, 

We simulate one charging step and deposit 4ym of lithium with 
JPP = 0.5 ms on the metal electrode. This leads to a charging time of 
5936s. Afterwards, we apply a relaxation phase of 3000s where J?PP is 


Zero. 


Again, we use COMSOL Multiphysics® to solve the coupled system of 
partial differential equations. We scale the field variables to improve con- 
ditioning of the system matrix and therefore increase numerical accuracy. 
The transition condition and are enforced with the concept 
of Lagrange multipliers [132], where scaling is done automatically by 
the simulation tool. We utilize a staggered solver and thus separate the 
transport equations from the deposition process. The first step solves for 
[w4 Cs yp? ®} and the second for {y*'”, Ce} including the Lagrange 
multiplier. The mechanics for the electrolyte is considered in both steps 
because it influences the ion-transport as well as the purely mechanical 
interaction with the electrode. If not otherwise stated, we use a struc- 
tured grid with quadrilateral elements, where the nodes on the interface 
between electrode and electrolyte align at time t = 0. Approximately one 
hundred thousand degrees of freedom were solved with the direct solver 
MUMPS in each time step. The backward differential method showed 
best performance for the numerical integration with respect to time. The 
integration order is automatically varied between one and five. The time 
steps are freely chosen by the simulation tool, but are restricted to a 
maximum time step At™®*. Numerical tests showed that At™** = 10s 


provides a good accuracy in a reasonable computing time. 
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Figure 8.2: Illustration of four hypothesis of trigger mechanism which might lead to 
inhomogeneous deposition. 


8.2 Analysis of Instability Triggers 


With the comprehensive description of electrode and electrolyte at hand, 
we now explore the growth of lithium metal during charging. Obviously, 
a planar interface as depicted in Fig. 8-1]in combination with homoge- 
neous transport properties yields homogeneous material deposition and 
stripping. To break the homogeneity, a trigger needs to be introduced. 
Fig. illustrates four hypotheses to be investigated in this context. 
Note that we will provide only a qualitative comparison based on edu- 
cated guesses of the parameter variations. A one-to-one comparison as 
done in non-dimensionalized settings, see e.g. [180], is not feasible in our 
case because variation of geometry, bulk properties, interface properties 
and boundary conditions are of different kind and cannot be easily related 


by one non-dimensionalized parameter. 


8.2.1 Four Hypothesis on Trigger Mechanisms 


Hypothesis 1 postulates an initial protrusion, e.g. due to surface rough- 
ness, as reason for a harmful protrusion growth into the separator. It 
influences the transport paths: The shorter the path, the higher the 
effective conductivity and the more likely an increased deposition pro- 
cess. In analogy to Chapter [6] we consider a Gaussian like protrusion, 
see Eq. (4.4), with an amplitude of AP = 4pm. We chose a support of 
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RP = 15m throughout this study. Other than in Chapter [6| we assumed 
the electrolyte to be stress free in the initial state. To resolve the curved 
shape, we utilize an unstructured grid with triangular mesh elements. 

Hypothesis 2 assumes an inhomogeneous ion supply at the cathode. 
Recall that the cathode is a porous composite with active storage particles. 
For example, particles close to the solid electrolyte that have lost their 
connection to the positive current collector are not functional, block 
the ion transport and thus increase the current density in the adjacent 
regions["] The morphology of the positive electrode is then mirrored on 
the negative electrode, due to local changes of current entering the solid 
electrolyte. We model this effect by 


(Jeppytrigger — (1 + w( X/ R38 )) FPP, (8.1) 


and thus double the applied current density at the center of the compu- 
tation domain. The support of the degraded zone is R@°® = 15 pm. 

Protrusion growth can also be triggered by inhomogeneities in trans- 
port properties, either in the bulk (Hypothesis 3) or at the interface 
(Hypothesis 4). Better transport properties in a region yield a higher 
deposition rate compared to the neighboring region and thus lead to 
instability. For Hypothesis 3, we increase the ionic conductivity and the 
diffusivity of the bulk with 


{r, D} 88" = (1 + 9w(X/R%8)) {r, D}. (8.2) 


by one order of magnitude in the center of the computational domain. 
Hypothesis 4 is similar to the adjustments of Eq. (6.5) and reads 


(KBV trigger = (1 4 2w(X/R%8)) KBV. (8.3) 


1 Active particles usually have a diameter of several micrometers which is 
the same order of magnitude as we use for the support in Eq. (8.1). 
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Figure 8.3: Four trigger mechanism for protrusion growth are investigated: (a) shows 
the protrusion height over time for the four hypotheses of Fig. (b) and (c) show 
the contour plot of hydrostatic pressure for hypotheses 1 and 3. The red lines indicate 
the initial interface shape. 


The exchange current density or, in other words, the interface conductivity 
is increased by a factor of three in the center of the interface and gradually 


declining to reference sidewards. 


8.2.2 Results and Discussion 


For this study, we use the solid polymer electrolyte PEO with a Young’s 
modulus of E = 10 MPa and transport properties as given in Section [4.2.1 
Recall that this material showed a high likelihood for morphological 
instability in the previous study of Chapter [6] 

Fig. [8.3h shows the growth of the protrusion height over time. We 
subtract the displacement measured at the side of the computational 
domain from total displacement of the electrolyte at the center in order 
to remove the translation of the interface due to homogeneous material 
deposition. 

The initial protrusion (Hypothesis 1) yields only a very small growth 
although the initial height of the protrusion is already significant. There- 


fore, one can conclude that an initial protrusion is not the main cause 
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of morphological instability. Nevertheless, we use this setup to under- 
stand material growth normal to the interface. The region close to the 
protrusion tip is under tension, see Fig. |8.3p, and the shoulder is under 
compression. This behavior is due to the broadening of the protrusion 
in the in-plane direction and results in quite strong compressive normal 
forces in the shoulder region and consequently pushes the separator up. 
In order to preserve compatibility, the perimeter as well as the tip have 
to pull on the electrolyte and are therefore under tension. The small pres- 
sure jump between electrode and electrolyte is due to the compatibility 
condition of Eq. (7.8) which requires only the normal traction force to 
align but not the lateral traction forces. Although the stress state with 
tension at the tip will always be present when a protrusion is growing, it 
is of second order when a planar interface is used as initial condition. 

Hypothesis 2 with an inhomogeneous current supply also shows only 
small protrusion growth. Ion transport in the in-plane direction smooths 
out the local current spike (J*??)'™88°", We note that this might change 
if either the separator thickness is reduced? or the transport properties 
of the polymer are degraded. 

Fig. [8.3f illustrates the interface for Hypothesis 3. Locally enhanced 
transport properties do not yield a pronounced protrusion growth. The 
modification of the interface conductivity (Hypothesis 4) results, however, 
in a significant increase of local deposition (Fig. |8.3h). With 4m the 
tip grows almost twice as fast as the neighboring regions. The growth 
speed even accelerates over time, meaning that the tip velocity increases 
the higher the protrusion is. This can be explained by the shortening of 
the transport paths. The effect of Hypothesis 1 and Hypothesis 4 then 
superpose. 

Fig. shows another interesting effect. There is a reorganization of 


lithium in the relaxation phase which yields to a decrease of the protrusion 


2 Note that Rosso et al. [35] measured protrusion growth for much thicker electrolytes. 
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height for Hypothesis 2-4 and a small increase for Hypothesis 1. Although 
no external current is applied, there is a small stripping current density 
at compression regions and a deposition current density at regions under 
tension. As discussed later, ions at the metal electrode prefer to be on 
a low pressure region. Since an ion as well as an electron conduction 
path is provided, there is no hurdle for a reorganization of the ions. 
The overpotential is, however, much smaller compared to doing charging 
which yields only a small decrease of protrusion height over time. 

The analysis of the trigger mechanisms suggests that Hypothesis 1-3 
are very unlikely to cause harmful protrusion growth, although the 
single influences were perturbed quite significantly. Hypothesis 4 with 
a perturbation of the interface conductivity has proven to be a likely 
candidate for instability. We do not claim that this is exclusive since 
purely chemical reactions might also contribute but are not depicted in 
our continuum description. Nevertheless, also experimental studies align 
with our outcome that the cause has to be related to the electrode-polymer 
interface. For example, Harry et al. [50] were able to correlate a protrusion 
growth to an impurity at the lithium surface. The impurity might increase 
the surface conductivity. Although this has not been experimentally 
proven, we will continue with Hypothesis 4 for the upcoming second 


stability analysis. 


8.3 Stability Analysis 


In the previous section, it was observed that an inhomogeneous interface 
conductivity, see Eq. (8.3), is the most likely trigger for protrusion growth. 
In the following study we stick to this trigger mechanism, vary the Young’s 
modulus and investigate the response of the resulting interface amplitude. 
We apply the analysis to both PEO and PS-PEO. 
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Figure 8.4: (a) The evolution of the protrusion height depending on the Young’s 
modulus for a PEO polymer. The contour plots on the right show the hydrostatic 
pressure in the electrode and electrolyte for (b) Æ = 10 MPa and (c) E = 300 MPa at 
t = 5900s. The red line indicates the initial surface at t = Os and the lines A and B 
are cross sections used to evaluate Fig. |8.7| Note that B is evaluated at the side of 
the computational domain. 


8.3.1 Lithium Metal and PEO Electrolyte 


We first analyze the PEO polymer electrolyte and vary the Young’s 
modulus between E = 10 MPa and E = 300 MPa. Fig. shows the 
protrusion growth over time. Note that the results for E = 10 MPa 
are equivalent to Hypothesis 4 in Fig. and recall that the protru- 
sion growth velocity increases with time due to the shortening of the 
transport path. Fig. shows the interface at t = 5900s. The tip 
experiences compression and the perimeter tension, similar to the study 
of Section [4.5] The hydrostatic pressure goes up to 3 MPa and is thus 
not very pronounced. If the Young’s modulus is increased, the effective 
protrusion height decreases significantly as seen in Fig. [8.4h. The rate of 
growth is identical in the beginning, but the more material is deposited, 
the slower is the protrusion growth for a higher Young’s modulus. The 
case of E = 300 MPa indicates an asymptotic behavior of the growth 
with a maximum protrusion height and therefore morphological stability. 


That means that even though the local interface conductivity is increased, 
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Figure 8.5: The contour plots show the displacement of the electrode (bottom) and 
solid electrolyte for (a) E =10 MPa, (b) E =100 MPa and (c) E =300 MPa. The red 
lines indicate the initial states of the interface. Note that the displacement of the 
electrolyte takes also translation due to deposition processes into account. 


the local stress state compensate this modification and yield to a homo- 
geneous deposition process. Again, we observe an amplitude decrease at 


the relaxation phase which is more pronounced the stiffer the material is. 


To understand the underlying effect of the decrease of protrusion height 
for increasing Young's moduli, we first consider the mechanical response. 
Fig. [8.5] illustrates the displacement of the electrode and the electrolyte. 
Note that the displacement of the solid electrolyte takes the translation 
due to deposition into account which explains the displacement of 3.2 ym 
for E = 10 MPa up to 3.6um for E = 300 MPa at the perimeter and 
t = 5900s. The difference between desired deposition (4m) and actual 
translation is explained by the protrusion growth. It requires additional 
material to create the protrusion which can not deposited onto the 
adjacent area. The displacement of the electrode is almost negligible 
for all three cases, especially if it is compared to the amplitudes of the 
protrusions. A simple push-down effect of the electrode based on elastic 
deformation to stabilize the interface can thus be excluded. 

We can explain the different protrusion heights with a change of the 
interface current densities. The current densities do not differ significantly 
for the three considered Young’s moduli at time t = 200s, see Fig. 
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Figure 8.6: The current density at the electrode-electrolyte interface at the beginning, in 
the middle and at the end of the charging process for (a) E =10 MPa, (b) E =100 MPa 
and (c) E =300 MPa. 


The Gaussian like shape of the current densities with a maximum of 
0.85 my comes from the modified exchange current density. When time 
passes and more material is deposited, we observe that for E = 10 MPa 
the current density increases with time, most likely due to the decrease 
of the related transport path. Contrary to that, the current density 
decreases for E = 100MPa and even more for E = 300 MPa as seen in 
Fig. and|8.6f. The smaller amplitudes of the current density explain 
the smaller protrusion growth for a higher Young’s modulus. 

Fig. shows the electrochemical potential w in the electrolyte 
and the electric potential ® in the electrode in the cross section along 
the Y-axis at X = 0 (line A in Fig. ) and X = 0.5W (line B in 
Fig. |8.42). 
stiffnesses thus we plot only the value for E = 10MPa. We observe 


The potentials at the side basically coincide for the three 


that ae is very similar for the three considered Young's moduli but 
differs strongly from aside, The electric potential is almost constant in 
the electrode. We note that the drop of electrochemical potential in the 
electrolyte is in the same order of magnitude compared to the potential 
jump along the interface with roughly 7 = 6mV. The concentration 
of the salt, again in the cross-section along the Y-axis at X = 0 (line 
A in Fig. B.4) and X = 0.5W (line B in Fig. [B.4F) at t = 1000s 
for the three different moduli, is shown in Fig. |8.7p. We observe a 
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Figure 8.7: The electro-chemo-mechanical response of PEO during charging at t = 
1000s for three Young's moduli: (a) Electric and electrochemical potential. Note 


that influence of E onto ® and wide is negligible and therefore not shown. (b) 


Concentration profile along cross section A (tip) and B (side), see Fig. . (c) Open 
circuit potential along the interface, here also for t = 5900s. 


concentration polarization in the out-of-plane direction along the Y-axis. 
At the perimeter, an increase of Young's modulus leads to a decrease of 
the concentration gradient. On the other hand, there is an increase of 
the concentration gradient or depletion of salt at the tip for an increasing 
Young’s modulus due to stress driven diffusion. For a further discussion 
of these effects we refer to Section [4.5] The open circuit potential ne: 
shown in Fig. for the three different Young’s moduli at t = 1000s 
and t = 5900s. The tip value of V$% decreases with respect to time 
and also with respect to the Young’s modulus. In this case, a negative 
open circuit potential provokes a decrease of the interface current density. 
This effect is especially significant for E = 300 MPa at t = 5900s with 
Vs% = —3mV and explains the reduction of interface current in Fig. 

In the relaxation phase at t > 5900s, we observe a decreasing protrusion 
height, see Fig. [8.4h. This is most pronounced for the solid electrolyte 
with E = 300 MPa. Fig. reveals the current density during relaxation 
at t = 6500s, t = 7500s and t = 8500s. There is a stripping current at 
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Figure 8.8: The figures show the relaxation process for PEO with E = 300 MPa for 
three different snapshots. (a): Interface current density, (b): open circuit potential 
and (c): hydrostatic mechanical pressure pM. 


the tip (¡BY > 0) and a deposition current (¡BY < 0) at the perimeter. 
This is triggered by the open circuit potential Via, see Fig. [8.8b, which is 
negative at the peak of the protrusion and positive at the side. Fig. [8.8 

shows the hydrostatic stress p™ for the three Young’s moduli along the 
electrode-electrolyte interface. There is compression at the tip yielding a 
negative value of V$4 (and thus stripping) and mechanical tension away 
from the protrusion (and thus deposition). The process is slow but will 
eventually yield to a completely flat surface in case if non-elastic effects 


are neglected. 


8.3.2 Lithium Metal and PS-PEO Electrolyte 


We redo the previous study, but now with the parametrization of a PS- 
PEO electrolyte, see Section 4.2.1] All other parameters are kept the 
same. In addition, we investigate a second system where we decrease 
the interface conductivity or exchange current density, respectively, by 
a factor of ten along the interface (KBV = 0.0589--A—). As we have 


Vm mol 
discussed in Section [6] the PS-PEO electrolyte requires higher stiffnesses 


for stability compared to PEO, therefore, we have chosen three Young’s 
moduli between E = 100 MPa and E = 1000 MPa. 
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Figure 8.9: The penetration depth over time for PS-PEO for three Young’s moduli 
and two interface conductivities. 


We observe a similar trend in the protrusion growth as observed for PEO 
in Fig. The higher the stiffness, the smaller the protrusion height. 
However, the protrusion height for the reference interface conductivity 
has values between 0.4 pm and 0.6 um and is therefore very small for all 
investigated Young’s moduli at the end of the charging process. Numerical 
studies showed that moderate changes of charging rate J?PP and the 
amplitude of the interface conductivity disturbance in Eq. (8.3), did 
not increase the instability significantly. By adjusting the prefactor of 
the interface conductivity KPY, we again observe a strong impact of the 
Young’s modulus on the protrusion height. Degrading the surface by 
dividing KBV by a factor of ten and therefore increasing the interface 
resistance by one order of magnitude yields a significant protrusion growth 
with values between 2.3 pm and 3.2 um, see Fig. In both cases, the 
influence of the mechanical stiffness is not as pronounced as for the PEO 
electrolyte. 

Fig. [8.10] reveals a detailed analysis of the electrochemical potential 
for E = 100 MPa at t = 1000s for the two prefactors of the interface 


conductivity KPV and 0.1K®Y. The graphs in Fig. and [8.10N are 
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Figure 8.10: The electrochemical potential for PS-PEO with E = 100 MPa at t = 1000s 
for two interface conductivities. (a) Contour plot of w+ for reference KBV. (b) Contour 
plot of w+ for 0.1KPV. The lines A-D indicate cross sections. Note that B is evaluated 
at the side of the computational domain. (c) Cross section in the direction A and 
B. (d) Cross section in the direction C and D. Note that the figures (a) and (b) are 
truncated to highlight the region around the protrusion growth. 


defined by the cross sections A, B, C and D, as depicted in Fig. [8.10p 
and Fig. [8.10p. Recall that the gradient of the electrochemical potential 
w4 is a major driving force for ion flux and is therefore a good indicator 
for the current in the electrolyte. 

In Fig. we see that the electrochemical potentials in the out-of- 
plane direction are shifted by roughly 50 mV if the interface conductivity 
is smaller but the slopes are similar. On the other hand, there are strong 


differences in the variation of w+ in the in-plane direction, see Fig. [8.10H1. 


For the standard interface conductivity, only a very small drop of w+ 


is observed in the center region (cross section C). Further away from 
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the electrode-electrolyte interface (cross section D) we observe an almost 
constant wy in the in-plane direction. This is different for the second 
case with 0.1K ®Y, where cross sections C and D show a notable drop of 


w4 in the region where the interface conductivity is increased. 


8.3.3 Discussion 


The results of protrusion growth at a metal electrode adjacent to a 
PEO and a PS-PEO electrolyte, respectively, extend the first study of 
Chapter [6] and give new insights into the morphological stability. Con- 
trary to the first study with a given deformation, the protrusion growth 
now has to be trigged with the help of an inhomogeneous interface con- 
ductivity. The two approaches are therefore not directly comparable, 
although the study of Section [6.3.5] with its modification of the interface 
conductivity is based on similar assumptions. 

We first point out that the shape and size of the protrusion chosen in 
Section [6] align with the outcome of the dynamic protrusion growth from 
this study. A bifurcation of the protrusion as indicated by the current 
density distribution of Fig. with the highest deposition rate at the 
perimeter to each side of the protrusion is not observed. This confirms 
the idea that the build-up of secondary tips is suppressed due to the 
electro-chemo-mechanical coupling and aligns with experiments which 
have shown no morphology change of the protrusion during charging 

The PEO material obeys a strong dependency of the growth speed on 
mechanical stiffness. This is in agreement with the results of Chapter [6] 
(Fig. [6.8p), where we have seen a characteristic change in stability 
at approximately E = 100 MPa. If the Young’s modulus is smaller 
than E = 100 MPa, then the growth accelerates, if it is larger, then an 


3 Note that repeated stripping and deposition can yield a porous structure in polymer 
electrolytes [50]. 
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asymptotic limit will be reached eventually. Both the modification of the 
interface conductivity and the shortening of the transport path will then 
be balanced by the mechanical stress at the protrusion peak. 

The deposition current has a maximum at the peak of the protrusion 
for all calculations and is hence different to the study of Chapter [6] 
which has shown a strong drop of the current density, partially even a 
stripping current, at the tip for high Young’s moduli (see Fig. |6.3). The 
mechanical stresses in the soft electrolyte (E = 10 MPa) are small and not 
able to lower the effective current density at the tip although the stress 
increases with protrusion depth. In fact, when the protrusion grows, the 
current densities increases at the peak of the protrusion, see Fig. [S-Gh, 
which eventually leads to an acceleration of growth as seen e.g. in Fig. 
for PEO and Fig. [8.9] for PS-PEO. If a Young's modulus above a 
critical value is chosen, then the variation of current density decreases 
while the protrusion is growing. This indicates a maximum height of the 
protrusion. The deposition rate at the tip equals the deposition rate at 
the flatter regions in this state. The protrusion growth is then stopped 
and a morphological stability is reached. 

PS-PEO has this pattern at around E = 500 MPa, which is slightly 
higher than the outcome of Chapter [6] However, the overall growth is 
small without a modification of KBV, even for small Young’s moduli, 
and a failure of the battery is not likely. This behavior changes when the 
electrochemical interface property such as KBV is degraded. PS-PEO 
then shows a similar amplitude of the protrusion as the PEO electrolyte 
without interface degradation. This indicates that a decrease of bulk 
transport and an increase of interface properties decrease the likelihood of 
protrusion growth. The correlation of bulk and interface properties with 


respect to morphological stability will be analyzed in the next section. 
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Figure 8.11: Problem description of the network model: (a) Axis symmetric section of 
an undeformed electrolyte-electrode system. The red area denotes a degraded interface 
between the electrode and electrolyte. (b) Sketch of the electrochemical potential 
w+ (isolines) and the current density in an electrolyte with degraded interfaces. (c) 
Idealized setting of ion transport in an electrolyte within a network model. The ions 
from the top are either transported directly through the bulk (vertical) or in-plane 
(horizontal) to the electrode. 


8.4 An Analytical Assessment of 
Morphological Stability 


We now want to rationalize the findings of the previous simulations such 
as the correlation of bulk and interface properties with the help of an 
analytical network model. We restrict the analysis to an undeformed 
geometry as shown in Fig. where the exchange current density is 
degraded locally at the electrode-electrolyte interface (red area, see also 
Hypothesis 4 in Section 8.2.1). The undeformed geometry corresponds 
to the initial state of the deposition process in Section [8.3] In this state, 
mechanical effects such as stress driven diffusion or change of length of 
transport paths due to deformation are negligible and will not be taken 


into account. 
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Fig. illustrates isolines of the electrochemical potential in the 
solid electrolyte as observed previously in Fig. |8.10p. The mass flux 
follows the gradient of w and is therefore perpendicular to the isolines. 
It is obvious that the in-plane contribution of the mass flux is directly 
related to the protrusion growth. The higher the current density in the 
in-plane direction, the more ions are transported into the degraded area 
and are thus deposited. This leads to the inhomogeneous growth and 
eventually to morphological instability. In the following, we approximate 
the continuous model of ion transport by a network model, separate the 
mass flux into a through-plane and an in-plane contribution in the region 
of the degraded area, see Fig. |8.11£, and deduce an expression for the 


harmful in-plane mass flux. 


8.4.1 Network Model of lon Transport 


To derive an analytical expression for the in-plane contribution of the 
ion flux, we consider an electrolyte with a locally degraded interface to 
the metal electrode. The domain of the electrolyte above the degraded 
zone is denoted as Bo|*° and the surrounding domain as Bo|"*', see 
Fig. [8-11f. A constant current density j®PP is applied at the top (cathode 
side) and the current densities at the interface are denoted as ¡BV|des 
and jPY |" for the degraded and non-degraded areas, respectively. Due 
to the degradation, we postulate a current density in in-plane direction 
Aj which is defined to be positive in the direction from Bo|'* to Bo|*°. 
Using the network model of Fig. [8.1 1f, we can now apply conservation 


of mass in the nodes and get 


¿3 ¡des Ajdes = joer Ajdee aL AAj 


‘ (8.4) 
Ar = jar ppret un AAj 
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with the degraded interface area A|d®, the remaining interface area A|"% 
and the area of the cylindrical surface of Bp|98 as A. Eliminating the 
in-plane current AAj in Eq. (8.4) yields 


ae Ey is = ¡PP(1 + a) (8.5) 
where a = AZ = AU _ı To obtain expressions for ¡BY|98 and 
= Altes — Aldes . P J 


j®V ref, we utilize a linearized version of the Butler-Volmer equation 
such as Eq. (5.66). The degraded interface has an exchange current 
density jo|@°®. Otherwise, the exchange current density is jo|"* which 
is proportional to K®Y, see Eq. (5.68). The electric potential of the 
electrode is given by ® = 0 (grounded). Furthermore, we consider a stress 


free state and thus have V{4 = 0. Eq. (5.69) then yields 1/48 = —w, |des 


and n= —w; |" and the two linearized Butler-Volmer equations are 
pPY [28 tee, ee 
Fal = —bjo| fw |£, 


with b = de They give a linear relationship between current densities 
and electrochemical potentials w}. Combining Eq. (8.5) and Eq. (8.6) 
yields 


l+a 


Joe |18 + ajo [Em pr = —79PP 5 


l (8.7) 


Subtraction of jo[YEw, |" on both sides then gives 


a : ann + Q 
Jo Elot — 49, pu) = er 


= (joldee a ajol**) e (8.8) 
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In the next step we combine Eq. (8.4), Eq. (8.8) and the linearized 
Butler-Volmer equation and obtain 


1 (a A ag 
w4 |des = oer = bolt © PP_-(1+(1-x)a) faas) (8.9) 


jol! 
Joldes 


example, if the interface conductivity in the degraded zone is higher 


with x= 1 — as a measure for the interface inhomogeneity. For 
compared to the side, then x > 0. In the next step, we assume a 
linear potential drop in the in-plane direction of the electrolyte bulk and 


compute the in-plane current density as 


[est _ w4 |388 
FAL 


Aj = „et (8.10) 


where x“ is the bulk conductivity and AL is a characteristic length 
between the Bo|" and Bo|*°8. Inserting Eq. in Eq. (8.10) yields 


eff A 
K A 
Aj = — 5 PP — (1+(1- ——Aj|. 8.11 
j= FF ç (1+ (1x) 4) Fa ) (8.11) 
Solving the implicit equation for Aj then yields 
` X f 
Age E (8.12) 
EALE + aay + Ae 


and thus gives the in-plane current Aj as a function of geometric proper- 


ties, the interface conductivity jo|"*, the bulk conductivity «° and the 


jor 


o [dez Assuming a circular degraded zone of radius Rs and 
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i el Â _ 2mRädeLey _ 2L% 
a hight of the electrolyte of L°Y one can deduce aje = (pas? = Ree 


and together with AL = dR@“ with d > 1 we finally obtain 


. X . 
Aj = —— - APP (8.13) 
Te Lely (1— 
fo Rod + ata + EY 


as a second approximation of the in-plane current density. 


8.4.2 Discussion 


Eq. (8.12) is a first order approximation of the in-plane current density 
which describes the material transport into the degraded zone. In general, 
morphological instability is found for x > 0 which is equivalent to 
jols > jo". If the exchange current density is homogeneous, then 
a stable deposition is expected. If the degraded interface has a much 
higher interface conductivity, i.e. x — 1, then protrusion grow is strongly 
promoted. In general, the higher the difference between the exchange 
current densities, the larger the protrusion growth. The applied current 
density j*PP enters linearly and indicates higher morphological instability 
issues for higher current densities. The linear correlation aligns with the 
findings of the numerical study in Section [6.3.4] 

Eq. (8.12) indicates a smaller in-plane flux for decreasing bulk conduc- 
tivity «° which is thus beneficial for morphological stability. Qualita- 
tively speaking, a small bulk conductivity penalizes in-plane ion transport. 
Note that this finding does not contradict the result of Section 
where an electrolyte with higher bulk conductivity showed lower require- 
ments for mechanical stiffness (see e.g. Fig. (6.6). The study in Section 
[6.3.3] correlates bulk transport properties with mechanical properties on 
a deformed geometry, whereas the analysis of Eq. deals with the 
initiation process. In other words, Eq. determines the deposition 
growth at time t = 0 (e.g. the initial slope in Fig. and not the 


growing pattern when a protrusion builds up (t > 0). However, both are 
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not decoupled as discussed in Section [6.3.5] (Fig. [6.8b), where a change of 
interface conductivity at the degraded interface also changes the require- 
ments on mechanical stiffness. The analytical model and the simulations 
agree that a degraded zone with jp[98 > jol* promote instability. Fur- 
thermore, the network models aligns with the numerical studies with the 
prediction that electrolytes with poorer electrochemical transport prop- 
erties such as PS-PEO are more robust against inhomogeneous interface 
conductivities. 

Another way to stabilize the system is to increase the reference interface 
conductivity jo". In this case, the drop in electrochemical potential 
along the interface is smaller and thus the driving force for the in-plane 


- ¡ref 
current decreases. The ratio D T that is found in the denominator of 


Eq. (8.12) aligns with the findings of Section In this study, we have 
observed almost the same protrusion growth for a system with PS-PEO 


(and an exchange current density decreased by a factor of ten) and for a 
system with PEO, which has roughly ten times better transport properties 
(see. Fig. |4.2.1). The analytical expression confirms therefore 
that interface properties and bulk properties cannot be decoupled in a 
morphological stability analysis. 

We further note that Eq. indicates a decrease of Aj for thick 
electrolytes. This effect, although not very pronounced, was observed in 
the parameter study of Section[6.3.3]for h > 50 pm. For small thicknesses, 
the mechanical effects dominate. We further note that the in-plane 
current density vanishes for a small degraded area (R*° — 0) and a 


large degraded area (Rd — 00). 


8.5 Recap of the Coupling Mechanisms 


The previous studies have shown several mechanisms which contribute 


to protrusion growth of a metal electrode in a solid state battery. In 
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what follows, the driving forces for protrusion growth are summarized 
and discussed. 

In Section [6.2] we have seen that the overpotential 7 at the SEI is a key 
ingredient in understanding the difference in deposition rates along an 
interface. Recall from Eq. and Eq. (5.61) that the overpotential 
consists of three components: the electrode potential ®, the electrolyte 
potential w, and the open circuit potential Voa. T he electrode potential, 
however, can be assumed to be constant due to the high electronic 
conductivity in the metal electrode. We assume that the electrode is 


grounded and thus set to ® = 0. The overpotential then reads 
n = —(w + VER). (8.14) 


Recall that a negative overpotential yields a deposition current. We 
emphasize that both the exchange current density jg and the overpotential 
determine the current density. In the following, we neglect the effect of 
the mechanical dependency of the exchange current density as we have 
seen in Section [6.2] that the overpotential is the dominant mechanism for 
inhomogeneous deposition. As shown in Fig. and Fig. [S-10F, the 
electrochemical potential is almost constant at the counter electrode (Y = 
L*y) although a protrusion growth is triggered. Therefore, we assume 
that the electrochemical potential at the separator-cathode interface 
w+l*P remains equal during protrusion growth. This also implies a 
constant current density j?PP|'P. 

Fig. [8.12] summarizes the effects which influence the overpotential and 
shall be discussed in the following. We consider the electrochemical and 
mechanical contributions on both w+ and V4. Quantitative statements 
with absolute values of the single effects would require knowledge of 
the material system, mechanical state, charging rate, etc. Therefore, 
this discussion shall serve as a qualitative study to gain an overview of 


the different mechanism and of which kind (stabilizing, destabilizing) 
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Figure 8.12: Processes that changes the electrochemical potential w+ and open circuit 
potential vor at the electrolyte-electrode interface clustered in three categories: 
Standard response (no inhomogeneity) in form of bulk conductivity (S1), concentration 
polarization in the bulk (S2) and at the interface (S3); Trigger mechanisms due to 
locally increased interface properties at area A4°8 (T1) and in-plane current flux (T2); 
Protrusion growth yields pressure in the electrolyte (P1), stress driven diffusion (P2), 
change of transport paths (P3) and normal stress at the interface (P4). The color 
code of the arrows is as follows: Black (initial response when current is applied), 
purple (changes due to concentration polarization), red (destabilizing process), green 
(stabilizing process) and blue (overpotential of each condition). 


they are. We cluster the effects in standard processes (S) specified by 
homogeneous ion transport, intermediate effects due to a trigger (T) and 


effects while protrusion growth (P). 


198 


8.5 Recap of the Coupling Mechanisms 


We start with the standard processes due to ion-transport in a homoge- 


neous electrolyte: 


e S1: At the initial state after applying a current, there is a drop 
in electrochemical potential in the electrolyte due to Ohm’s law 


inversely proportional to the ionic conductivity kK. 


e S2: A concentration polarization occurs for a galvanostatic current 
as discussed in Section [4.4] The ionic conductivity then decreases 
(see discussion in Section |3.5.4), which increases the drop in elec- 


trochemical potential in the bulk. 


e S3: The change of concentration also has an impact on swelling 
and osmotic pressure which eventually lead to mechanical stress. 


The open circuit potential —V¿2|s is therefore not zero anymore. 


Note that these steps occur regardless whether the interface is degraded or 
not and are therefore not relevant to the considerations of morphological 
stability. This state now defines the overpotential n|$d at a flat, not 
degraded interface. 

As a trigger mechanism, we follow the previous findings and degrade 


the interface A4°® with an increase of interface conductivity: 


e T1: A lower electrochemical potential is necessary to transfer the 


same amount of current (see discussion in Section 18.4). 


e T2: The drop in w+ leads to a change of the ion transport path, see 
Fig. [8.10h. Following the highest decent, more ions will enter the 
area of Ads, which ultimately leads to an increase of the current 
density in this region, see Eq. (8.12). The increase of current 
density results, however, in an increase of the overpotential, see e.g. 

` . deg . à 
Eq. (5.66). We denote this overpotential as nlg., if the interface is 
still flat. 
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The increased current density at the degraded zone compared to side 
regions then leads to a protrusion growth, which affects the overpotential 


as follows: 


e Pl: The protrusion growth yields mechanical stresses in the solid 
electrolyte. The hydrostatic pressure pM increases w, due to Eq. 
(3.77) and influences the overpotential such that more deposition 
takes place at Ads, 


e P2: The hydrostatic pressure also has the effect that the ions are 
depleted in the region around the protrusion tip (see Section 4.5), 
which decreases wą and therefore stabilizes the interface current 


density. 


e P3: If the protrusion grows, the effective bulk conductivity increases 
due to a shorter transport path, which then promotes instability 


(see e.g. Fig. |8.8h). 


e P4: Similar to the hydrostatic pressure of the electrolyte, we obtain 
positive normal stresses at the metal electrode. This lowers the open 
circuit potential Ve; decreases the overpotential and is therefore a 


stabilizing effect. 


Ultimately, if the material properties of a solid electrolyte are chosen 
such that n]@°® at time t decreases with protrusion height, than the 
deposition current will also decrease and morphological stability might 
be reached. Ideally, only a small or even infinitesimal protrusion is 
required to equilibrate stabilizing and non-stabilizing effects. That might 
be the case for ceramic electrolytes, which are single-ion conductors. 
They do not feature the effects Pl and P2 due to the missing swelling 
mechanism. The stabilizing effect of P4 is then dominating even for small 
deformations[] We refer to McMeeking et al. [91] and Klinsmann et al. 


4 Ceramic electrolytes struggle, however, with lithium growth in grain boundaries. 
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for a more detailed study of stability in ceramic electrolytes. The 
coupling mechanism in polymer electrolytes is, however, more complicated 
and showed to be very sensitive to small changes in mechanical stiffness 
as discussed in Section [6] Deducing morphological stability maps for 
these systems require therefore a comprehensive electro-chemo-mechanical 


formulation as presented in this thesis. 


8.6 Summary 


With the help of the comprehensive model of coupled ion transport, 
reaction kinetics, deposition and electrode response, we were able to 
model the growth of a protrusion in a binary solid electrolyte. We 
investigated within our continuum mechanical framework four different 
trigger mechanism which can cause harmful protrusion growth. The 
analysis indicates that a heterogeneous interface conductivity (modeled 
by a local increase of the exchange current density) is most likely to 
initiate a protrusion growth with a deposition rate big enough to penetrate 
the solid electrolyte. We studied the impact of the Young’s modulus 
on the protrusion growth for the two parameterizations of the polymer 
electrolytes PEO and PS-PEO. Both showed a characteristic change 
of protrusion growth at E = 100 MPa and E = 500 MPa, respectively, 
and align quite well with the results of Section [6.3.5] Below this value, 
we have seen an acceleration of growth whereas an asymptotic limit 
is indicated for higher Young’s moduli. We have confirmed that the 
suppression mechanisms do not stem from a deformation of lithium but 
from a decrease of the interface current density for higher Young’s moduli. 
Furthermore, the model indicates a self healing effect due to residual 
stresses after penetration. The influence of the mechanical stress on the 


open-circuit potential of the electrode-electrolyte interface leads to a 
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stripping current at the protrusion tip and deposition at the perimeter. 
This effect is more pronounced the higher the Young’s modulus is. 

The simulation offered insights into the correlation of bulk transport 
properties and interface properties. With the reference interface conduc- 
tivity used in this work, we observed only a small protrusion growth for 
PS-PEO. With a decrease of the interface conductivity, we could trigger 
a harmful protrusion growth again. This response could be confirmed by 
an analytical network model, which describes a required in-plane current 
density for protrusion growth as a function of inhomogeneity, geometry, 
applied current, bulk conductivity and exchange current density. Accord- 
ing to this model, a low applied current, a low bulk conductivity, a high 
exchange current and a homogeneous exchange current density are all 
stabilizing. This aligns with the findings of Section |6.3.5} which states 
that PS-PEO is more robust against interface imperfections compared to 
PEO and of Section [6.3.4] which indicates a linear correlation of applied 
current and mechanical stiffness. Finally, we summarized the influence 
of the different coupling mechanism in a binary solid electrolyte on the 


morphological stability. 
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The scope of this work was the development of a model for a polymer 
based solid state battery with a metal electrode in order to get insights 
into the morphological interface stability that is a major concern in this 


kind of next generation energy storage system. We presented therefore 


e a transport theory tailored for binary solid electrolytes, 
e a theory for reaction kinetics including mechanical effects, 
e a description of deposition and stripping for metal electrodes. 


The novel transport model combines the idea of finite strain continuum 
mechanics and the widely utilized electrochemical description of batter- 
ies, namely the concentrated solution theory. Herein, both anions and 
cations are mobile and, during charge and discharge, form a concentra- 
tion gradient along the electrolyte. Besides migration and diffusion, we 
considered stress driven diffusion in the binary solid electrolyte, a direct 
consequence of swelling of the host material. Using the electroneutral 
limit we identified the required input parameter by those accessible by 
experiments. The electro-chemo-mechanical coupling mechanism were 
discussed using several exemplary simulations. 

Another important aspect of this work is the assessment of the influence 
of mechanics on the reaction kinetics between a solid electrolyte and a 
metal electrode. Following the concept of transition state theory, we 
elaborated the influence of mechanical stress on the electrochemical equi- 
librium (Nernst equation) and non-equilibrium (Butler-Volmer equation). 
The generic concept using energetic states of electrode, electrolyte and 
energy barrier is able to describe a broad range of interface combinations 


consistent with literature. We put a special focus on a metal electrode 
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where the normal stress, not the hydrostatic stress, was shown to be the 
variable of interest for the thermodynamics of deposition and stripping. 

The third aspect required to model a solid state battery is a consistent 
description of the metal electrode. Cations are deposited on and stripped 
from the metal electrode and therefore we introduced (in addition to the 
deformation map) a second mapping describing the stress free change of 
material shape. With the formulation of a compatible interface between 
solid electrolyte and metal electrode, we have gathered all elements to 


describe a solid state battery. 


In this work, we have carefully analyzed the morphological stability of a 
binary solid state battery to determine the impact of transport properties 
and mechanical stiffness on the likelihood of undesired protrusion growth. 


We were able to conclude that 


e an electrolyte with higher mechanical stiffness is in general more 
prone to morphological stability. 

e at the same time, an electrolyte with better bulk transport proper- 
ties has less requirement on the Young’s modulus. 

e the influence of mechanical stress on electrochemistry causes stabil- 
ity and not, as widely believed, a purely mechanical push down. 

e instabilities are most likely caused by inhomogeneous interface 
properties such as the interface conductivity. 

e materials with better bulk transport properties are more sensitive 


to inhomogeneous interface properties. 


We executed two studies on morphological interface stability. The first 
assessment assumed an infinitely stiff electrode penetrating the solid 
electrolyte. Varying the Young’s modulus led to a characteristic change 
of the current density along the deformed interface. Above a critical 
Young’s modulus, the current density at the tip was below average which 
we interpreted as stable. Several parameter studies showed that an 


increase of transport parameters (by means of diffusivity, conductivity or 
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transference number) lowers the requirement on the mechanical stiffness. 
The predictions of the required Young’s moduli for PEO and PS-PEO 
electrolytes align with experimental findings. 

With the second study on morphological interface stability, we explic- 
itly modeled the material deposition during charging. We found that 
protrusion growth is most likely caused by an inhomogeneous interface 
conductivity. Again, we observed a lower deposition rate for stiffer elec- 
trolytes. This is caused by a decrease of the interface current density 
rather than a push down effect of the metal electrode. There is a tran- 
sition of accelerated to decelerated growth at approximately the same 
Young’s modulus as predicted by the first study. Furthermore, we ratio- 
nalized findings of the numerical simulations with an analytical network 
model namely that solid electrolytes with poor transport properties are 
more robust against inhomogeneous interface conductivity. 

In conclusion, we have seen that a polymer based electrolyte with a 
certain balance of stiffness and electrochemical transport properties in 
combination with a sufficiently homogeneous interface conductivity is a 


promising concept to enable stable lithium metal electrodes. 


The presented model of a solid state battery can be extended by various 
aspects. One should overcome the assumption of pure elastic response 
and use an elastic-visoplastic material model as recently parametrized for 
lithium metal [178]. In the same course, an adequate material model for 
the solid electrolyte going beyond an elastic stability analysis is needed. 
Further, anisotropy of electrolytes such as PS-PEO with a laminar struc- 
ture [151157] and the metal electrodes [14] should be examined. Depending 
on the charge-discharge profiles, different protrusion morphologies were 
observed by Harry et al. [50] and should be investigated further. 
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A Nomenclature 


Kinematics 


[material frame; spatial frame; unit; first appearance; description| 


un 


un 


“|B «|B B 


mR Ble Bie B 


A a = 


time 


time domain 


mapping from material to spatial frame 


host velocity 

velocity of species a 

spatial coordinate 

gradient operator 

divergence operator 

normal vector 

tangent vector 

volume element 

volume element of the host 

area element 

control volume w.r.t. to the host 
control volume w.r.t. species a 
surface of control volume Po 
body of the host material 
displacement vector 

deformation gradient 

rotation and right stretch tensor 
Jacobian of deformation gradient 
elastic and swelling part of F 


elastic and swelling part of J 
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A Nomenclature 


Kinematics 


[material frame; spatial frame; unit; first appearance; description] 


Y y Eq. scalar or vector field 

Y Y Eq. vector field defined on a surface 

Y gy Eq. scalar quantity defined on a surface 
I 1 identity matrix 

w 1 Eq. Wendland function 


Kinematics of Deposition and Stripping 


[symbol; unit; first appearance; description] 


cet m mapping from geometry to material frame 

Xy m geometric coordinate 

d m geometry displacement field 

Boo body describing geometry change but no deformation 
St deformed interface between electrode and electrolyte 


Fields [Mechanics] 


[material frame; spatial frame; unit; first appearance; description] 


BM pM => Eq. mechanical body forces 

BE bE = Eq. electrical body forces 

BS bf “Ss Eq. chemical body forces 

BY © Eq. body force w/r to component inertia 
Ba => Page extended specific body force 

TM +% == Eq. mechanical traction force 

TE tE == Eq. electrical traction force 

Te t© Y, Eq chemical traction force 


B 
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A Nomenclature 


Fields [Mechanics] 


[material frame; spatial frame; 


unit; first appearance; description] 


P 

pM 
pE 
pc 


I, I2, Is 
c p Je 
a eo 
AS 
e 


<b 
A las ae jaja jefa jaja jo 


Eq. 
Eq. 


Fields [Electrochemistry] 


w 


jan 


= 
a 
© 


= 
y 
=) 


p 
y 
2 


A 
w 
© 


ot 
a 
N 


total Piola & Cauchy stress tensor 


mechanical Piola stress tensor 


Maxwell stress tensor 


stress tensor due to chemical forces 


normal compressive Cauchy stress 


deviatoric Cauchy stress tensor 
total hydrostatic pressure 
pressure due to mechanical strain 
osmotic pressure 

angular momentum 

total torque 

invariants of F 

prefactors defined in Eq. 
parameter from swelling model 


absolute temperature 


[material frame; spatial frame; unit; first appearance; description] 


Na 

N 
ce Ca 
P& Pa 


o 
Prot Ptot 


Pda Pla 
Pg Pq 
Ja da 


5 pal 
e 


al 


gB 
o, 


[IS 


3 


amount of atoms of species œ 
number of species 
concentration of species a 
partial mass density 

total mass density 

partial charge density 

total charge density 


total molar flux of species a 
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A Nomenclature 


Fields [Electrochemistry] 


[material frame; spatial frame; unit; first appearance; description] 


J aa 4 charge flux of species a 
JEJ e at electric and chemical part of Ja 
Jtot Itot == total mass flux 
Jq ia A, total ionic current density 
J E J T A, electric and chemical part of Jq 
(J “ye 4 steady state current density 
Ra Ta ze concentration sink or source 
Rag Tae AS charge source 
D d as electric displacement 
e Y electric field 
P P a electric polarization 
WSS As surface charge 
RH r = E heat source 
JE j! -z heat flux 
Eint, grin J interal and kinetic energy 
IM. DP LE 2 power supply due to mechanics, heat 
and electrostatic 
ea = specific internal energy of species a 
w 5 total Helmholtz energy 
Vo = specific Helmholtz energy 
a YE, pe = mech., elec. and chem. part of db 
Na ee specific entropy of species a 
ôy un dissipation density 
Ho nn chemical potential of species a 
us, pM mn chemical and mechanical part of pa 
® V electric potential 
Wa E electrochemical potential of species a 
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A Nomenclature 


Interface Kinetics 


[symbol; unit; first appearance; description] 


et el 
ke k kely 


Jo 


sia 


et 


J> 


>E 
B 


3 
e 


|- <q < <Á EJ EJ EJ > 


vm mol 


current density by Butler-Volmer kinetics 
rate constant 

exchange current density 

prefactor in exchange current density 
reference value of exchange current density 
anodic and cathodic current 

open circuit potential w.r.t. ® 

open circuit potential w.r.t. w+ 
overpotential 

molar Gibbs energy 

energy barrier 

symmetry coefficient 

prefactor to relate barrier and bulk energies 
molar Gibbs energy incl. electric potential 
success frequency of forward /reverse reaction 
area density 

energy due to deposition and stripping 
symmetry coefficient (mechanics) 

interface resistance 


prefactor for metal electrodes 
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A Nomenclature 


Parameter 


{material frame; spatial frame; 


unit; first appearance; description] 


Ma 
Qa 
Qe 
Bs 
Za 
ez Er 
cret 
Me, Mag 
M?,M?  My,M_ 
Mon, 
Kap, Kag 
fa 
OS 
K? K 
(k?) (1555 
Ke 
to ta 
D* D 
Ne 
NH) ANH 
E,G 
V 
AP 
RP 
A 
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3 
a 


N 


mJ 


mJ 


Nal vn 


molar mass of species a 
molar volume of species a 


molar volume of species a (de- 


formed) 

factor to account for non- 
ideality in swelling 

charge number 

relative electric permittivity 
reference concentration 


mobility matrix between 


species a and 8 

mobility of cation and anion 
cross mobility 

ambipolar mobility matrix 
drag coefficient 

fugacity or activity coefficient 
thermodynamic factor 

ionic conductivity 

steady state conductivity 
mass conductivity 
transference number 
diffusion coefficient 
Newman number 

lamé constants 

Young’s and shear modulus 
Poisson ratio 

amplitude of protrusion 
radius of protrusion 


amplitude of heterogeneity 


A Nomenclature 


Macroscopic Properties 


[symbol; unit; first appearance; description] 


F N Eq. 
Wı,W2 J Eq. 
AE J Eq. 
A m? Eq. 
MLi kg Eq. 


Network Model 


force 

potential energy (total) 
electrochemical energy (total) 
electrode-electrolyte interface area 


mass of the electrode 


[symbol; unit; first appearance; description] 


Aj + Eq. 
Alte? m? Eq. 
a 1 Eq. 
1 

b y Eq. 
x 1 Eq. 
d 1 Eq. 
Constants 


in-plane current density 
cylindrical surface of degraded zone 


ratio of reference and degraded area 


prefactor 


prefactor 


measure for interface resistance inhomogeneity 


[symbol; value; unit; first appearance; description] 


F 964885 As Faraday constant 

€0 8.854 - 10-12 E electric permittivity of vacuum 

kg 1.38- 10728 x Boltzmann constant 

h 6.626-10-34 Js Planck constant 

R 8.314 ki universal gas constant 

Indices 

etr electrode iso isotropic 

ely electrolyte ref reference 

M metal eff effective 

+ Cation PEO poly(ethylene oxide) 
Anion PS-PEO block copolymer 

BIC Binary conducting BV Butler-Volmer 

ME Metal electrode app applied 

AP Active particle dep deposition 
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